NASA 

Contractor Report 4342 


AVSCOM 
Technical Report 90-C-028 


Local Synthesis and 
Tooth Contact Analysis 
of Face-Milled Spiral 
Bevel Gears 


Faydor L. Litvin and Yi Zhang 


GRANT NAG3-964 
JANUARY 1991 


r (MASA-CR-4342) LOCAL SYNTHESIS AND TOOTH N91-19439 

CONTACT ANALYSIS OF FACE-MILLED SPIRAL BEVEL 
GEARS Final Report {Illinois Univ.) 186 p 

CSCL 131 Unclas 

Hl/37 0332290 


NASA 



AVIATION R&T ACTIVITY 



NASA . AVSCOM 

Contractor Report 4342 Technical Report 90-C-028 


Local Synthesis and 
Tooth Contact Analysis 
of Face-Milled Spiral 
Bevel Gears 


Fay dor L. Litvin and Yi Zhang 
University of Illinois at Chicago 
Chicago, Illinois 


Prepared for 
Propulsion Directorate 
USAARTA-AVSCOM and 
NASA Lewis Research Center 
under Grant NAG3-964 


NASA 

National Aeronautics and 
Space Administration 

Office of Management 

Scientific and Technical 
Information Division 


1991 




TABLE OF CONTENTS 


SECTION PAGE 


1 LOCAL SYNTHESIS of GEARS(GENERAL CONCEPT) 1 

1.1 Introduction 1 

1.2 Basic Linear Equations 1 

2 PINION AND GEAR GENERATION 12 

2.1 Pinion Generation 12 

2.2 Gear Generation 13 

2.3 Gear Machine-Tool Settings 13 

3 GEAR GEOMETRY 16 

3.1 Gear Surface 16 

3.2 Mean Contact Point and Gear Principal Directions and Curvatures 20 

4 LOCAL SYNTHESIS OF SPIRAL BEVEL GEARS 23 

4.1 Conditions of Synthesis 23 

4.2 Procedure of Synthesis 23 

5 PINION MACHINE-TOOL SETTINGS 28 

5.1 Introduction 28 

5.2 Head-Cutter Surface 28 

5.3 Observation of A Common Normal At The Mean Contact Point For Surfaces E p , S 2 , 

and Ei 31 

5.4 Basic Equations For Determination of Pinion Machine-Tool Settings 33 

5.5 Determination of Cutter Point Radius 37 

5.6 Determination of mpi = and Xgi 39 

xtap 


..JNftNMfllUl SH 


111 




BLANK ftor FILMED 



6 TOOTH CONTACT ANALYSIS 


48 


6.1 Introduction 48 

6.2 Gear Tooth Surface 50 

6.3 Pinion Tooth Surface 51 

6.4 Determination of Transmission Errors 61 

6.5 Simulation of Contact 62 

7 V and H Check 69 

7.1 Determination of V and H values 69 

7.2 Tooth Contact Analysis for Gears With Shifted Center of Bearing Contact 73 

A Generation With Modified Roll 74 

A.l Introduction 74 

A. 2 Taylor Series for The Function of Generating Motion 74 

A. 3 Synthesis of Gleason’s Cam 78 

A. 4 Cam Analysis 82 

A. 5 Determination of Coefficients of The Taylor Series 92 

A. 6 Selection of Cams and Cam Settings 95 

B Description of Program and Numerical Examples 97 

References 100 


iv 

» 4 ‘'- 



NOMENCLATURE 


[A] 

A* 

A m 

®ij ( * > / = 1 » 2 ) 3 ) 

a, b 

bG 

c 

C 

D 

E 

F 

6CX 
24 DX 
120 EX 
Emi 

e/>4 
e« , e^ 

e «p2, e qp2 

Hg,V g 

h-m 

i 

j 

4 ° (* = 1 , 2 ) 


Augmented matrix of linear equation system 

Distance between the shifted center of bearing contact and 

pitch apex 

Mean cone distance 

Coefficients of basic linear equations 

Half-long and short axes of contact ellipse 

Gear dedendum 

Clearance 

Coefficient of the second order of Taylor series of generation motion 

Coefficient of the third order of Taylor series of generation motion 

Coefficient of the forth order of Taylor series of generation motion 

Coefficient of the fifth order of Taylor series of generation motion 

Third order parameter of generation motion 

Forth order parameter of generation motion 

Fifth order parameter of generation motion 

Blank offset in generation of gear i 

Principal directions of surface Ej 

Principal directions of surface E 2 

Principal directions of gear surface in system S p 2 

Gear horizontal and vertical settings 

Mean whole tooth depth 

Tilt angle 

Swivel angle 

Sum of principal curvatures of surface 1 or 2 


v 



K?\ 4 / (* = 1 . 2 ) 

[-t'yx] 


W*2l(^l) 

M 

[My*] 


Ni (* = 1,2) 

n p 2 

n 

Oi 

O2R 

P 

PW 

9* 

Rep 

RaG 

Ru2 

r c 

r P 2 

n 

r h 

J.°f) 

T h 

ri(9 F ,<t>F) 

F2(Sg,(P p ) 


Principal curvatures of surface 1 or 2 

Matrix of coordinate transformation from system to system 5 y 
for free vectors 

Derivative of transmission ratio 
Mean contact point 

Matrix of coordinate transformation from system S x to system S y 
for position vectors 

Number of teeth of pinion (i = 1) or gear (i = 2) 

Unit normal vector of gear cutter surface in system S p 2 
Common unit normal at point of contact 
Pitch cone apex of gear i 
Root done apex of gear 

Percentage of amount of shift along the pitch line over face width 

Point width of gear cutter 

Cradle angle for gear i 

Point radius of pinion head cutter 

Gear ratio of roll 

Gear nominal cutter radius 

Gear cutter tip radius 

Position vector of gear cutter surface in system S p 2 

Position vector of tooth surface of gear i represented in system Si, 

is equivalent to [r,-] 

Position vector of mean contact point in system Sh 
Position vector of pinion cutter center in system Sh 
Position vector of pinion in system Si 
Position vector of gear in system S 2 


vi 



S a > Sf, 
Smi 

Sci 

Sh 

Si 

(•sg>#g) 

(sf,8f) 

Sri 

s x 

AT 
f?M 


Vr W (t = 1,2) 

V t ? (i = 1,2) 

(* = 1 , 2 ) 


Vi 


12 


^P 2 ) 

v m2 

{Xo \ zP) 


X Gi 

(XL,RL) 

V G ,H G 

Zr 

S J 

(p) (p) 

A/g , /Vq 

Kj , Kh 

Imi 


Coordinate systems originated at point of contact between Si and E 2 
Coordinate system rigidly connected to the cutting mac hin e of gear i 
Movable coordinate system rigidly connected to the cradle of cutting 
machine for gear i 
Fixed coordinate system 

Coordinate system rigidly connected to gear i 
Surface coordinates of gear cutter surface 
Surface coordinates of pinion cutter surface 
Radial setting of gear i 

Auxiliary coordinate system identified by subscript x 
Cam setting 

Sliding velocity of contact point in the motion over surface E i 
Transfer velocity of contact point in the motion with surface E i 
Projection of 0*^ upon e„ and e q 
Relative velocity at contact point 

Relative velocity in the process for gear generation represented in system 

Coordinates of center of the arc blade 

Sliding base for generation of gear i 

Machine center to back for generation of gear i 

Parameters determining mean contact point 

Vertical and horizontal adjustments for the gear drive 

Distance of gear root cone apex beyond pitch cone apex 

Principal curvatures of surface E 2 

Principal curvatures of surface E 2 

Principal curvatures of surface Ei 

Machine root angle for generation of gear i 

vi i 



Ti 

7 « 
a 

a G ,ap 

P 

(\,6 F ) 


Vi 

(v,() 

8 

<t>i 

ti 

<t>F, 4>p 

V’G? 

<r( 12 ) 

Si 

Yi F 

£ P 

(i = l,2) 

u( F \u}( p ) 


d pz ) 


u>, 


m2 


Pitch angle of gear i 
Root angle of gear i 
Cam guide angle 

Cutter blade angles for gear and pinion respectively 
Radius of circular arc 

Surface coordinates of the surface of revolution generated 
by circular arc blade 

Direction angle of contact path on surface Sj 
Unit vectors along long and short axes of contact ellipse 
Elastic approach 
Gear dedendum angle 

Angle of rotation of gear i in the process for generation 

Rotation angle in meshing of gear i between the gear (2) and the pinion (1) 

Rotation angles of cradle in the process for pinion and gear 

generation, respectively 

Gear spiral angle 

Surface coordinates of gear tooth surface at mean contact point 
Angle formed between principal directions ef and e„ (in meshing 
and generation ) 

Surface of gear i 
Pinion generating surface 
Gear generating surface 

Angular velocity of surface E t - (in meshing and generation ) 

Angular velocity of the cradle in the process for pinion and 
gear generation, respectively 

Relative angular velocity in the process for gear generation represented 


vi i i 



in system S m 2 

u> ( iI ’ 1 ) Relative angular velocity in the process of pinion generation 

£(*) Angular velocity of gear i 

uJ(*i) Relative angular velocity between gear i and gear j 


IX 



SUMMARY 


Computerized simulation of meshing and bearing contact for spiral bevel gears and hypoid gears 
[1,2] is a significant achievement that could improve substantially the technology and the quality of 
the gears. This report covers a new approach to the synthesis of face-milled spiral bevel gears and 
their tooth contact analysis. The proposed approach is based on the following ideas proposed in [3] 
(i) application of the principle of local synthesis that provides optimal conditions of meshing and 
contact at the mean contact point M and in the neighborhood of M; (ii) application of relations 
between principle directions and curvatures for surfaces being in line contact or in point contact. 

The developed local synthesis of gears provides (i) the required gear ratio at M; (ii) a localized 
bearing contact with the desired direction of the tangent to the contact path on gear tooth surface 
and the desired length of the major axis of contact ellipse at M; (iii) a predesigned parabolic function 
of a controlled level (8-10 arc seconds) for transmission errors; such a function of transmission errors 
enables to absorb linear functions of transmission errors caused by misalignment [3] and reduce the 
level of vibrations. 

The proposed approach does not require either the tilt of the head-cutter for the process of 
generation or modified roll for the pinion generation. Improved conditions of meshing and contact 
of the gears can be achieved without the above mentioned parameters. The report is complemented 
with a computer program for determination of basic machine-tool settings and tooth contact anal- 
ysis for the designed gears. The approach is illustrated with a numerical example. 

The contents of the following sections cover the following topics: 

(1) . Basic ideas of local synthesis of gears and the mathematical concept of this approach 
(Chapter 1). The local synthesis discussed in this chapter is applicable for all types of gears and 
provides the optimal conditions of meshing and contact at the mean point of tangency of gear tooth 
surfaces. 

(2) . Methods for generation of the pinion and the gear and basic machine-tool settings that are 


x 



necessary for gear generation (Chapter 2). 

(3) . Determination of geometry of geaT tooth surface, the gear mean contact point and the 
principal directions and curvatures at this point (Chapter 3). 

(4) . Application of basic principles of local synthesis for spiral bevel gears (Chapter 4). 

(5) . Determination of pinion machine-tool settings considering as given:(i) the gear geometry, 
and (ii) the conditions of meshing and contact at the mean contact point obtained from the local 
synthesis (Chapter 5). 

(6) . Computerized simulation of meshing and contact (Tooth Contact Analysis) for spiral bevel 
gears that have been synthesized in the previous chapters (Chapter 6). 

(7) . Analysis of the shift of bearing contact caused by the misalignment of gears (Chapter 7). 

(8) . The theory of modified roll (variation of cutting ratio in the process for generation) and 
mechanisms used for application of modified roll (Appendix A). 

(9) . Description of developed computer programs and numerical examples that illustrates the 
application of those programs. 
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1 Local Synthesis of Gears (General Concept) 


1.1 Introduction 

The main goals of local synthesis are to provide: (i) contact of gear tooth surfaces at the mean point 
of contact of gear tooth surfaces, and (ii) improved conditions of meshing within the neighborhood 
of the mean contact point. The local synthesis is the first stage of the global synthesis with a goal 
to provide improved conditions of meshing for the entire area of meshing. The criteria of conditions 
of meshing are the transmission errors and the bearing contact. The principles of local synthesis 
that are discussed in this chapter for face-milled spiral bevel gears can be applied for other types 
of gears as well. 

1.2 Basic Linear Equations 

Consider two right-handed trihedrons 5 a (ey,e/ l ,n) and Sb(e a ,e q ,n ) (Fig. 1.2.1). The common 
origin of the trihedrons coincides with the contact point M, the n-axis represents the direction of 
the surface unit normal, ej and are the unit vectors of the principal directions of surface Ex, e s 
and e q represent the principal directions of surface £ 2 ? and cr( 12 ) is the angle formed between ej 
and e s (measured clockwise from e a to ef and counterclockwise from ej to e, ). In reference [4] 
three linear equations were derived that relate the velocity v^ 1 of the contact point over surface Ex 
with the principal curvatures and directions of contacting surfaces and the transfer components of 
velocities. These equations are: 


(1) , (1) 

a n vl + a i2 Vq = oi 3 
(1) . (1) 

OX2^« + 022 Vq — 0 2 3 " 

(1) 1 l 1 ) 

0x3^8 + 0 23 X^ =033 


( 1 . 2 . 1 ) 


Here (see the designations in [4] ) 
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<*13 

= 

<*31 = -K.X 12) - 

6j( 12 ^ne 8 


<* 22 

= 

faq 

-K f sin 2 cr( 12 )- 

- Kh COS 2 1 

cr( 12 ) 

a 23 

= 

’ <*32 = -KqV { q 12) - 

^ l2 ' > ne q J 


<*33 

= 

Kf g 


2 

n w (12) 


+ 

{ 

m 2l(** X *2) 

' (* 

'-*)} 

v?> 

= 

• e g 





= 






( 1 . 2 . 2 ) 


- n (u! (1 ) x x 


Equations (1.2.1) and (1.2.2) can be applied for two cases where: (i) surfaces Si and S 2 are 
in line contact, and (ii) the surfaces are in point contact. The instantaneous line of contact is 
typical for the case when the gear tooth surface (£ 1 ) is generated by the tool surface (£ 2 )- The 
instantaneous point of contact is typical for gears with localized bearing contact. 

Line Contact 

When the gear tooth surfaces are in line contact, the direction of velocity tJr 1 ^ can be varied, and 
equations (1.2.1) can not provide a unique solution for the unknowns uj; 1 ^ and Ug 1 ^. This results 
in that the rank of the augmented matrix 


M = 


a 11 

<*12 

<*13 

<*12 

<*22 

<*23 

<*13 

<*23 

<*33 


(1.2.3) 


must be less than 2. This requirement yields 



2 \ 
a 12 — a ll a 22 

ana23 = anai3 

<112033 = 013023 . 

Equivalent equations are 



033 


012 


013023 

033 


022 = 


x 23 

033 


(1.2.4) 


(1.2.5) 


Using equations (1.2.5) and (1.2.2) we obtain equations that will enable us to determine 
Kf and Kh for Ei considering as given k„ and K q for surface £ 2 - The equations are: 


tan2cr( 12 ) = 


2Ql3<Z23 

O23 — 0 2 g + ( k s ~ ^9)033 


(1.2.6) 


_ __ 2a 13 a 2 3 

K * Kh a 33 sin2tr( 12 ) 


(1.2.7) 


Kf ~j“ Kfi — ( Kg "b Kg) 


2 2 
a 13 ~t~ O23 

033 


( 1 . 2 . 8 ) 


Equation (1.2.6) provides two solutions : (r[ 12 ' 1 and = cr[ 12 ^ + 7r/2 and both of them can be 
used for computations of Kf and Kh that are represented by equations (1.2.7) and (1.2.8). Fig. 1.2.2 
shows the orientation of two couples of unit vectors ( i = 1,2), with respect to unit vector 

e». The magnitude of principal curvature for the direction with collinear vectors e j 1 and is the 
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same (k^ = k^) although the notation for the unit vectors has been changed. Similarly, we can 
say that 

Knowing the angle a^ 12 \ and the unit vectors e a and e q , the principal directions on surface Ei 
can be determined with the following equations, 


.-O) 


•j=«/ = COS CT 


( 12 )e - 


sm cr 


( 12 ): 


(1.2.9) 


5<1> 


II = e h 


sincr^ 12 ^e. 


+ cos & 


( 12 ); 


(1.2.10) 


Point Contact 

In the case of instantaneous point of contact, the direction of motion of the contact point over 
the surface is definite, equations (1-2.2) for the unknowns can provide a unique solution for the 
unknowns and Uq 1 * and the rank of matrix [A] is 2. This yields that 


an 

ai2 

<*13 



a 21 

a22 

<*23 

= 0 

(1.2.11) 

<*31 

a 32 

<133 




Equation (1.2.11) yields the following relation 

f(ns,K q ,Kf,K, h ,a( u \m' 12 ) = 0 ( 1 . 2 . 12 ) 


Our goal is to determine and er^ 12 ) (the principal curvatures 

provide at the mean contact point (i) a certain direction of the tangent 


and directions of Si) and 
to contact path on surface 
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S 2 j (ii) a desired length of the major axis of instantaneous contact ellipse, and (iii) a parabolic 
function of transmission errors. For these purpose we have to derive extra equations in addition to 
equation (1.2.12) 

Determination of m' 21 

The derivative is the second derivative of function <^> 2 (</> 1 ) that is taken at the mean 

contact point; <f>i and (f> 2 are the angles of rotation of gears 1 and 2. In the case of an ideal gear 
train, function <f> 2 [<f>i) is linear and is represented by 

4>2 = 4>i^r ( L2 - 13 ) 

iv 2 

However, due to misalignment between the meshing gears the real function <f> 2 ((/>i) becomes a 
piecewise periodic function with the period equal to the cycle of meshing of a pair of teeth (Fig. 
1.2.3). Due to the jump of angular velocity at the junction of cycles, the acceleration approaches 
to an infinitely large value and this can cause large vibration and noise. For this reason it is 
necessary to predesign a parabolic function of transmission error that can absorb a linear function 
of transmission error and reduce the jump of angular velocity and acceleration [3]. This goal(the 
predesign of a parabolic function) can be achieved with certain relations between the principal 
curvatures of contacting surfaces . 

Fig. 1.2.4 shows the predesigned transmission function for the gear convex side (Fig. 1.2.4(a)) 
and gear concave side (Fig. 1.2.3(b)). Both functions-^ (<?!>;[) and 0 2 ^(<£i)- are in tangency at the 
mean contact point and have the same derivative m 2 j , at this point. 

Consider now that the predesigned transmission function is represented as 

fo - ^ 0) = F{4>\ ~ (1-2-14) 


5 



Here: (f>^ and are the initial angles of rotation of gears 1 and 2 that provide the tangency 
of gear tooth surfaces at the mean contact point M . 

Using the Taylor expansion up to the members of second order, we obtain 




2 d<j>i 


= ™2l(<£l - 4 >( l ) ) + ^2l(01 - <^l 0) ) 2 


(1.2.15) 


where m 2 i(4>x) is equal to N 1 /N 2 at the mean contact point and m' 21 is the to be chosen constant 
value: positive for the gear concave side, and negative for the gear convex side. The synthesized 
gears rotates with a parabolic function of transmission errors represented by 


= 2 m 2i(<^i - 4 0 ') 2 


(1.2.16) 


where 


< Ux - < — 

Nl ~ ^ ~ N x 

Equation (1.2.16) enables the determination of m 21 considering as known the expected values 
of transmission errors. 

Relation between Directions of Paths of Contact 

We recall that velocities and vi" 2 ^ are related by the equation [4], 



(1.2.17) 


Directions of velocities v,- 1 ^ and v T 2 ' 1 coincide with the tangents to the contact path that form 


angles 7 ? 1 and r/ 2 with the unit vector e s (Fig. 1.2.5). Equations (1.2.17) yield 
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,( 2) = v w + v w 


(1.2.18) 




+ v 


( 12 ) 

q 


According to Fig. 1.2.5 


Vg l ) = tan t]i 


Third equation of system (1.2.1) and equations (1.2.18) and (1.2.19) yield 


tan 7 ?! = 


-G3i4 12) + ( a 33 + a 3 ii>i 12 ^) tan 772 


0-33 + Cl32 { v q 


~ V ' 3 


tan 772 ) 


(1.2.19) 


( 1 . 2 . 20 ) 



033 

a 13 + a 23 tan hi 


( 1 . 2 . 21 ) 



Q33 tan 771 
«i 3 + a 23 tan 771 


( 1 . 2 . 22 ) 


Prescribing a certain value for 772 (choosing the direction for path of contact on E 2 ), we can 
determine tan 771, vi ^ and . We recall that coefficients 031,032 and 033 do not depend on the 
to-be determined principal curvatures Kf and and o( 12 ). 

Relations between the Magnitude of Major Axis of Contact Ellipse, Its Orientation and 
Principal Curvatures and Directions of Contacting Surfaces 

Our goal is to relate parameters and Kh. of the pinion surface Ei with the length of the 

major axis of the instantaneous contact ellipse. This ellipse is considered at the mean contact point 
and the elastic approach S of contacting surfaces is considered as known from the experimental 
data. The derivation of the above mentioned relations is based on the following procedure 
Step 1: Using equations (1.2.2), we obtain 
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a n + a 22 = = K s 

ail - a 22 = 92-91 cos2<r( 12 ) 

(an - a 22 ) (i) 2 + 4a 2 2 = g\ - 2g x g 2 cos 2 cr( 12 ) + g\ J 


(1.2.23) 


Step 2: It is known from [4] that 


a = 


(1.2.24) 




- \Jg\- 2g 1 g 2 cos2<r + g\ 


(1.2.25) 


Equation (1.2.25) yields 


[( a n + a 22 + 4A] 2 — (an — a 22 ) 2 + 4a 2 2 


(1.2.26) 


Step 3: We may consider now a system of three linear equations in unknowns an,ai 2 and a 22 

(i) , (i) 

v\ an + v\ 'a X2 = ai 3 

(i) , (i) 

vl an + Vq ’ a 22 = a 23 

an + a 22 = K-£ j 

Step 4: The solution of equation system (1.2.27) for the unknowns an,ai 2 and a 22 allows to 
express these unknowns in terms of ai 3 ,a 23 , Ks,vi^ and Then, using equation (1.2.25) we 
can get the following equation for Ky 


(1.2.27) 
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Here: 


K ^ - ("? + "D 

2 A — (ni cos 2tji + n 2 sin 2??i ) 


(1.2.28) 




a 13 ~ a 23 tan2 >7l 
(1 + tan 2 7ji)a 33 


a 2 


a i 3 tan 771 + a 23)( a i 3 + a 2 3 tan 771 ) 
(1 + tan 2 7?i)a 33 


A = 


6 


(1.2.29) 


The advantage of equation (1.2.28) is that we are able to determine K s knowing the major axis 
2a of the contact ellipse and the elastic approach S. 

Step 5: The sought for principal curvatures and directions for the pinion identified with k/,k/, 
and cr( 12 ) can be determined from the following equations 



- 


(1.2.30) 


tan2<r( 12 ) = 


2022 

92 ~ ( a ll — a 22) 


2ri2 — A”s sin 27/1 
02 - 2ni + AT S cos 2^ 


(1.2.31) 


2a!2 2n2 — ATs sin2r?i 

sin2<r( 12 ) “ sin2o-( 12 ) 


(1.2.32) 


(1) , 

_ (1) + <7 1 

*/ = */ = -*-2 


_ (1) 4 1} - 91 

K h = ^ 


(1.2.33) 


(1.2.34) 
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Step 6: The orientation of unit vector ej and is represented with equations (1.2.9) and 
(1.2.10). The orientation of the contact ellipse with respect to e) is determined with angle 
(Fig. 1.2.6) that is represented with the equations 


cos 2a^ = 
sin2a^ = 


9 \ - 92 cos2(t( 12 ) 

[g\ - 2gi02 cos 2 cr( 12 > + g\) 2 

gi sin2cr( 12 ) 

(ffi - 29192 cos 2 cr( 12 ) + gj)i 


(1.2.35) 

(1.2.36) 


The minor axis of the contact (26) ellipse is determined with the equations 


6 = 



(1.2.37) 



i4 2) + \J 9 1 - 2^152 cos2er + 


92 


(1.2.38) 


Local Synthesis Computational Procedure : 

The following is an overview of the computational procedure that is to-be used for the local 
synthesis. 

The input data are: k s , n q , e s , e q , f( M ), £( 12 \ #( 12 ) and 6. The to-be chosen parameters are: 
Tj 2 ,m' 21 and 2 a. The output data are: «/,«/,, o^ 12 ), e/ and eh- 
Step 1: Choose r? 2 and determine tjj from equation (1.2.20) 

Step 2: Determine and from equations (1.2.21) and (1.2.22) 
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Step 3: Determine A from equation (1.2.29) 


Step 4: Determine K s from equation (1.2.28) 

Step 5: Determine and k, h by using the set of equations from (1.2.30) to (1.2.34) 

Step 6: Determine the orientation of the contact ellipse and its minor axis by using equations 
from (1.2.35) to (1.2.37) 

1.3 Conclusion 

The contact of tooth surfaces is considered for two cases: line contact and point contact. For line 
contact, the principal directions and curvatures of one surface can be determined in terms of the 
other’s knowing the relative motion between the two . For point contact, we proposed an approach 
for local synthesis of spiral bevel gears which enables: (i) to provide a limited level of transmission 
errors, (ii) optimal direction for the path of contact on gear surface S 2 ,and (iii) the guaranteed 
length of the major axis of contact ellipse. 

The output data obtained from the procedure of local synthesis are: Kf , Kh, ej and e),. 
The machine-tool settings for the generation of the gear tooth surfaces must be carefully chosen to 
guarantee the above mentioned conditions of local meshing and contact. 
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2 Pinion and Gear Generation 


2.1 Pinion Generation 

To describe the pinion generation we will use the following coordinate system (Fig. 2. 1.1): (i) S rn i~ 
a fixed coordinate system that is rigidly connected to the cutting machine; (ii) -a movable 
coordinate system that is rigidly connected to the cradle and performs rotation with the cradle 
about the Z m j- axis; initially, S c i coincides with S m i (Fig. 2. 1.1 (b)); angle <f>p determines the 
current position of S c i (Fig 2.1.1 (c)) : (iii) Coordinate systems S a and Sb that are rigidly connected 
to the cradle and its coordinate system S c i; systems S a and Sb are used to describe the installment 
of the head-cutter on the cradle. Angle q\ determines the orientation of S a with respect to S c i', 
(iv) Coordinate system Sp that is rigidly connected to the head-cutter (not shown in Fig.2.1.1); 
the head-cutter in the process for generation performs rotation with the cradle (transfer motion) 
and relative motion with respect to the cradle about an axis that passes through O a \ (v) Auxiliary 
coordinate systems Sd and S p are used to describe the installment of the pinion on the cutting 
machine (Fig.2.1.1 and Fig.2.1.2); the pinion axis forms angle 7 ml with axis X d that is parallel to 
I ra i. (vi) A movable coordinate system Si that is rigidly connected to the being generated pinion; 
the pinion rotates about the axis X p and <j>\ is the current angle of pinion rotation (Fig.2.1.2). 

Henceforth, we have to differentiate the parameter of motions that are performed in the process 
for generation and the parameters of installment of the head cutter and the pinion on the cutting 
machine. 

In the process for generation the cradle of the cutting machine with the mounted head-cutter 
performs rotation with angular velocity a>( F ) (Fig.2.1.2). The head-cutter performs rotational 
motion with respect to the cradle but this motion is not related with the process for generation 
and just provides the desired velocity of cutting. The being generated pinion performs rotational 
motion with angular velocity a^ 1 ) (Fig.2.1.2) that is related with . 

The parameters of installment of the head-cutter are: (i) the swivel angle j (Fig.2.1.1) and the 
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tilt angle i that is the turn angle of S t about Yj, (Fig. 2. 1.3); S T i = |O c O m i| is the radial setting; q\ 
is the cradle angle. 

The parameters of installment of the pinion are: F m i-the shortest machine center distance (Fig. 
2.1.1, Fig. 2. 1.2); root angle 7 ml ; sliding base Xr\\ machine center to back Xgi- 

2.2 Gear Generation 

While describing the gear generation, we will consider the following coordinate systems: (i) S m 2 
that is rigidly connected to the cutting machine; (ii) S C 2 that is rigidly connect to the cradle, (iii) S p 2 
that is rigidly connected to the head-cutter and 5 c2 ; (iv) Sd 2 that is an additional fixed coordinate 
system rigidly connected to S rn 2 ; and (v) S 2 that is rigidly connected to the being generated gear. 

The cradle performs rotation about the Z m 2 axis with angular velocity uJb') (Fig.2.2.1). The 
initial and current positions of coordinate systems S c 2 and S p 2 with respect to S m 2 are shown in 
Fig.2.2.1 (a) and Fig.2.2.1 (b), respectively. 

Coordinate system Sd 2 (it is rigidly connected to S m 2 ) is used to describe the installment of 
the gear at the cutting machine (Fig. 2. 2. 2(a)). In the general case apices O 2 R and O 2 of the gear 
root cone and pitch cone do not coincide. Apex O 2 R is located on axis X m 2 of the cutting machine. 
The origin 0d2 of $d 2 coincides with the apex O 2 of the gear pitch core. Axes Xd 2 and X m 2 form 
angle 7 m2 which is the gear machine root angle. 

Coordinate system S 2 is rigidly connected to the gear that in the process of generation performs 
rotation about Xd 2 with angular velocity (Fig. 2. 2. 2(b)). Angle 4>2 is the current angle of 
rotation of gear 2. 

2.3 Gear Machine Tool Settings 
Gear Cutting Ratio 

Fig.2.3.1 shows the sketch of the gear with noncoinciding apexes of the root and pitch cones. 
In the process for generation the pitch line 0 2 P is the instantaneous axis of rotation. It is evident 
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that the angular velocity of rotation in relative motion, u 5 ( p2 ) , must lie in the plane that is formed 
by vectors Q ^ and u>( 2 ) (Fig.2.3.2) 


j(P 2 ) = ajW - £( 2 ) 

(2.3.1) 

The cutting gear ratio is: 


D |u;^ 2 )| COS^g COs(T 2 - 72 ) 

RaG - |u5(p)| "" sinr 2 “ sin T 2 

(2.3.2) 

Gear Settings 

Fig.2.3.3 shows the installment of the head- cutter. We designate the 
0 2 P (Fig. 2. 3.1, Fig.2.3.3) by A m . Then we obtain (Fig.2.3.3) 

mean pitch cone distance 

Hg - A m cos 8 g - R U 2 sin V>G 

(2.3.3) 

Vg ~ Ru 2 cos rj) G 

(2.3.4) 

S T 2 = {Hh + V£)* (Sr2 = \O m2 O p2 \) 

(2.3.5) 

. Vg 

q 2 = sm — 
Or 2 

(2.3.6) 

Here: tpG is the spiral angle on the root cone, R u2 is the mean radius of the head cutter. The 
sliding base | 0 m 2 0 2 | is 

Xb 2 = Zr sin 7^2 

(2.3.7) 
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Here: 7 m 2 is the same as the gear root cone angle 72. and Zr is the distance between O2R and 
0 2, which are the apexes of the root cone and the pitch cone, respectively. 
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3 Gear Geometry 


3.1 Gear Surface 

The gear tooth surface is the envelope to the family of generating surfaces. We recall that the 
cradle carries the head-cutter that is provided with finishing blades. The blades are rotated about 
the axis of the head-cutter and generate two cone surfaces. Fig.3.1.1 shows one of the cones. 

The fa mil y of a generating surface (the cone surface) is generated in S 2 while the cradle and 
being generated gear perform related rotations, about the Z m 2 -axis and X 2 - axis (Fig.2.2.2). 

The derivation of the gear tooth surface is based on the following procedure: 

Step 1: We represent the cone surface and its unit normal in system S p2 (Fig.3.1.1) as follows 


i.e. 


f P 2 = 


(r c — sin olg) cos Oq 
(r c — sg sin etc) sin #g 
— sg cos ac 
1 


n P 2 = 


N p 2 

\Npt\ 



dr p 2 

ddc 


x 


fp2_ 

ds G 


n P 2 


cos c*g cos 6q 
cos olg sin 9q 


smac 


(3.1.1) 


(3.1.2) 


(3.1.3) 


Here: sq and 6q are the surface coordinates; ac is the blade angle; r c is the radius of the 
head-cutter that is measured at the bottom of the blades. It is evident (Fig. 3. 1.2) that 


r c = R u 2 ± 


PW 

2 


(3.1.4) 
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Here: R u 2 is the nominal radius, PW is the so called point width; the positive sign in (3.1.4) 
corresponds to the gear concave side and the negative sign corresponds to the gear convex side. 

Equations (3.1.1) and (3.1.3) represent both generating cones with ag > 0 for the gear convex 
side and ag < 0 for the gear concave side. 

Step 2: The family of generating surfaces that is generated in S 2 is represented by the following 
matrix equation 

r 2 @G > tfip) = \M2d 2 \lMd2mi\\MTn2Cz\\Mc2p2\ r P'2 (3.1.5) 

Here (Fig.2.2.2, Fig.2.2.1): 

10 0 0 

0 cos <f >2 sin 0 

0 — sin <f> 2 cos (j> 2 0 

0 0 0 1 

cos 7 m2 0 sin q m2 --XB 2 sin 7 m2 

0 10 0 

_ sin 7 m 2 0 cos 7 m2 -XB 2 COS 7 m2 

0 0 0 1 





(3.1.6) 
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1 0 0 S r2 cos g 2 


[^C2Pi] 


0 1 0 S r 2 sin g 2 

0 0 1 0 

0 0 0 1 


[-^771202] — 


cos <f) p — sin <f> p 0 0 

sin <f>p cos <f> p 0 0 
0 0 10 

0 0 0 1 


(3.1.8) 


(3.1.9) 


The machine root angle 7 m2 in equation (3.17) is equal to gear root cone angle 7. 

Step 3: The derivation of the equation of meshing is based on the equation 

*W-^ 2 2) = 0 (3.1.10) 


The subscript “m2” means that vectors in equation (3.1.10) are represented in coordinate system 
Sm 2 ; «m2 is the unit normal to the generating surface; = v^\ - v^\ is the relative (sliding) 
velocity. Vector n m2 is represented by the matrix equation 


U m2 — [f'rri2P2 3^'p2 


- COS OCQ COS ($G + <t> p) 

- cos c*g sin(^G + 4> p ) 

sin ac 


(3.1.11) 


where [ Tm 2 p 2 ] is the 3x3 submatrix [Mmjpj]. 

We consider that the axes of rotation of the cradle and the gear intersect each other (Fig. 2. 2. 2(a)), 
thus 
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(3.1.12) 


where 


v' p2 ' > = xf — (3(») 

‘m2 r™ 2 m 2 J A ' m 2 — “Vj 


m 2 ^ p m 2 


0>, 


m 2 2) = [- COS 72 o - sin 72 )] 3 


(3.1.13) 


We assumed that \uj^ l2 | = 1 in equation 
following relation 


(3.1.13). Equations from (3.1.10) to (3.1.13) yield the 


Here 


A{6 G , <t> P ) 

B{0gAp) 


(3.1.14) 


A(6 g , 4> p ) = n m2a .[-i4, 1 (sin72 - -^— )] + 7 i TO2W [Xb 2 COS72 + A 2 (sin 7 2 - -^— )] 

■K-aG -ttaG 

+nm 2 zAiC0S'f2 (3.1.15) 


B{dG,<t>v) = ~ n m 2 x sinac sin(# G 4- <j> p ) + n m2y [sin7 2 — — — ) sin a G cos(6 G + 4> P ) 

tt-aG 

- cos a G cos 72] 4- Timiz cos ~f 2 sin a G sin (0 G + <t> p ) (3.1.16) 


A x = r c sm(6 G + <t> p ) - S r2 sin(g 2 - <f> P ) (3.1.17) 


A 2 = r c cos (9 g + 4> P ) - S r 2 cos( 5 2 - <j> P ) (3.1.18) 


Step 4: Equations (3.1.5) and (3.1.14) considered simultaneously represent the gear surface in 
three- parametric form but with related parameters. Since parameter s G in equation of meshing 
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(3.1.14) is linear, it can be eliminated in equation (3.1.5), and then the gear tooth surface will be 
represented in two-parametric form, by the vector function r 2 (9 G ,<f> p ) ■ 

3.2 Mean Contact Point and Gear Principal Directions and Curvatures 

The mean contact point M is shown in Fig.2.3.1. Usually, M is chosen in the middle of the tooth 
surface. The gear tooth surface and the pinion tooth surface must contact each other at M . 

The procedure of local synthesis discussed in section 2.1 is directed at providing improved 
conditions of meshing and contact at M and in the neighborhood of M . The location of point 
M is determined with parameters XL and RL (Fig.2.3.1) that are represented by the following 
equations 

XL = A m co S r 2 + (b G — sinr 2 (3.2.1) 


RL — A rn sin T 2 - (b G - C )cosF 2 

Z 


(3.2.2) 


Here: A m is the pitch cone mean distance; h m is the mean whole depth; b G is the gear mean 
dedendum; c is the clearance Equations (3.2.1), (3.2.2) and vector equation r 2 (6 G ,4> p ) for the gear 
tooth surface allows to determine the surface parameters 6q and <j) p for the mean contact point 
from the equations 


* 2 (0c,$) = XL 


(3.2.3) 


Yi(0b<t>;) + zi(eh,<t>* p ) = (RL) 2 


(3.2.4) 
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Gear Principal Directions and Curvatures 


The gear principal directions and curvatures can be expressed in terms of principal curvatures 
and directions of the generating surface (see chapter (13) in [4]), that is the cone surface. 

Step 1 : The cone principal directions are represented in S p 2 by the equations (see (3.1.1)) 

dr^ 

(3.2.5) 


4p 2 = -II s - = [— sin#c COS0G 0] r 

I UT P2 I 


36g 


dr. 


P2 


r\ 

4 p 2 = = [- s i R a G cos 6q - sin ac sin 6q - cos a G ] r 


P2 


(3.2.6) 


8sg 


The superscript “p” indicates that the cone surface S p is considered. Unit vector e£) 2 is directed 
along the cone generatrix and unit vector ef P2 is perpendicular to e$p 2 . The unit vectors of cone 
principal directions are represented in 2 by the equations 


e<p) 

8TT12 


[- sin(0G + <f> P ) cos (9 g + 1 i> p ) 0] 3 


(3.2.7) 


4m 2 = [- sinaGCOs(^G + <P P ) - sin a g sin (0 G + <t> P ) - cos a G ] T (3.2.8) 


The cone principal curvatures are: 

(3.2.9) 


Step 2 : The determination of principal curvatures and directions for gear tooth surface S 2 is 
based on equations from (1.2.6) to (1.2.8). The superscript “2” in these equations must be changed 




cosac 
- sg sin qg 


and 


4 p) = 0 
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for “p” and superscript “1” for “2”. The second derivative of cutting ratio, m' 21 = is zero 

because the cutting ratio is constant. The principal curvatures of the gear tooth surface will be 
determined as Kf and The principal directions on gear tooth surface will be represented in by 
ej and and they can be determined from equations (1.2.9) and (1.2.10). To represent in 52 the 
principal directions on gear tooth surface S 2 and its unit normal we use the matrix equation that 
describe the coordinate transformation from 5 m2 to 5 2 . This equation is 

a 2 = [T 2 ti 2 ][T(i 2 m2] a m2 (3.2.10) 

Here: a m2 stands for vectors n m2 , e f TIL2 and e) im2 , and a 2 stands for n 2 , efy an d 
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4 Local Synthesis of Spiral Bevel Gears 


4.1 Conditions of Synthesis 

The basic principles of local synthesis of gear tooth surfaces discussed in Section 1 will enable us 
to deter min e the principle curvatures and directions of the being synthesized pinion. Thus, we will 
be able to determine the required machine-tool settings for the pinion. While solving the problem 
of local synthesis, we will consider as known: 

(i) The location of the mean contact point M in a fixed coordinate system, and the orientation 
of the normal to gear surface £ 2 - 

(ii) The principle curvatures and directions on £2 at M . The local synthesis of gear tooth 
surfaces must satisfy the following requirements: 

(1) The pinion and gear tooth surfaces must be in contact at M. 

(2) The tangent to the contact path on the gear tooth surface must be of the prescribed direction. 

(3) Function of gear ratio m 2 i((/>i) in the neighborhood of mean contact point must be a linear 
one, be of prescribed value at M and have the prescribed value for the derivative m' 21 (<f>i ) at M. 
The satisfaction of these requirements provides a parabolic type of function for transmission errors 
of the desired value at each cycle of meshing. 

(4) The major axis of the instantaneous contact ellipse must be of the desired value (with the 
given elastic approach of tooth surfaces). 

4.2 Procedure of Synthesis 

We will consider in this section the following steps of the computational procedure: (i) representa- 
tion of gear mean contact point in a fixed coordinate system Sh ; (ii) satisfication of equation of 
meshing of the pinion and gear at the mean contact point ; (iii) representation of principle directions 
on gear tooth surface £2 in 5/ t ; (iv) observation of the desired derivative (v) observation 

at the mean contact point of the desired direction of the tangent to the path contact on gear tooth 
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surface ; (vi) observation at the mean contact point of the desired length of the major axis of the 
contact ellipse; (vi) determination of principal directions and curvatures on pinion tooth surface 
Si at the mean contact point. 

Step 1: We set up a fixed coordinate system S h that is rigidly connected to the gear mesh 
housing (Fig.4.2.1(a)). In addition to S h, we will use coordinate systems 5 2 (Fig.4.2.1(a)) and S 1 
(Fig. 4. 2.1(b)) that are rigidly connected to gears 2 and 1, respectively. We designate with (f>' 2 and 
4>\ the angles of rotation of gears being in mesh. We have to emphasize that with this designation 
</>l(i = 1, 2) we differentiate the angle of gear rotation in meshing from the angle <(> l of gear rotation 
in the process of generation. 

The orientation of coordinate system S/ t is based on following considerations: (i) The axes of 
rotation of the pinion and the gear in a drive of spiral bevel gears intersect each other. Taking into 
account the possible gear misalignment, we will consider that the pinion-gear axes are crossed at 
angle F and the shortest distance is E. (ii) We will choose that Xh coincides with the pinion axis 
and Oh is located on the shortest distance (Fig.4.2.1(a)). (iii) Considering as given the shaft angle 
T, we will define jh~ the unit vector of Yh - as follows 

v % X O/l 

Jh = 7? ~ 

\%h X ah 

where 3h is the unit vector of gear axis that is parallel to plane ( Xh,Yh ). 

The coordinate transformation from S 2 to Sh is based on matrix equation 

= [M w ][Md2]r 2 (0G,^p) (4.2.2) 

where Sj (Fig.4.2.1) is an auxiliary fixed coordinate system. The unit normal to S 2 is repre- 


(4.2.1) 
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sented in S h as 


Here (Fig. 4.2.1) 


* 4 2) = [ Lhd\[Ld2}n2{0G,<t>p ) 


[M d2 \ 


10 0 0 
0 - cos <f > j sin <f > 2 0 

0 — sin <f >2 — cos <f> 2 0 

0 0 0 1 


[Mhd] = 


cosT 0 sinT 0 


0 E 


— sin T 0 cos r 0 


0 0 0 


(4.2.3) 


(4.2.4) 


(4.2.5) 


where T is the shaft angle. 

Equations (4.2.3), (4.2.2) and (4.2.3) enable to represent in Sh the position vector and unit 
contact normal at M by 


^h\ e G^* P ^2) 


4 2) (0G><&>^2) 


(4.2.6) 


where {0Q,<t>p) are the surface coordinates for the mean contact point at £2 ; the angle <f>' 2 of rotation 
of gear 2 will be determined from the equation of meshing (see below). 
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Step 2: The equation of meshing of pinion and gear at the mean contact point is 



4 12) = /(<&.*;. *2) = 0 


(4.2.7) 


Here (Fig. 4.2.1) 


4 I2) = [(4 1 ’ - 4 2) ) x 4 M) i - (e x 44 

4 2) = i- 1 0 of (rf'^i) 

4 2) = 4 [cos I’ 0 - sinr] r 

h n 2 

since at point M the angular velocity ratio is 

w ( 2 ) _ Ni 
uW ~ N 2 


(4.2.8) 


(4.2.9) 


(4.2.10) 


(4.2.11) 


Substituting equations (4.2.3), (4.2.8)- (4.2.11) in equation (4.2.7), we can solve equation (4.2.7) 
for <f > ' 2 . Usually equation (4.2.7) yields two solutions for <j>' 2 but the smaller one, say (cf > ' 2 )*, should 
be chosen. 

Step 3: We consider as known the principal curvatures and directions on X 2 a t any point of S 2 , 
including the mean contact point (see section 3). To represent in Sh the principal directions at the 
mean contact point, we use the matrix equation 

4 M) = [Lhd][L d 2 }a 2 (4.2.12) 
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where a 2 is the unit vector of principal directions on S 2 that is represented in ,5V The following 
steps of computational procedure are exactly the same that have been described in section 1.2. 
This procedure permit determination of the pinion principal directions and curvatures at the mean 
contact point. 
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5 Pinion Machine-Tool Settings 


5.1 Introduction 

We consider at this stage of investigation as known: 

(i) the common position vector and unit normal at the point of contact point M of £2 
and £i 

(ii) pinion surface principal directions and curvatures at M . 

The goal is to determine the settings of the pinion and the head- cutter that will satisfy the 
conditions of local synthesis. We consider that the pinion surface and the generating surface are 
in line contact. Henceforth, we will consider two types of the generating surface: (a) a cone 
surface, and (b) a surface of revolution. We consider that each side of the pinion tooth is generated 
separately and two head-cutters must be applied for the pinion generation. 

5.2 Head-Cutter Surface 
Cone Surface 

The cone surface is generated by straight blades being rotated about the zjr-axis (Fig. 5.2.1(a)). 
The £p equations are represented in coordinate system Sf that is rigidly connected to the head- 
cutter as following: 


r F = 


{Rep + sp sina^cos 6p 
{Rep + sjesin ap)s'mOp 
-Sp cos cup 


1 


(5.2.1) 


Here: sp and 6p are the surface coordinates; ap and R cp are the blade angle and the radius of the 
cone in plane zp = 0. The blade angle ap is standardized and is considered as known. Parameter 
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ap is considered as negative for the pinion convex side and ap is positive for the pinion concave 
side. The point radius R cp is considered as unknown and must be determined later. 

The unit normal to pinion tooth surface is represented as 


ftp = 


Np 


and 


Np 


dfp 

86 p 


x 


drp 

dsp 


( 5 . 2 . 2 ) 


i.e. 


) 


np = — [cos ap cosdp cosaps'mdp sinair] T 


( 5 . 2 . 3 ) 


The principal directions on the cone surface are: 


drp 

ef = = [- sin0p cos 6p 0] T 

W' 


drp 

JF) _ 

dsp 


[sin apcos6p sin ap sin 6 p 


— COS Ctp] 3 


( 5 . 2 . 4 ) 


( 5 . 2 . 5 ) 


The corresponding principal curvatures are 



cos ap 

Rcp + sp sina^ 


and 



( 5 . 2 . 6 ) 


Surface of Revolution 
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We consider that the head-cutter surface E F is generated by a circular arc of radius p by 
rotation about the 2o _ &xis that coincides with the 2 iP-axis of the head-cutter (Fig. 5.2.1(b)) and 
(Fig. 5. 2. 1(c)). The shape of the blade is represented in S 0 by the vector equations 


0 o N = 0 o C + CN = ( X { 0 c) + pcos\)i 0 + ( Z i c) + psm\)k 0 (5.2.7) 


Here: (Xq C \Zq C ^) are algebraic values that represent in S a the location of center C of the arc; 
p = \C N\ is the radius of the circular arc and is an algebraic value, p is positive when center C 
is on the positive side of the unit normal. ; A is the independent variable that determines the 
location of the current point N of the arc. By using the coordinate transformation from S 0 to Sf 
(Fig. 5. 2. 1(c)), we obtain the following equations of the surface of the head-cutter: 


r F = 


(X^ + p cos A) cos Of 
(X^ + p cos A) sintfjr 
+ p sin A 


1 


where A and Op are the surface coordinates (independent variables). 

The surface unit normal n F is represented by the following equations 


(5.2.8) 


njr — 


N f 

\n f \ 


, - dr F dr F 

“ d Nr = ~W r x -ex 


(5.2.9) 


Then we obtain 
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(5.2.10) 


np = — [cos A cos dp 


cos A sin#F 


sin \} T 


The variable A at the mean contact point M has the same value as the standardized blade angle 
ap. The principal directions on the head-cutter surface are 


dfp 

ef = - = [— sin#jr cos 6 p 0] r 

^ d6p ^ 


e II ~ 


drp 

~d\ 

I — I 

1 d\ 1 


[sin A cos dp sin A sin 6p — cos A] 3 


(5.2.11) 


(5.2.12) 


The principal curvatures are 


(F) _ 

I — 


cos A 
+ pcos A 


and 



(5.2.13) 


The radius R cp of the head-cutter in plane (Fig. 5.2.1) can be determined from the equations 


= x!;> + P 


1 


y( c ) 

i-(^y 


(5.2.14) 


5.3 Observation of a Common Normal at the Mean Contact Point for Surfaces E p , 
S 2 ) S p and Si 

We consider that at the mean contact point M four surfaces- E P ,E 2 ,E p and Si- must be in 
tangency. The contact of E p and S 2 at M has been already provided due to the satisfication of 
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their equation of meshing (3.1.10). Our goal is to determine the conditions for the coincidence 
at M of the unit normals to Ejp, E p and S 2 . The tangency of Si with the three above mentioned 
surfaces will be discussed below. 

We will consider the coincidence of the unit normals in coordinate system 5 m i . To determine the 
orientation of coordinate system 5/ t with respect to S m i , let us imagine that the set of coordinate 
systems Sh,S\ and S 2 (Fig. 4. 2.1) with gears 1 and 2 is installed in S m i with observation of following 
conditions (Fig. 5. 2. 2): (i) axis X} L of Sh coincides with axis x p of S p ; (ii) coordinate system Si 
coincides with S p and the orientation of Sh with respect to 5 1 is designated with angle 4>h = (<^i)° 
where 4>h is the to be determined instalment angle. Angle <j>h will be determined from the conditions 
of coincidence of the unit normals to Ejp, E 2 E p and E^ The procedure for derivation is as follows: 

Step 1: Consider that the coordinate system Sh with the point of tangency of surfaces E 2 and 
E p is installed in S m We may represent the surface unit normal in S m i by using the following 
matrix equation (Fig. 5.2.2). 




COS 71 

0 

— sin 7 i 


' 1 

0 

0 

0 

1 

0 


0 

COS <f>h 

— sin (f>h 

sin 71 

0 

cos 71 


0 

sin <f>h 

COS <f>h 

_ 




_ 


- 



(5.3.1) 


The unit vector has been represented by equation (4.2.3). 

Step 2: The unit vector to the surface of Ejr of the head-cutter that generates the pinion 
has been represented in Sp by equation (5.2.3) for a cone and equation (5.2.10) for a surface of 
revolution. Axes of coordinate systems Sf and S rn i have the same orientation and 


d F ) - 


n. 


ml 


= n\ 


(5.3.2) 
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Equations (5.3.1), (5.3.2), (5.2.3) and (5.2.10) yield the following equations 


cos Op 


n 


( 2 ) , . 

xl [ + sm a F sm 7 i 
cos 71 cos a p 


(5.3.3) 


COS 4 >! l 


( 2 ) , ( 2 ) 
<*1 %h + Q 2 n zh 


{n. 


S ) 2 + ("S ) 2 


sin <f)h = 


(2) 
“1 n yh 


( 2 ) 
“2 n zh 


(-S ) 2 + (»S ) 2 


(5.3.4) 


Here: 


ai = - cos ap sin#Jr a-i = cosaf-siny! cos#]? — sinaj? cos 71 (5.3.5) 

The advantage of the proposed approach is that the coincidence of the unit normals to surfaces 
T,p, T 2 , Up and £1 can be achieved with standard blade angles and without a tilt of the head-cutter. 

5.4 Basic Equations for Determination of Pinion Machine-Tool Settings 

At this stage of investigation we will consider as known: K / 1 an d • It is 

necessary to determine: K\ F \K^\cr( lF ),R cpt E Tn i,XBG> and m! Fl . Here: k ^\ and ($ ml 
are the principal curvatures and unit vectors of principal directions on the pinion surface that are 
taken at mean contact point M\ and are the position vector of M and the contact normal 
at M. The subscript “ml” indicates that the vectors are represented in S m i . Designations 
and indicate the principal curvatures of the surface of the pinion head-cutter that are taken 
at M , The angle <r^ 1F ) is formed by the unit vectors and e'f" ) of principal directions on £j and 
Tp ; R cp is the cutter ’’point radius” (Fig. 5. 5.1) that is measured in plane zp = 0 and is dependent 
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on k^\ E m i and Xqi are the pinion settings for its generation (Fig. 2. 1.1 and Fig.2.1.2); mp i, 
which is equal to and m Fl are the cutting ratio and its derivative. 

We recall that the pinion surface curvatures tt'p and have been determined in the process 
of local synthesis. Vectors e^, have been determined in system 5/ t . To represent these 
vectors in S m i we have to apply the coordinate transformation from Sh. to S m i similar to equation 
(5.3.1). 


®ral — [L T nip\{Lph]ci'h 


(5.4.1) 


where a represents that principal directions of the pinion surface and ej l j h , the position vector 
of mean contact point a rnl represents the corresponding vectors 6/^1 > anc ^ r ini ■ 

Now our goal, as it was mentioned above, is to determine k\ F \k^\c r( 1F \ E m i, Xai, mpi and 
m ' F r We recall that vectors e j ] and £$ are known from the local synthesis, and e\ F} and e'f I > 
become known from equations (5.2.4) and (5.2.5) for straight blade, and from equations (5.2.11) 
and (5.2.12) for curved blade, after the coincidence of the contact normal to surfaces S 2 , £ p and 
is provided. Thus parameter <r( lF ) can be determined from the equations 


sin 


?(M) /=0) 

' V/ml 


ml 


JF) \ 

X O 


cn „ -(IF) _ ,d F ) 

cos - e Iml e Iml 


(5.4.2) 


According to Fig. 5.5.1, since the Z^i-axis is parallel to Zp-axis, surface parameter sp for the 
cone surface at mean point can be determined as: 


Sp 


7 ml I 


cos ap 


(5.4.3) 
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(F) 

Parameter Kjj is equal to zero for a cone surface of the head-cutter and it must be chosen 
for a head-cutter with a surface of revolution. Then, the number of remaining parameters to-be 
determined becomes equal to five and they are: k^ F \ E m \,XGi,m\F and m' Fl . 

It will be shown below that we can derive only four equations for determination of the unknowns 
of the output data. Therefore one more parameter has to be chosen, and this is m^-the modified 
roll. Usually, it is sufficient to choose m Fl = 0, but the more general case with m' Fl ^ 0 is 
considered in this report as well. 

The to be derived equations are as follows, 


n ml V rr,^ — O 


ml 


(5.4.4) 


an&22 



(5.4.5) 


^ 11^23 = £12013 


(5.4.6) 


012 033 = 013023 


(5.4.7) 


Equation (5.4.4) is the equation of meshing of the pinion and head-cutter that is applied at 
the mean contact point. Equations from (5.4.5) to (5.4.7) come from the conditions of existence of 
instantaneous line contact between £1 and Ep- The coefficients a# in equation (5.4.5)-(5.4.7) are 
represented as follows, 

on = - K/^cos 2 cr( li? ) - Kpj sin 2 <7 (1F) (5.4.8) 


35 



(5.4.9) 


a\2 


a 21 


K 


(1) _ Jl) 


- K 


II 


sin 


2^ 1F ) 


a \z 


a 31 


K i v i 


u> 


(IF) n g) 


022 


= — K sin 2 <T^ lF ) — if.CV rr>K 2 /r( li? ) 


V/ 


77 


cos a K 


023 = o 3 2 = - u?( 1F ) n 


«33 = + 


n u> (lF) 


a<» «<'»] + [» m , a£> <£’] + (s ■ »i 2> ) 


Vectors in equation of meshing (5.4.4) can be represented as follows 


ofi = [ cos 7i 0 sin 7 i] r (|w (1) | = 1) 


u>. 


S~r-v> o u 

ICn 


^ap 


where i? ap , which is equal to , is the ratio of roll. 

m Fi 


d lF ) - .d 1 ) 


U>, 


ml 


(V. 


ml 


tv. 


(F) 

ml 


d 1F ) - d') _ d F ) 


ml 


v tr 


tr 


(5.4.10) 

(5.4.11) 

(5.4.12) 

(5.4.13) 

(5.4.14) 

(5.4.15) 

(5.4.16) 

(5.4.17) 
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(5.4.18) 


-4i) 


= U). 


;( x ) x J.M) 

mi * ml 


-Ymi sin 71 

X ml sin 7i - Z m i cos 71 

y ml COS 7! 


-(F) 

% 


u;. 


(F) 

ml 


X 


(r 


s(") 

ml 



^rnl^Fl “1“ -®ml^Fl 

Xm\m F \ + m F \X G \ cos 7a 
0 


(5.4.19) 


5.5 Determination of Cutter Point Radius 

Step.l: Equations (5.4.3), (5.4.6) and (5.4.7) yield the following expression for k'P 


( F ) _ 

I ~ 


4 1} 4V + 4/ w 


( F ), 


COS 2 <j ( 1F i 


+ 


4? - 4 1} sin 2 cr( lF ) - 4V 


4V sin 2 <r( ljr >) 
cos 2 cr(l F ) 


(5.5.1) 


Step 2: According to Meusnier’s theorem, the cutter radius R rn at the mean contact point is 
(Fig.5.5.1) 


cos a F 

= W^ 


(5.5.2) 


As shown in Fig.5.2.1. the cutter point radius can be determined for a straight blade cutter as 
follows, 


Rep — Rm s j? sin ot F 


(5.5.3) 


For the arc blade, the location of the center of the arc can be determined in S 0 by following 
equations, 
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= Rm - p cos a F 


(5.5.4) 


Zi c) = Rm - psina F 


(5.5.5) 


Knowing X^ and Zo C \ we can determine the point radius for the arc blade by equation (5.2.14). 
In order to find the position vector of the center of the head cutter, we define the following two 
vectors in S m i as shown in Fig. 5.5.1. 


= n cos ap — Cjj ^ sin ap 


(5.5.6) 


cos 9p 

— sin Of 

0 

0 ‘ 


' xi c) ' 


Xo° ) cos dp 

sin 9p 

cos Op 

0 

0 


0 


xi°^ sin 6p 

0 

0 

1 

0 


zi c) 


zi c) 

0 

0 

0 

1 


1 


1 



(5.5.7) 


where, p (°) is a unit vector directed from the blade tip M a to the cutter center Of, and p is 
a position vector directed from Of to the arc center C. Referring to Fig. 5. 2.2 and Fig. 5. 5.2, the 
position vector of the cutter center Of with respect to Oh, 7^°^ can be determined in system S m l 
as follows, 

For straight blade: 


J° F ) _ JM) 


ml 


S*F + Rcp^ o) 


(5.5.8) 
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For arc blade : 


4° F) = + P^rnl ~ ^ (5.5.9) 

It can be verified that the Z rn \ component of is zero, since equations (5.4.3) and (5.5.5) 
are observed. It is worth to mention that 0 F and s F are the surface coordinates where the contact 
is at mean point. The values of 6 F and s* F will serve as the initial guess in tooth contact analysis. 

5.6 Determination of rriF i = tv— , Emi and Xqi 

■* t ap 

The determination of cutting ratio R ap , settings E m i and Xgi is based on application of equations 
(5.4.2), (5.4.4) and (5.4.5). 

Initial Derivations 

It is obvious that equation of meshing (5.4.4) is satisfied at point M if the relative velocity E Fl } 
lies in plane that is tangent to the contacting surfaces at M . Thus, if velocity v^ Fl ^ satisfies the 
equation, 

* n ) = vW*P+vW$P (5.6.1) 

it means that equation of meshing (5.4.4) is also satisfied. Assuming that vectors of equation 
(5.6.1)) are represented in coordinate system S m i, we obtain 
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V 


(Fl ) (F) 

I e ImlX 


+ V 


(FI )(F) 

II e IIm\X 


d FX) - 

V ml ~ 


„(F1)(F) 


+ V 


(FI) (F) 

II e IImlY 


V 


(FI) (F) 

I e ImlZ 


+ V 


(Fl).(F) 

II e IImlZ J 


For further derivations we will use the following expressions for a 13 and 023 . 


013 = 4 F) 4 F1) + -^lUi + ^12 
“23 = K n^ v n^ + Muti + M22 , 


Here, 


^11 - n mlA'e /ml y - n m iYe Im IX 

M\2 = — COS 7l[ n ml y e /m\z - -«miz4miy] 

,, (F) (F) 

■““21 - n mixe HmlY - n miye// m ix 

M 22 = — COS 71 [“miy 4 /miZ - - n ^i^ e //iiy] 
ti = m^i — sin 71 


Using equations (5.6.2), (5.6.3) and (5.4.16), we obtain 


u 


(Fl) (F) 

/ e ImlZ 


+ V 


(Fl) (F) 

II e IImlZ 


+ Vml cos 71 = 0 


Following Derivations 


(5.6.2) 


(5.6.3) 


(5.6.4) 

(5.6.5) 


(5.6.6) 
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Step 1 : Expressions for v \ Fl ^ and v^\ Equations (5.6.6) and (5.4.4) represent a system of 
two linear equation in unknowns v ( ' ] Fl ' > and 4/ The solution of these equations for the unknowns 
yields: 


(Fl) 
Vj ’ 


— L 2 \t\ + L 22 - 


(5.6.7) 


v 


(Fl) 

II 


— L\\t\ + L\2 


(5.6.8) 


Here, 


j _ e JJmlz( a ll-^21 ~ a 12-^ll) 

fl 12 K I e IImlZ + a ll K H e ImlZ 


(5.6.9) 


_ e IlLlzi a ll M 22 - “ 12 -^ 12 ) - au^u’Ym! COS 7! 


L 2 2 = 


(FU 


ai24 F) 4/iiz + a n4/ } 4miz 


(5.6.10) 


'ii 


-e rmiz( a ll-^21 — ai 2 Mn) 

ai,« (F) e (F) + ai i K (F) e (F) 
“12 K I e IImlZ + “11 K II e ImlZ 


(5.6.11) 


r _ -4mlz( a H' M 22 - duMn) - an4 F) ^ml COS 71 


(5.6.12) 


Step 2 : Expression for 

Substituting the above equation in equation (5.6.2) , we obtain 
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Xuh + X 12 


= 


X2ltl + X 22 

X 3 lt\ + X32 


Here: 

A"n = L 2 i ej*2ix + L ^ e nLix 
X12 = L 2 ie\2\x + L ^ e nhix 
X 21 = Lneflr + l ^ € uLiy 

X 22 = L 22^J aY + L^ulnl Y 

A r 3i = L 21 e\2\z + L ^ e nliiz 

X32 = L22 ef 2 iz + L ^ e nLiz , 

Step 3 : Expression for i^T 1 ^ 

Equations (5.4.16) and (5.4.18) yield 

Xntj + 

^ 21^1 + -^23 
-^ 31^1 + X33 _ 

Here, 


4F) 

= 


X \3 = X\2 E m i sin 7 

X23 = X22 + X m \ sin 7 - Z ml cos 7 > 


A r 33 = A 32 + Yml cos 7 


(5.6.13) 


(5.6.14) 


(5.6.15) 


(5.6.16) 
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Step 4 : Expressions for triple products in equation (1.2.2) for a 3 3 


[n — E\\t\ + Eiiti -(■ E\z (5.6.17) 

where 

En = n m iyXu + n m ixX2i 

E 12 = n m iyXi2 - n m ixX 22 - n m izX2i COS 71 + n m iyA r 3i cos 7 j (5.6.18) 

ElZ = -(n m lzX22 — nmiyXzz) COS 71 

[iirf 1) $P} = Y n t 1 +Y 32 (5.6.19) 

where 


Y 21 = -n m ix X 21 sin 71 + n m iy(A"ii sin 71 - A 31 cos ) + n rnl zX 2 i cos 71 
E 22 = ^rniA '-^23 sin 7 i + n mi y (X 13 sin 71 - A33COS71) 


(5.6.20) 


[n w (r) v^] = Ynti + y 


12 


(5.6.21) 


where 


^11 = -n m iy{Xmi sin 7 ! - Z m 1 cos 71 ) - n m iyY m i sin 71 
Y u = sin 71 EH 


(5.6.22) 


Step 5 : Expression for the last term in equation (1.2.2) for a 33 . 

We have to differentiate between two derivatives: m 2 i and rn Fl . The first one, m' 21 , is apphed 
to provide a parabolic function of transmissions errors for the case of meshing of the generated 
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pinion and the gear. Such a function is very useful because it will allow to absorb linear functions 
of transmission errors caused by the gear misalignment. The other derivative, m' Fl , means that 
the cutting ratio in the process for pinion generation is not constant and it is just an additional 
parameter of machine-tool settings. 

In the approach proposed in this research project it is not required to have modified roll. 
However the use of such parameter in the more general case with m' F1 ^ 0 is also included to offer 
an extra choice. After some derivations, we obtain 

• iP) = Zntl + Z u ti + ^is (5-6.23) 

Here 

Z\\ = ( 2 C)(n m ixXn + + -^miz-^31) 

Z12 — {2C)[n m ixXi3 + n mlY X 23 + n m izX33 + sm^i^mixAn + N tji \yX2i + ^miz-^31)] ’ 
Z13 = (2C) sin 7i(n m ix-X’i3 + n m \Y X32 + n mizX33) 

(5.6.24) 

where 


2 C = 


(m F 1) 2 


(5.6.25) 


Step 6: Final expression for 033. 

Using the expressions received in steps 4 and 5, we obtain the following expression for (I33 


U33 — Z\t\ -f Z 2 t -f Z3 


(5.6.26) 
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where 


Z\ = k\ ^ L\\ + — Eu + Z\\ 

Zi = 2 k^ F ^ L21L22 + 2 kjj)Z/iiI<i 2 — E12 — Y21 + Y11 + Z12 
Z3 = L\2 + K-^iP L\ 2 — E13 — Y22 + Y12 + Z\2 

Step 7: New representations of coefficients (213 and 023- 
Equations (5.6.3), (5.6.7) and (5.6.8) yield 

dl3 = N 2 \t\ + N 22 

> 

®23 = -^11*1 + Ni 2 


Here; 

Nn — K’ffLu + M21 
N12 = L12 + M22 

N 21 = + Mu 

N22 = ^22 + M12 

Step 8 : Derivation of squared equation for <i 
Equations (5.6.30) and (5.4.5) yield 


ai tj diti -f 03 — 0 


(5.6.27) 

(5.6.28) 

(5.6.29) 


(5.6.30) 


(5.6.31) 


(5.6.32) 
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where 


ai = a 12 Z\ - N 2 \N\\ 

a 2 = o,\ 2 Z 2 — {N 2 \N12 + N 22 N\\) > 

a 3 = a-i 2 Zz ~ N 22 N\ 2 


(5.6.33) 


Solving equation (5.6.32), we obtain 


-a 2 ± \ al — 4aia 3 

** = ^ ( 5 ' 6 ' 34 > 

There are two solutions for t\ and we can choose one of them. If the tilt and the modified roll 
are not used, it can be proven that in this case ai becomes equal to zero and equation (5.6.32) 
yields 


h = (5.6.35) 

a 2 


knowing t \ , the ratio of roll may be easily determined as 


(5.6.36) 


According to equations (5.4.19) and (5.6.15), the blank offset and machine center to back can 
be determined by 


mpi = tj + smqi 
1 


R, 


■ap 


Eml 


Tmi mpi + X\\t\ + -X13 
m Fl 


(5.6.37) 
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(5.6.38) 


Xgx 


X 21^1 + X23 - ^Tnimpi 

mp 1 cos 7! 


Knowing E m 1 and Xqi , we may represent the position vector of the center of head-cutter with 
respect to the cradle center as follows, 


A°f) _ AOf) 


' ml 


= n" + 


X Gl cos 71 
-Eml 

Xqi sin 7! + X G1 


(5.6.39) 


In practice, the position of the center of the head cutter is defined by radial setting S T 1 and 
cradle angle q it which may be determined by the following equations, 


sn = yjix^y + iv^r } 

?i = sin(-^) 


(5.6.40) 


Since the cutter center Op must he in the machine plane, the component Z^ F ^ must be zero. 
Thus, the sliding base Xp\ may be determined as, 


X m = -Xqi sin 71 


(5.6.41) 


47 



6 Tooth Contact Analysis 


6.1 Introduction 

The tooth contact analysis (TCA) is directed at simulation of meshing and contact for misaligned 
gears and enables to determine the influence of errors of manufacturing , assembly and shaft deflec- 
tion. The basic equations for TCA are as follows: 

^{QfAfA i ) = ( 6 . 1 . 1 ) 

nu\0F,<f>F,<l>'i) = n^\e G ,(j) p ,(i) 2 ) ( 6 . 1 . 2 ) 

Equations (6.1.1) and (6.1.2) describe the continuous tangency of pinion and gear tooth surfaces 
Ei and E 2 . The subscript h indicates that the vectors are represented in fixed coordinate system 
Sh- The superscripts 1 and 2 indicate the pinion tooth surface E a and gear tooth surface E 2 , 
respectively. Vector equation (6.1.1) describes that the position vectors of a point on Ei and a 
point on E 2 coincide at the instantaneous point of contact M; vector equation (6.1.2) describes 
that the surface unit normals coincide at M. 

Parameters 6p and 4>f represent the surface coordinates for Ei; 6g and (f> p are the surface 
coordinates for S 2 . Parameters </>[ and <j>' 2 represent the angles of rotation of the pinion and gear 
being in mesh. 

Two vector equations (6.1.1) and (6.1.2) are equivalent to five independent scaler equations in 
six unknowns, which are represented as 


fi(0p, <S>f i 4> i) @Gi 0p, <^ 2 ) = o (i = 1,2, •■•,5) 


(6.1.3) 
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The continuous solution of equations (6.1.3) means determination of five functions of a param- 
eter chosen as the input one, say <f>' x . Such functions are: 

0f{4> i)> 4>f(4>i), 6g(<1> i), 1 ), </>2(<f>i), (6.1.4) 

In accordance with the theorem of Imphcit Function System Existence [4], solution (6.1.4) exists 
if at any iteration the following requirements are observed: 

(i) There is a set of parameters 

P{0f, <t>F , 0 Gt <p p , <f> 2 ) (6.1.5) 


that satisfies equations (6.1.1) and (6.1.2) 

(ii) The Jacobian that is taken with the above mentioned set of parameters and with (fi[ as an 
independent variable, differs from zero, i.e. 

# 0 ( 6 . 1 . 6 ) 

YFi “Gj Ypt r 2/ 

The solution of the system (6.1.3) of nonlinear equations is based on application of a subroutine, 
such as DNEQNF of the IMSL software package. The first guess for the starting the iteration process 
is based on the data that are provided by the local synthesis. 

The tooth contact analysis output data, functions (6.1.4), enable to determine the contact path 
on the tooth surface, the so called line of action, and the transmission errors. 

The contact path on pinion tooth surface is determined in Si by the following functions 
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n(o F) (j> F , $[) , Q F {<j > x ) , i ) ( 6 . 1 . 7 ) 


Similarly, the contact path on gear tooth surface is represented by functions 

M e G, <t > P , 4>i) , 6g(4>x) . 1 ) ( 6 . 1 . 8 ) 


Function ^ 2 (^ 1 ) re l ates the angles of rotation of the gear and the pinion being in mesh. Devi- 
ations of from the theoretical linear function represent the transmission errors (see section 

6.4). TCA is accomplished by the following procedure: (i) derivation of gear tooth surface, (ii) 
derivation of pinion tooth surface, (iii) determination of transmission errors, and (iv) determination 
of bearing contact as the set of instantaneous contact ellipses. 

6.2 Gear Tooth Surface 

The gear tooth surface E 2 and the surface unit normal have been represented in S 2 by equations 
(3.1.5) and (3.2.10), where 6q is the parameter of generating cone and <p p is the rotational angle 
of the cradle. Coordinate system S 2 is rigidly connected to the gear. To represent the gear tooth 
surface S 2 and its unit normal in fixed coordinate system Sh we can use the following matrix 
equations: 


4 2) (^G, 4>p, </>•>) = [M h 2{<l>2W2{6G,<j>p) 


(6.2.1) 


<t>p , </4) = [Lh2(<t>'2)}n2{°G, <j>p) 


( 6 . 2 . 2 ) 
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6.3 Pinion Tooth Surface 


We will consider two cases for generation of pinion tooth surface: (i) by a cone, and (ii) by a surface 
of revolution that is formed by rotating curved blades. 

Generation by a Cone Surface 

Step 1: We recall that the generating cone surface and the surface unit normal has been repre- 
sented in Sf by equations (5.2.1) and (5.2.3). 


r F = 


( R cp + s F sin a F ) cos 9p 
(R cp + sf sina.p) sin 0 F 
— sp cos ap 


1 


(6.3.1) 


rip - 


— cos ap cos 6p 

— cos ap sin Op 


— sp cos ap 


(6.3.2) 


where Sp and Op are the surface coordinates. 

Step 2: During the process for generation the cradle with the mounted cone surface performs a 
rotational motion about the Z mi -axis and a family of cone surfaces with parameter <j>p is generated 
in S m i . This family is represented in 5 m i by the matrix equation 


^ml(sF,0p,<j)F) — [M m ici(<£F)Kl(sF,0F) 


(6.3.3) 
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where 


r c i = fp + [ S T i cos q x - S r i sin 91 0 ] T 


(6.3.4) 


The position vector f c 1 represents a point of the cone surface in coordinate system S&i S T i and 
qi are the settings of the head-cutter center Of in SVni- 
Matrix [M m i c i] is (Fig. 2. 1.1) 


[Afmlcl] — 


cos <f>F sin <j>F 0 0 

-sin (j>p cos <f>F 0 0 

0 0 10 

0 0 0 1 


(6.3.5) 


The unit normal at a point of the generating surface is represented in Smi by 


i 4>f) = [ Lmlcl{<t>F)}n c i{0F ) 


(6.3.6) 


where n cl = tip- 

We recall that the generating cone surface is a ruled developed surface and the surface unit 
normal does not depend on sp (Parameter sp determines the location of a point on the cone 
generatrix.) Matrix [T m ici] is the 3x3 rotational part of [i TO ici] and is represented as follows, 


[^mlcl] — 


cos <f>p sin <f>p 0 

— sin^jr cos <f>p 0 

0 0 1 


(6.3.7) 
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Step 3 : Equation of meshing of the head-cutter cone with the pinion tooth surface. The equation 
of meshing is considered with vectors that are represented in S ml - Thus: 


n, 


(n . d lF ) 


ml 


v m\ 


(6.3.8) 


Here: is the sliding (relative) velocity represented as follows 


-(IF) / — (1) -(Fk - , s -(i) 

v rnl ~ ( w ml w ml ) r ml + Rml X U> m j 


(6.3.9) 


While deriving equation (6.3.9), we have taken into account that vector of angular velocity 
of pinion rotation does not pass through the origin 0 m l of S m i ; Rmi represents the position vector 
that is drawn from 0 rn l to a point of line of action of j Zmi can be represented as (Fig.2.1.2): 

Rmi = [A'gi cos 7i - Em\ A r Gisin7 1 ] T (6.3.10) 


vectors uK 1 ) and Q ( ' F '> are represented in S ml as follows 


4! = [cos7i 0 sin 71 f (|<2 {1) = 1|) 


(6.3.11) 


d F ) 


U>, 


ml 


R, 


-[0 0 if ( 


O' 


( F ) 


ap 


Ri 


ap 


U ) 


( 1 ) 


(6.3.12) 


Equations from (6.3.8) to (6.3.12) yield 
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(6.3.13) 


T\(6fAf) 

T 2 (0fAf) 


Here: 


Ti = X nm i(-£ m i sin 71 - Ai(sin7i - m F i)) + Y nm i(X B i cos 71 + ^(sin 71 - m Fl )) 

+ Z nml (E rn 1 COS71 + A x cos 7^ (6.3.14) 

T 2 = ^nmi(sin 71 - m F1 ) sina F sin(0 F + (£ F ) - Y„ m i[(sin7i - m F 1) sina F cos(0 F + <f) F ) 

- cos a F cos 7^ - Z nrnl cos 7! sina F sin(# F + <^ F ) (6.3.15) 


where 


A\ = sin(0 F + <£ F ) + 5 r i sin(-gi + 4>f) 
A 2 = R cp cos (Of + 4>f) + Sri cos (-qi + 4>f) 


(6.3.16) 


Step 4: Two-parametric representation of surface of action 

The surface of action is the set of instantaneous lines of contact between the generating cone 
surface and the pinion tooth surface that are represented in the fixed coordinate system S m \ ■ The 
surface of action is represented by equations (6.3.3) and (6.3.13) being considered simultaneously. 
These equations represent the surface of action by three related parameters. Taking into account 
equation (6.3.13) , we can eliminate s F and represent the surface of action in two-parametric form 

by 
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r m 1 = ^thI^fAf) 


(6.3.17) 


The common normal to contacting surfaces has been already represented in two-parametric 
form by equations (6.3.6). 

Generation by a Surface of Revolution 

Step 1: The shape of the blades is a circular arc (Fig. 5. 5.1) and such blades generate a surface 
of revolution by rotation about the head-cutter axis. 

The position- vector of the center of the generating arc is represented in S m l by the equation 


mi( e F,<l>F ) = [■Wmici3{p <c) + [ S r i cos q x - S rX smq 1 0] r } 


(6.3.18) 


where, 


[M, 


771 1 Cl I — 


cos <j>p sin <t>p 0 0 
-sin^jr cos <f>p 0 0 

0 0 10 


0 


0 0 1 


(6.3.19) 


and has been expressed by equation (5.5.7). 

Step 2: We will need for further transformations the following equations 


£/mi = [^rnlcljCjr/ = 


— sin(#F + <j>p) 
COS (6p + 4>f) 
0 


(6.3.20) 
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and 


— cos(0jr + <f>p ) 

T m i = - sin(0 F + <Pf) (6.3.21) 

0 

Here: is the unit vector of principal direction I on the head- cutter surface and f ml is a unit 

vector that is perpendicular to ej mi and the axis of the head-cutter (Fig. 6.3.1) 

Step 3 : To simplify the equation of meshing we will represent it by the following equation 

' 4T’ C) = 0 (6-3-22) 

where is the relative velocity of the center of the circular arc that generates the head-cutter 

surface of revolution. The proof that (6.3.22) is indeed the equation of meshing is based on the 
following considerations: 

(i) The relative velocity for a point of the head-cutter surface is represented by equation (6.3.9), 
given as 

Cf } = ("ml - "ml) X *Vnl + Rml X "ml (6.3.23) 

We can represent position vector f m \ for a point M as 

Fml = + pfirni (6.3.24) 
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where p is the radius of the arc blade. 

While deriving equation (6.3.24), we have taken into account that a normal to the head-cutter 
surface passes through the current arc center C; the sign of p depends on how the surface unit 
normal is directed with respect to the surface. 

Then, we may represent the equation of meshing as follows 




ml 


n, 


ml 


{(^ml - “2) X [(nm Hr pn„h + [Rmi X w^)]} • 

K^ml - “ml) X [( r ml + {Rml x £ (l) )] ■ nJJ 

4 1 f’ C) ' Ci' = 0 (6.3.25) 


Thus, equation (6.3.22) is proven. 

Step 4 : It follows from equation (6.3.22) that vector v^' C) belongs to a plane that is parallel 
to the tangent plane T to the head- cutter surface (Fig.6.3.2). This means that if vector is 

translated from point C to M it will he in plane T. The unit vector lies in plane T already. 
Then, we may represent the unit normal n ml by the equation 




ei mi X 

|e/ mi x v£f ,C) | 


(6.3.26) 


where 


is represented as follows, 




(6.3.27) 
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v m 1 


= O' 


( 1 ) 


+ [ ~ X G\ cos 7 ! E ml - X G1 sin Ti]^} 


(6.3.28) 


.T<1* C ) _ flO) _ 


"ml 


"ml 


ml 


(6.3.29) 


The advantage of vector equation (6.3.26) is that the surface unit normal at the point of contact 
is represented by a vector function of two parameters only, 9p and <pj?\ this vector function does 
not contain the surface parameter A. 

The order of co-factors in vector equation must provide that the direction of n m i is toward the 
axis of the head-cutter. The direction of n m \ can be checked with the dot product 

A = n ml • f mi (6.3.30) 


The surface unit normal has the desired direction if A > 0. In the case when A < 0, the desired 
direction of n m \ can be observed just by changing the order of co-factors in equation (6.3.26). 

To determine parameter A for the current point of contact we can use the equation, 


cos A — ■ T m i 


(6.3.31) 


Step 5: Our final goal is the determination in S m i of a position vector of a current point of 
contact of surfaces Ejr and Ei. This can be done by using the equation, 


r ml{@F > $f) — Pm\ P n ml 


(6.3.32) 
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where p is the radius of the circular arc. 

Finally, the pinion tooth surface may be determined in Si as the set of contact points. Thus: 

r x (6 F ,4> F ) = [Mi p ][M pmi ]f mi (9 f , 4> f ) (6.3.33) 


The unit normal to surface Si is determined in S i with the equation 


ni(0 F ,<t> F ) = [Li p ][Lp mi ]n mi (0 F ,(f) F ) (6.3.34) 


Here: f m i(0.F,^F) and n m i(0 F ,<f> F ) have been represented by equations (6.3.17) and (6.3.6) for 
straight blade cutter and by equations (6.3.26) and (6.3.32) for curved blade cutter. Here (Fig.2.1.2): 


[Af, 


p m i j — 


cos 7i 0 sin 71 

0 10 

- sin 71 0 cos 71 

0 0 0 


-X G i sin 71 

-®mi 

■(Xgi sin 71 + Xbi) 
1 


(6.3.35) 


[Mi P ] = 


10 0 0 

0 cos <f> 1 sin 4>i 0 

0 — sin <f> 1 cos (f> 1 0 

0 0 0 1 


(6.3.36) 
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where <f>i is the angle of the pinion rotation in the process for generation. Angles d>i and <f> F (the 
angle of rotation of the cradle) are related as follows: 

(i) in the case when the modified roll is not used and R ap is constant, we have 


<f>l — Rap4*F 


(6.3.37) 


(ii) when the modified roll is used, <j>i is represented by the Taylor’s series 


4>i = f(<f>F ) = Ra P (<t>F ~ C<j> 2 F - D<f> 3 F - E<t> 4 F - F<f> 5 F ) (6.3.38) 


where C, D, E and F are the coefficients of Taylor’s series of generation motion (see Appendix 
B). 

Step 7: The tooth contact analysis, as it was mentioned above, is based on conditions of tan- 
gency of the pinion and gear surfaces that are considered in the fixed coordinate system Sh (see 
section 6.1). To represent the pinion tooth surface and the surface unit normal in Sh we use the 
matrix equations 


= [M hl ]^\e F ^ F ) (6.3.39) 


nP = [L hl )ng](e F ,<l> F ) (6.3.40) 


Here: 
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1 0 


0 


0 


[M hl ] 


0 cos <f>\ - sin 4>\ 0 

0 sin<^ cos^j 0 

0 0 0 1 


where <j> [ is the angle of rotation of the pinion being in mesh with the gear. 

6.4 Determination of Transmission Errors 

The function of transmission errors is determined by the equation 


(6.3.41) 


H*\) = [4 - (< 4 ) 0 ] - 


JVj 


(6.4.1) 


Here: (<^) 0 ( i = 1,2) is the initial angle of gear rotation with which the contact of surfaces Si 
and £2 at the mean contact point is provided. Linear function 


f - rt)°l 


(6.4.2) 


provides the theoretical angle of gear rotation for a gear drive without misalignments. The 
range of <f>' 2 is determined as follows 


(4> 0 -£< A <(*>• + £ 


(6.4.3) 
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The function of transmission errors is usually a piecewise periodic function with period equal 

2tt 

t o fa = — (i = 1,2) (Fig.6.4.1). The purpose of synthesis for spiral bevel gears is to provide that 

™ i 

the function of transmission errors will be of a parabolic type and of a limited value S' (Fig.6.4.1). 

The tooth contact analysis enables to simulate the influence of errors of assembly of various 
types, particularly, when the center of the bearing contact is shifted in two orthogonal directions 
(see section 7). 

6.5 Simulation of Contact 

Mapping of Contact Path into a Two-Dimensional Space 

It was mentioned above that the contact path on the pinion and gear tooth surfaces is determined 
with functions (6.1.7) and (6.1.8), respectively. For the purpose of visualization , the contact path 
on the gear tooth surface is mapped onto plane ( X C ,Y C ) that is shown in Fig.6.5.1. The X c -axis 
is directed along the root cone generatrix and Y c is perpendicular to the root cone generatrix and 
passes through the mean contact point (Fig. 6.5.1). 

Consider that a current contact point N* is represented in 5 2 (Fig.6.5.2) by coordinates: 
-^ 2 (^ 2 ), RL'{4>' 2 ) where (j>' 2 is the angle of rotation of the gear and RL' = \EN\ = (Y 2 + Z 2 )% . 
Axis X 2 belongs to plane (X C ,F C ) (Fig.6.5.2). While mapping the contact path onto plane (X C ,Y C ), 
we will represent its current point N* by N that can be determined by coordinates X 2 and RL', 
where RL' = \EN\ = |.EiV*| (Fig.6.5.3). The coordinates of mean contact point M, XL and R, 
have been previously determined by equations (3.2.1) and (3.2.2). Drawing of Fig.6.5.3 yield 


O c N = 0 C 0 2 + 0 2 E + EN (6.5.1) 


Here: 
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(6.5.2) 


0 r 0 2 = 0 C K + K0 2 


K0 2 - -RL cos(jk - j2% (6.5.3) 

where 7* is determined by: 

RL 

7k = tan 1 (—) (6.5.4) 

Equations from (6.5.1) to (6.5.8) yield 


0 C K = —\0ji0 2 \ sin y 2 t c (6.5.5) 

0 2 E = X 2 cos 72 ? c - X 2 sin 72 (6.5.6) 

EN = RL‘ (sin ~) 2 i c + cos 7 2 / c ) (6.5.7) 

O c N = Xj c + Y c ] c (6.5.8) 


X c - X 2 (<j> 2 ) cos 72 + RL ( 4> 2 )sinj 2 - [(XL) 2 + (RL) 2 ] 2 cos(7fc -72) 

, . , , (6-5.9) 

Y e = X 2 (<f> 2 ) sin 72 + RL (<f> 2 ) COS72 - Zfisin72 

Contact Ellipse 
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Theoretically, the tooth surfaces of the pinion and the gear are in point contact. However, due 
to the elastic deformation of tooth surfaces their contact will be spread over an elliptic area. The 
dimensions and orientation of the instantaneous contact ellipse depend on the elastic approach 8 of 
the surfaces and the principal curvatures and the angel <r( 12 ) formed between principal directions 
efp and e'p of the surfaces. The elastic approach depends on the magnitude of the applied load. 
The value of 8 can be taken from experimental results and this will enable us to consider the 
determination of the instantaneous contact ellipse as a geometric problem. Usually, the magnitude 
8 is taken as 8 — .0.00025 inch. 

In our approach the curvatures and principal directions of the pinion and the gear are determined 
with the principal curvatures and directions of the generating tools and parameters of relative 
motion in the process for generation. 

Gear Tooth Principal Curvatures and Directions 

The procedure for determination of gear tooth principal curvatures and directions was de- 
scribed in section 1.2. Knowing functions 0 p (<f>' 2 ),(f> p (<f>' 2 ) from the TCA procedure of computa- 
tion, we are able to determine the position vector r m2 (6 p {<j> 2 ) , 4> p {(j>' 2 )) an d the surface unit normal 
f° r an instantaneous point of contact. The principal directions and curvatures 
for the generating surface can be determined from equations (5.2.4), (5.2.5) and (5.2.6). The pa- 
rameters of relative motions in the process for generation can be determined with equations (3.1.12) 
and (3.1.13). 

Pinion Tooth Principal Curvatures and Directions 

As it was mentioned above, the pinion tooth surface can be generated by a cone or by a surface 
of revolution. The derivation of principal curvatures and directions on the pinion tooth surface 
is based on relations between principal curvatures and directions between mutually enveloping 
surfaces of the head- cutter and Si of the pinion. The procedure of derivation is as follows: 

Step 1: We represent in S m l the principal directions on the head- cutter surface S^- using the 
following equations 
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O' = UI) 


(6.5.10) 


^ F \ - 
mu 


- f L le^ 

— L-°Tnx cj J e j Cl 


Step 2 : Parameters of relative motion in the process for pinion generation have been represented 
by equations (5.4.14) to (5.4.19). The derivative of cutting ratio, m' F1 , is equal to zero for the case 
when the modified roll is not used, and can be determined when the modified roll is applied as 
follows (see the Appendix) 


> _ d 2 <f>p ) 

mF1 ~ d<j>\ ~ [f{<t> F )f 


(6.5.11) 


where, 


/W) = Rap( 1 - 2 C4>f - 3 D<f> 2 F - 4 E4? f - 5F<f> F ) 


f"{<j> F ) = -R ap {2C + 6D(j> F + 12 E4 > 2 f + 20^4) 


(6.5.12) 


Step 3 : Now, since the principal curvatures and directions on Ef are known and the relative 
motion is also known, we can determine for each point of contact path the principal curvatures k / 
and kh of the pinion tooth surface Ex, the angle cr^ 1 ) and the principal directions on 

Ex- We use for this purpose equations (1.2.6) to (1.2.10). The principal directions on Ex can be 
represented in coordinate system Sh by the matrix equation (Fig. 5. 2. 2), 
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(6.5.13) 


Orient ation and Dimensions of the Instantaneous Contact Ellipse 
Knowing the principal directions and principal curvatures for the contacting surfaces at each 
point of contact path, we can determine the half-axes a and b of the contact ellipse and angle a* 1 ) 
of the ellipse orientation (Fig. 6. 5. 4). The procedure of computation is as follows [4]: 

Step 1: Determination of a and b 





- 1gig 2 cos 2 a + g% 




- K^ ] + \Jg\ - cos 2<r + g\ 




(6.5.14) 

(6.5.15) 

(6.5.16) 

(6.5.17) 


where, 


Kf = K ( p + K { /p 9i = K ( p - K { /p ( i = 1 , 2 ) 


(6.5.18) 


Step 2: Determination of (Fig.6.5.4) 
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(6.5.19) 


sin = [ n h 4ft 4? ] 


cos cr^ 12 ) = 4l ) ' 42 


7/i C Ih 
Ih e Ih 


Step 3 : Determination of at 1 ) 

Angle at 1 ) determines the orientation of the long axis of the contact ellipse with respect to 4/ 
(Fig. 6. 5. 4) and is one of the angles determined by the following equations, 


tan2at 1 ) 


<72 sin 2 crt 12 ) 

< 7 i - <72 cos 2 <r( 12 ) 


(6.5.20) 


Step 4 : The orientation of unit vectors ff and ( of long and short axes of the contact ellipse 
(Fig. 6. 5. 4) with respect to the pinion principal directions is determined with the equations 

fjh = 4V cos a^ 1 ) - 4 % sina^ 1 ) (6.5.21) 

4 = 4V sin + 4// cos (6.5.22) 

Step 5 : In order to visualize the contact ellipse we represent its axes of contact ellipse in plane 
(X C ,Y C ) (Fig. 6. 5.1), using the following equations 


m = [L2h]vh 


C*2 = [^2h]Ch 


(6.5.23) 


where 
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cos r o — sin r 


[£ 2 /i] = 0 — cos <(> 2 - sin <f> 2 0 10 

. 0 sin </>2 - cos <f> 2 J j_ sin T 0 cos T 


( 6 . 5 . 24 ) 


Axes of the contact ellipse form with the gear axes the following angles 


Ar, = arccos(7?2 • h) 
= arccos(£2 ■ *2) 


( 6 . 5 . 25 ) 


The unit vectors of axes of contact ellipse form in plane (X C ,Y C ) the following angles with the 
X c -axis (the generatrix of the root cone) : 


Ti = A v - 72 r 2 = - 72 


68 



7 V and H check 


The purpose of the so called V and H check is the computer aided simulation of the shift of 
the bearing contact to the toe and to the hill of the gear. The gear quality is judged with the 
sensitivity of the shape of the contact pattern and the change in the level of transmission errors to 
the above-mentioned shift of contact. 

7.1 Determination of V and H values 

Fig. 7. 1.1 shows the initial position M of contact point (it is the mean contact point), and the new 
position M* of the contact point). The shift of the contact pattern was caused by the deformation 
under the load. Coordinates XL and RL determines the location of M . For the following derivations 
we will use the following notations. 

(i) PF = A — A* is the shift of the center of bearing contact, where F is the tooth length 
measuring along the pitch line; p is an algebraic value, that is positive when A* < A and the shift 
is performed to the toe as shown in Fig. 7. 1.1. Usually, p is equal to 0.25. 

(ii) 8q and ac are the gear dedendum and addendum angle. 

(iii) PD = bo and P* D* = 8q are the gear dedendums that are measured in sections I and I*. 

(iv) h m — BD and h* = B* D* are the gear tooth heights. 

(v) r 2 is the pitch cone angle 

The determination of V and H for point contact M* is based on the following procedure. 

Step 1 : Determination of XL* and RL*. 

Fig. 7. 1.1 results in : 


h* = h m — pF(tanSG + tanac) 


(7.1.1) 
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b(} = ba — pF tan Sq 


(7.1.2) 


where b G = P*D* and ba = PD 


We assume that M*D* 


h* 


c , „ , „ h -f c 
— and MD = — - — 


where c is the clearance. 


Taking into account that 


0 2 M* = 0 2 P* + P*M * 


(7.1.3) 


we obtain 


XT* = A*cosT 2 + P*M* sinl? 2 = .A* cosr 2 + (6^- 


h* + C 
2 


) sinT 2 


(7.1.4) 


EX* = A* sin T 2 + P*M* cosr 2 = A*sinr 2 - (6£ - 


h* + C 
2 


) cos r 2 


(7.1.5) 


The surface coordinates ( Oq ,< j>p) can be determined by solving the following two equations, 


x 2 (8* g ,<i>;) = xl* 


(7.1.6) 


[Y2(eh,4>;)} 2 + [Z 2 (ehj;)} 2 = (RL*f 


(7.1.7) 
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Step 2: Determination of V and H 


We introduce the shift of the bearing contact in coordinate system Sh by V and H that are 
directed along the shortest distance between the pinion and gear axes, and the pinion axis, respec- 
tively (Fig.7.1.2). V is positive when the gear is shifted apart from the pinion in Y/, direction, H is 
positive when the pinion is withdrawn. It is obvious that 


T-* 

s* 

+ 

CM 

'TC- 

ll 

* 

CM 

Te: 

(7.1.8) 

[4 :) r = 4 1} + H Jh 

(7.1.9) 

Here: ( i = 1,2) is the position vector for the initial point of contact, (i = 1,2) is the 

position vector for the shifted contact point; ii u jh and k/ L are the unit vectors of coordinate axes 

s h . 

Equations of tangency at the new contact point provide 


(7.1.10) 


(7.1.11) 


Gear surface coordinates 6* G and <j>* can be determined from equations (7.1.6) and (7.1.7). 
Equations (7.1.10) and (7.1.11) yield 

v = {Y^\e^r p A 2)}* - (7.1.12) 
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(7.1.13) 


i)r - [Z l h\Pr*4M i)]* = 0 


(7.1.14) 


sin <f> x 


cos 4>i 


-n (2) n (1) + n (2) n (1) 
n hY n lZ + n hZ n lY 


(»ly)’ + Mi')’ 




n (2) n (1) + n (2) n (1) 
n hY n lY J_ n hZ n \Z 

(ng ) 2 + (nS ) 2 


(7.1.15) 


4a'( 0 G.C.^2)]* - n Wc(0F»0Fi0l)]* = 0 


(!), 


(7.1.16) 


Equations from (7.1.12) to (7.1.16) represent a system of five independent equations in six 
unknowns: V, H, </) 2 ,<I)[,0f and <f>p. The sixth independent equation, that is required for the 
solution of unknowns, can be derived based on the condition that the equation of meshing must be 
satisfied with the designed gear ratio, i.e., 


4° • 4 12) = WoAhMxfFitpMH) = o 


(7.1.17) 


In solving the above system, we first solve a sub-system composed of equations (7.1.14), (7.1.16) 
and (7.1.17) for <f>' 2 ,4>F and 9p, and then calculate the values of V and H directly, by equations 
(7.1.12), (7.1.13) and (7.1.15). 
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7.2 Tooth Contact Analysis for Gears with Shifted Center of Bearing Contact 


After the determination of parameters V and H, the tooth contact analysis for gears with shifted 
center of bearing contact can be performed similarly to the analysis described in sections 6.4. and 
6.5. The initial guess for the first iteration in the procedure of computations is provided by the set 
of six unknowns obtained in section 7.1. 
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Appendix A 


Generation with Modified Roll 


1 Introduction 

Modification of roll or sometimes called modified roll means that the cutting ratio is not constant 
but varied in the process for generation. The variable cutting ratio-the variable ratio of roll- can 
be provided by a cam mechanism of the transmission of the cutting machine or by the servo-motors 
of a computer controlled cutting machine. According to the developments of Gleason, the TCA 
program can analyze the process for generation up to members of the fifth order. However, due to 
the limitations caused by application of cam mechanisms only the parameters up to the third order 
are controllable in the process for generation. 

The modified roll is an additional parameter for the synthesis of spiral bevel gears. In our 
approach the synthesis of spiral bevel gears can be performed, as it was mentioned above, with a 
constant cutting ratio. However, we consider in this section the application of modified roll as well 
to provide a broader point of view on synthesis of spiral bevel gears. 

2 Taylor Series for the Function of Generation Motion 

According to the practice of Gleason, the kinematic relation between the angles of rotation of the 
workpiece and the cradle is represented by a Taylor’s series up to fifth order. To the knowledge of 
the authors, Gleason has never published any materials related to the kinematics of the modified 
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roll. However, Professor Zheng had done a good job in deciphering Gleason’s mechanisms for 
modified roll and represented the kinematic relations in his valuable book [5]. 

Consider that the angles of rotation of the pinion and the cradle are related by a nonlinear 
function 


4>i = f(<t>F) f eC K (K > 3) 


(A.l) 


We assume that <j>i = 0 at <pp = 0 and represent f(<f>F ) in the neighborhood of 4>p — 0 by the 
Taylor series as follows, 


<f>i — f ( 0 )4>f + / ( 0 )4>‘f + ■ • • 


(A.2) 


Taking into account that 


d<t> i 

d<f>F 


= f ( 4>f ) 


(A.3) 


We obtain 


/>) = 


U) 


(i) 




|0^.=:O — A 


up 


(A.4) 


where R ap is the ratio of roll. 

Without loosing generality of the solution, we can take u/ 1 ) = 1 and then obtain 
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(A.5) 


f (4 >f) 


d<j>F 

dt 


= 1 


Differentiation of equation (A.5) yields 


/'W)( 


d<pF 

dt 


f = f\<t>F) 


d 2 4>F 

dt 2 


(A.6) 


Equation (A.6) yields 


«2_ _ ) 
Up ~ f'(fo) 


(A.7) 


where 02 = 


d 2 <f>F 


is the angular acceleration of the cradle. 


Equation (A.7) with new designations can be represented as follows 


2 C 


a 2 


77, 


'/ ( 0 ) 


ap 


(A.8) 


Similar differentiation of higher order of equation (A. 3) yields: 


<pi = R ap {4>F - C4> 2 f - D<j> z F - E<f>p - F<f>p) 


(A.9) 


Here: 
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2C = -i/”(0) 

p 


6 - D = -Tr f '"(°) 

11 ap 


24£=--J-/ ,V (0) 

ap 


120F=--±-f V (0) 

ap 


Unfortunately, function /(</>f) cannot be represented in explicit form for certain cutting ma- 
chines, for instance, for the Gleason spiral bevel grinder. For such a case we will consider the 
following auxiliary expressions 

_ d3( t>F RrY - 03 

U 3 1,3 > bG -A 3 

at A wp 

(A.10) 

„ di< t>F 0 , nY a 4 

“4- *4 ' UDX ~ U * 

(A.ll) 

as = ~f , 120FX = 4 

at® 

(A.12) 

Then, differentiating equation (A. 6) and taking <^p = 0, we may obtain the following equations 

6 D = 6 CX - 3(2C*) 2 

(A. 13) 
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24 E = 24 DX + (2C)[15(2C) 2 - 10(6CX)] 


(A. 14) 


120F = 120EX - 15(2C)(24DX) + 105(2C) 2 [6CA - (2C) 2 ] - 10(6CX) 2 (A.15) 


The procedure for determination of coefficients C,D,E and F for the Taylor’s series (A. 9) when 
function f\{(f)F ) cannot be represented in explicit form is as follows: 

Step 1: Differentiate the implicit equation that relates <j>p and fa up to five times and then find 
u>p, 02 , a 3 , a 4 , a 5 in terms of fa and fa at fa = <f>p = 0. 

Step 2 : Considering fa = <f>p = 0, find 6CX ,24DX,120EX by equations (A. 10) - (A. 12). 

Step 3 : Find 2C,6D,24E,120F by equation (A. 8), (A. 13)- (A.15). 

3 Synthesis of Gleason’s Cam 
Introduction 

Gleason’s cam mechanism, as shown schematically in Fig. A. 3.1, is an ingenious invention that 
has been proposed and developed by the engineers of the Gleason Works. The mechanism trans- 
forms rotation of the cam about O q into rotation of the cradle about O c . The rotation of the cam 
about O q is related with the rotation of the pinion being generated, but the angles of cam rotation 
and pinion rotation, (f> q and fa, are related by a linear function when there is no cam settings. 

To authors’ knowledge, the engineers of the Gleason Works have not published the principles 
of synthesis and analysis of this mechanism. However, H. Cheng [6], Zheng [5] have made good 
contributions to the deciphering of this mechanism. The following is a systematic representation 
of synthesis and analysis of Gleason’s mechanism. 

The purpose of cam synthesis is to obtain the shape of the cam, considering that the angles of 
rotation of the cam and the cradle are related by a linear function, <f> c (fa). However, this function 
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can be modified into a nonlinear function by changing the location of the designed cam with respect 
to 0 q and the orientation of the cam guides that are installed on the cradle. Fig. A. 4.1 shows the 
settings of the cam mechanism with the designed shape : (i) the cam is translated along the line 
O c O q an amount AT ; (ii) and then, the cam guides are rotated about the cam rotation center 
and formed angle a with O c O q . It is obvious that the cam mechanism with the settings AT and 
a will transform rotation about O q to O c with a nonlinear function between the angles of rotation 
of the cam and the cradle. The deviation of this function from a linear one depends on settings of 
the cam mechanism and will be discussed in section A. 4. 

Coordinate Systems 

While considering the synthesis of the cam mechanism, we will use three coordinate systems: 
the movable coordinate systems S c and S q that are rigidly connected to the cradle and the cam, 
and Sf that is the fixed coordinate system (Fig. A. 3.2). 

Equation of Meshing, Contact Point in S c 

Assuming that the transformation of motion is performed with constant ratio of angular veloc- 
ities and in the same direction, we can determine the location of instantaneous center of rotation, 
I, in coordinate system Sj by using the equation (Fig. A. 3. 2) 


Wg _ E + r u 
u c r u 


(A.16) 


Where, E in the distance between the cradle center O c and the cam rotation center O q , r u is the 
so-called pitch radius of the cam. 

The location of instantaneous point of contact M on the guides can be determined by using 
the theorem of planar gearing [4]. According to this theorem the common normal to the guides 
and the cam at the point of their contact must pass through the instantaneous center of rotation 
I. Thus, contact point M and the unit normal at M are represented in S c as follows 
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r c — [-b - u 0 1 } T 

(A. 17) 

n c = [ 1 0 0 } T 

(A.18) 

(A.19) 

Here: b is an algebraic value ( b is positive if the left side of guides is considered and b is negative 
if the right side of guides is considered); u is a variable parameter that is determined with the 
equation 

u = (E + r u ) cos B c - E 

(A. 20) 

Equations (A. 17) and (A. 20) yield 


f c (0c) = [ -b /(0 C ) 0 1 ] r 

(A.21) 


where 

/(0c) = E-{E + r u ) cos 6 C (A. 22) 

Shape of the Cam 

The shape of the cam is a planar curve that is represented in S q by the matrix equation 
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f q (e c ) = [M qp ][M pf }[M fc ]f c (9 c ) (A. 23) 

Here: coordinate system S p is an auxiliary fixed coordinate system (Fig. A. 3. 2). Matrices in 
equation (A. 23) are represented as follows 


COS 9q 

— sin 9 q 

0 

0 

sin 6 q 

COS 9 q 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 


(A. 24) 


\M P f) = 


10 0 0 
0 10 E 
0 0 10 
0 0 0 1 


(A. 25) 


\M f J = 


cos e c 

sin 9 C 

0 

0 

- sin 9 C 

cos e c 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 


(A. 26) 


The normal to the cam shape is represented by the matrix equation 


^q{9c) — 


(A. 27) 
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Here: [L p f] is the identity matrix and is the (3 X 3) submatrix of the respective matrix [M]. We 
consider that the shape of the cam and its normal depend on the generalized parameter 9 q only 
since 


0 C = ( 


E + r u 




(A. 28) 


The final equations of the cam and its normal are represented as follows 


r 


g 


- 6cos (l^T^) + ( r »+ E ) co <^ 0 g) sin (wh d g) - Esm6 q 

- bsin (ETE: e «) “ ( r « + E ) cos (7Z+E e <i) cos (ET7Z 4 * * * * 9 g) ~Ecose q 

0 

1 


(A. 29) 


Tiq — 


COS ( E+r u Qq) 


(A. 30) 


4 Cam Analysis 

The cam analysis is directed at the determination of function 0*(6*) for a cam and guides with 

modified settings. The analysis is based on simulation of tangency of the designed cam with the 

cradle guides taking into account the settings of the cam and the guides. 

Coordinate Systems and Coordinate Transformation 

Coordinate systems Sj, S e and S c are rigidly connected to the guides and the cradle (Fig. 
A. 4. 1(a)). The guides after rotation about O c form angle a with the i/ c -axis. 
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Coordinate systems S q , S n and S m are rigidly connected to the cam. The settings of the cam 
with respect to S m are determined by A T and angle a. 

The cam and the cradle perform rotations about 0 m and Of, respectively (Fig. A. 4.2). The 
conditions of continuous tangency mean that the designed cam and the guides have a common 
normal and a common position vector at every instant in Sf. 

A current point N of the guide is determined in Sf with the equation (Fig. A.4.1 and Fig. 
A .4.4): 

r 4 f 1] = [M} c ]{M ce }[M ed )f d (A.31) 


where 


r d = [-b - A 0 1 f (A. 32) 


The unit normal is determined in Sf as follows 

= [£/c] [Lce][L e d]nd (A. 33) 


where 


n d =[l 0 0 ] T (A. 34) 


Here: 
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[M%] = 


cos 6* 

— sin 6* 

0 

0 

smO* 

cos 0* 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 


[M ce ] = 


10 0 0 
0 1 0 -E 
0 0 1 0 
0 0 0 1 


[M ed ] = 


cos a sin a 0 0 

— sin a cos a 0 0 

0 0 10 

0 0 0 1 


A current point of the cam and the unit normal at this point are represented in 
equations (Fig. A. 4. 1(b), Fig. A. 4.2). 


_ 




nf = [iy[L^,]K„][LyS, 


(A. 35) 


(A. 36) 


(A.37) 


Sj by the 


(A. 38) 


(A. 39) 
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Equations (A, 38) and (A. 23) yield 



4 2) = [ M / P ][^][CJ[M;p, p ][M p/ ][M /c ]f c (A.40) 

where 

r c = [ ~b - u 0 1 ] T (A. 41) 

Matrices [M qp ], [M p /j and [ M/ c ] have been represented by equations (A. 24), (A. 25) and (A. 26), 
respectively. Matrices [M/ p ], [M^], [M^J and [M* g ] are represented as follows (Fig. A.4.1(b), 
Fig. A. 4.2): 

' 1 0 0 0 

0 10 -E 

[M%] = 

0 0 10 
0 0 0 1 

cos 0* - sin 0* 0 0 

sin 0* cos 0* 0 0 

IK J = 

0 0 10 

0 0 0 1 
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[M ed ] = 


cos a sin a 0 0 

- sin a cos a 0 0 

0 0 10 

0 0 0 1 


[M, 


* 1 
nq J 


10 0 
0 1 0 
0 0 1 
0 0 0 


0 

-AT 

0 

1 


Equations (A. 39) and (A. 27) yield 




(A.44) 


(A. 45) 


(A. 46) 


Here: 


[n c ] = [l 0 0] 


(A. 47) 


Matrices [L*] and [X] are 3x3 submatrices of matrices [M*] and [M j. 

Equations of Tangency 

The tangency of cam and guides with modified settings is represented by equations 


^ 1] {6:,e c ,a) 


rf\9* q ,9 qt A T,a) 


(A.48) 
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n { p(0*,6 c ,a) = nf\0* q J q ,,a) 


(A. 49) 


We recall that vector equations (A.48) and (A. 49) yield a system of only three independent 
equations in four unknowns: 6*,6*,6 C and A; setting parameters AT and a are considered as given; 
0 q and 0 C are related with equation (A. 28) and 6 q is considered as a generalized parameter. Our 
goal is to determine the function that relates angles of rotation of the cam and the cradle, <f>* and 
</>*, and the parameters of settings a and AT, i.e. the function 

F(0*,0 c *,AT,a) = O (A. 50) 

Equality of Contact Normal-Satisfaction of Equation (A. 49) 

Equations (A.49), (A. 33) and (A.46) yield 

\L%)\A) = [T^][A][T qp ][T /c ] (A. 51) 


Then we obtain 


\B][A) = [A][c] 


(A. 52) 


Here: 
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[B] = [L%][L^] = 


cos(0* - 0* c ) sin(0* - 61) 0 

- sin(0* - 6*) cos (0* -6*) 0 


(A. 53) 


[C] = [£ W ][I/ C ] = 


cos(# 9 - 6 C ) - sin(6 g - # c ) 0 

sin(0 9 - 0 C ) cos (6 q — 6 C ) 0 

0 0 1 


(A. 54) 


[A] = [L,i\ = [£* ] = 


cos a sin a 0 


— sin a cos a 0 


(A. 55) 


Matrices (A. 52) are rotational matrices that describe rotation about axes of the same orienta- 
tion. This means that we can change the order of co-factor matrices and 


[B)[A] = [C)[A] 


(A. 56) 


This yields that 


[<?]- WJ = [A] 


(A.57) 


[cr i [B) = [i] 


(A. 58) 


where [I] is a unitary matrix and 
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[B] = [C] 


(A. 59) 


Equation (A. 59) yields 


1 


K = 


-{Oq - 0 C ) = 


£ + r u 9 


since 




E 

E + r u 




Equality of Position Vectors-Satisfaction of Equatio n (A. 48) 
Equations (A. 48), (A. 31) and (A. 38) yield 


[M;^-\M%]^{M%]{M ce ){M ed }v d = [MU[M^ q ][M qp ][M pf ][M fc }f c 


After transformations we obtain 


[Q}r<l = [S]f c 


Here: 


(A. 60) 


(A.61) 


(A. 62) 


(A.63) 
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[Q} = 


[S] = 


an 

012 

0 

014 

021 

022 

0 

024 

0 

0 

1 

0 

0 

0 

0 

1 

an 

a 12 

0 

“14 

021 

022 

0 

«24 

0 

0 

1 

0 

0 

0 

0 

1 


(A. 64) 


(A. 65) 


The rotational 3x3 submatrices of [£?] and [S] are equal due to the equality of contact normals 
(see equation (A. 52). The elements of [<5] and [.5] are represented by 


an = cos 77 ai 2 = sin 77 021 = — ai 2 a 22 : an 


(A. 66) 


where 77 = 0* - &* + a , and 

an = E[- sin(ff* - ff*) + sin 0*} , a 2 4 = -E[- cos((9* - 6*) + cos & q ] (A. 67) 


a ^4 = — I5[sin(0g — a) — AT sin a] , a^ = E[cos(0 g - a) — AT cos a] (A. 68) 

Matrix equation (A. 62) yields the following system of two linear equations 


012(7/ - A) + 014 - a* lA = 0 , a 22 (7/ - A) + 024 - a^ = 0 (A. 69) 
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Eliminating (u - A), we obtain 


(ai 4 - 0^4)022 — (a.24 - a 24 ) a n = 0 


(A.70) 


Equations (A.70), (A. 66), (A. 67) and (A. 68) results in 


r AT 

F(e* q J' c ,AT,a) = sin[-~(0* - <£)] + — sin(0* - 0* c ) - sin (0* - a) - sin a = 0 (A.71) 


Equation (A.71) represents in implicit form the displacement function for the cam mechanism 
with settings a and AT. It is easy to be verified, that equation (A.71) with AT = 0,a = 0 
represents the linear function , 


6 q _ 


E 4- r u 
r u 


6c 


(A.72) 


For Gleason’s grinder, E is equal to 15 inch. According to Gleason’s practice, the sense of 
rotation of the cradle and the cam is opposite to the assumption in the derivation in this report. 
Without loss of generality, by substituting 0* = -0* and 0* = -0* with E = 15 in equation (A.71) 
we obtain the final expression of the relation between 6* and 0* as follows, 

AT r 

sin(0* + a) - sina + — sin(0* - 0* q ) + sin^(0* - 6*) = 0 (A.73) 
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5 Determination of Coefficients of the Taylor’s Series 


The determination of the coefficients of the Taylor’s Series for generation motion with modified 
roll is a lengthy process. Gleason provides its customers with computer program which can select 
the cams with settings and analyze the effects of the modified roll. However Gleason’s program 
is a black box with no explanation for the determination of the coefficients of the Taylor’s Series. 
Valuable contribution to the understanding of Gleason’s program has been made by C.Q. Zheng 
[5). For reader’s convenience, a series of derivations are represented in this section, which coincide 
with the equations in [5] except some printing errors. 

In the process of generation, the cam rotates at a constant angular velocity . Without loss of 

de* 

generality, we assume that the cam rotates with unitary velocity, i.e. = 1. Using the procedure 

dt 

discussed in section A. 2, we differentiate equation (A. 73) five time as follows, 

AT v v 

u>* c cos(0* + a) + (a>* - 1)[— cos(0 c * - 0* q ) + ~ cos ^ sin(0 c * - 0* q )] = 0 (A.74) 


a 2 cos(0* + a) + ( u > *) 2 sin(0* + a) 

T T AT 

+az[ li cos ii (0 * " *!) + TjT cos( ^ - e « )] 

7* v A\ T 

-K - !)’[( jf) 2 sin j£(«S - flj)] + ^ sin(e: - flj)] = 0 (A.75) 


a 3 cos(0* + a) - 3a 2 (w*) 2 sin (6* + a) - (w*) cos(0* + a) 
t v AT 

+a 3 [^ cos ^(e: - e* q ) + ^ c OS (e; - 0 *)] 

v v AT 

-3 - IJKiJ ) 2 sin ^ - 01) + ^ sin(0* - 0*)] 
v v AT 

-K - ^(0 C * - 01) + ^ cos(0l - 0*)} = 0 (A. 76) 
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►O 3f 


a 4 cos (0* + a) - 4 a 3 u>* sin(0* + a) - 3 a\ sin(0* + a) 
-6(u>*) 2 a 2 cos (9* + a) + (u>*) 4 sin(0* + a) 

-4«>K - 1)[(^) J "in - »',) + H cos(9; - 9,*)] 

-3 <4iCf h ) 2 si- ^(*: -«;) + — sin(« 0 * - *;» 

v v f\ 7 1 

+“4iil “>• ifw - «;) + -jj- “>(«; - o',)] 

-««jK - 1 ) 2 Hj|) Scos - «,*) + co S («; - «*)j 
K* - l)*!!^) 4 sin - «J) + ^ sin(»; - «,*)) = 0 


a 5 cos(6»* + a) - [10a 3 (w*) 2 + 15 ajto* - ( u>*f } cos(0* + a) 

-[5 a 4 w* + 10a 2 a 3 - 10a 2 (w*) 3 ] sin(0* + a) 

7 * v AT 

+ “ 4 ii cos - e ;> + -jj- «»«: - *,*)] 

-5a 4 « - 1)[(^) ! sin - 9,*) + sin(9; - *;)] 

- 10 °sK - i) 2 [(^) s cos ^(9; - 9,*) + cos(9; - «;)] 

f /\T 

si" ,!(«• -9J)+ _ sin(9* - 9„*)] 

-15«tK - 1)1(^) 3 n°s jf (<E - «,*) + ^ cos (9* - 9J)J 

+10«j(“>: - 1) 3 K^) 4 ™ ^(9 C * - 9J) + ^ sin(9* - 9,*)] 

+K* - 1)'[(§) S oos y^(9* -ty+jf nos(9; - 9J)] = 0 


— 0 , we can determine w*, o 2 , a 3 , < 24 , as from above expressions as follows, 


u> 


k 

C 


r u + AT 

15 cos a + r u + AT 


(A. 77) 


(A. 78) 


(A.79) 
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0-1 


(A. 80) 


15 sin a 

15 cos a + r u + A T 


*\2 


K) 


a 3 = 


r 3 AT 

3a 2 a>* sin a + «) 3 cosa + (u>* - 1) 3 (^ + — ) 


cos a + 


r u + AT 
15 


AT, 


04 = 


6a 2 «) 2 cosa + [4a 3 u>* + 3a| - (w*) 4 ]sina + 6a 2 (u>* - l)^:^ + "^g") 


cos a + 


r u + AT 
15 


(A.81) 


(A. 82) 


®5 


cos a + 


r^ + A r{[ 1Qa3 (^) + 15a ^ ~ ( w c) ] COSQ 


15 


-|-{5a4a>* + 10d2 a 3 — 10a2(u>*) ] s i nc * 4* [10ci3 — 1) 

+15a|K"-l))4 + f)-K-mji + ^)} 


(A. 83) 


Using equation (A. 8) and (A. 10) - (A. 12), we obtain 


Rac — — - — 1 + 


o-'r 


15 cos a 
r u + AT 


(A. 84) 


2 C 


Ra 


Ra 


tana 


(A. 85) 


6CA = 


(1 — R l 3 r 3 

1 + 3(2C) tana + l , g . ^ (?g + ^T) 


1 + 


15 cos a v 15 3 
r„ + AT 

15 cos a 


(A. 86 ) 


2 ADX 


L_^{6(2C) cos a + [ 4 ( 6 CX) + 3 ( 2 C ) 2 

cos a + - — — 

15 

r 3 AT 

+6(2C)(1 -A ac ) 2 (^ + —)} 


1 ] sin a 


(A. 87 ) 
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120 EX = ^^{[lO(6CX) + 15(2C) 2 -l]cosa 

cos a + — 

15 

+[5(24 DX) + 10(2C)(6CX) - 10(2C)]sina 

+[10(6CX)(1 - R ac f + 15(2C) 2 (1 - £ ac ](^ + ^f) 
r 5 AT 

-(1 - Racfi^s + Jf)} (A. 88) 

Knowing M ac ,2C,6CX,24DX and 120 EX, we can determine QD,24E and 120.F by equations 
(A. 13) - (A. 15). 

6 Selection of Cams and Cam Settings 

In order to provide the desired low transmission errors and bearing contact, the ratio of roll R ap 
and second ratio of roll (2c), which are determined by the local synthesis, must be applied for the 
grinder. Due to the structure of Gleason’s grinder, the ratio of roll, R ap is related to R ac as follows, 


Rac — TTlil Tic Rap 


(A. 89) 


Here, R ac is the transmission ratio between the cam and the cradle, as determined by equation 
(A. 84). m c is a fixed gear ratio and is equal to 1 in Gleason’s grinder; m; is the gear ratio from the 
workpiece to the cam and is determined as, 


n 

mi = — 
Ui 


(A. 90) 
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here, n is the number of teeth of the workpiece and n; is the index internal, i.e. the gear tooth 
number skipped over in indexing. 

Prom equation (A. 80) and (A. 89), we obtain 


a 




Raa ~ 1 


(A.91) 


r u + AT = 


15 cos a 
Rac-l 


(A.92) 


The cams and their pitch radii r u are tabulized. A cam with pitch radius closest to (r„ + AT) 
calculated by equation (A.92) should be selected. After the cam with pitch radius r u is selected 
the corresponding setting, AT, can then be determined as: 


AT = 


15 cos a 
Rac-l 


- r u 


(A. 93) 


In some cases, it is also necessary to control 6 CX. In order to satisfy R ac ,2C and 6 CX, the 
value of rii can be used together with AT and a. Since 7i; must be an integral number it is difficult 
to obtain an accurate solution. But by careful selection of cams and index interval n,, a practical 
engineering solution is often achievable. 

When the cam and its settings are selected, it is then necessary to determine the coefficients of 
the Taylor’s Series of the generation motion and carry out the TCA to see how the higher order 
coefficients (i.e., 6 £>,241? and 120T) affect the transmission errors and bearing contact. If the result 
of TCA are satisfactory, then the gears can be ground by the selected cam and cam settings. 
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Appendix B 


Description of Program and Numerical Example 


Input and Output of Program 

The research project is complemented by a computer program, which can be used for the 
determination of machine tool settings through the method of local synthesis and simulate the 
transmission errors and bearing contact through TCA. The input data to the program include four 
parts. 

Part 1. Blank Data 
TNI : pinion number of teeth 
TN2 : gear number of teeth 
C : shaft offset ( zero for spiral bevel gear ) 

FW : width of gear 
GAMMA : shaft angle 
MCD : mean Cone distance 
RGMA1 : pinion root cone angle 
B1 : pinion spiral angle 
B2 : gear spiral angle 
RGMA2 : gear root cone angle 
FGMA2 : gear face cone angle 
PGMA2 : gear pitch cone angle 
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D2R : gear root cone apex beyond pitch apex 
D2F : gear face cone apex beyond pitch apex 
ADD 2 : gear mean addendum 
DED2 : gear mean dedendum 
WD : whole depth 
CC : clearance 

DEL : elastic approach (experiment datum) 

Part 2. Cutter Specifications 
RU2 : gear nominal cutter radius 
PW2 : point width of gear cutter 
ALP2 : blade angle of gear cutter 

Part 3. Parameters of Synthesis Condition 

FI21 : derivative of transmission ratio, negative for gear convex side and positive for gear concave 
side. The range is -0.008 < FI21 < 0.008 . 

KD : percentage of the half long axes of contact ellipse over face width. K D = 0.15 - 0.20 . 
ETAG : direction angle of contact path. For right hand gear, —80° < ET AG < 0° for gear convex 
side and —80° < ETAG < 0° for gear concave side; For left hand gear, 0° < ETAG < 80° for 
gear convex side and —80° < ETAG < 0 for gear concave side. When ETAG is close to zero, the 
contact path is along the tooth height, when the magnitude is increased, the contact path will have 
bias in and reach almost longitudinal direction if ETAG is close to 90 degrees. 

GAMAl : pinion machine root angle, which is the same as the pinion root angle if no tilt is used. 
RHO : radius of the arc blade if curved blade is used, which can be any values when curved blade 
is not used. 

C2 : second order ratio of roll if modified roll is used. C2 is zero without modified roll. 
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ALP1 : pinion cutter blade angle. ALPl is positive for gear convex side and negative for gear 
concave side. ALPl can be the same as ALP2. For better result, it is suggested for pinion concave 
side the magnitude of ALPl is smaller than ALP2 and for pinion convex side, the magnitude of 
ALPl is larger than ALP2.(As shown in the example). 

TN1I : number of teeth skipped over indexing. TN1I is only used in modified roll, the ratio between 
TN1I and TNI must not be an integer. 

Part 2. Control Codes 

JCL : JCL control V and H check, JCL — 1 means no V-H check. 

JCH : For right hand gear, set JCH = 1, for left hand gear set JCH = 2. 

JCC : For straight blade, set JCC = 1, for curved blade set JCC = 2. 

TL1, TL2 : Extra points on contact path, both should be less or equal than 2. 

The program output includes: (1) the machine-tool settings for gear and pinion; (2) the trans- 
mission error; (3) the contact path; (4) the length and orientation of the long axes of the contact 
ellipse; and bearing contact at toe and heel position. 

Numerical Example 

The model used in this report is the spiral bevel drive with the shaft angle of 90 degrees. In 
the numerical example, modified roll and curved blade for generation of gears were not used since 
favorable results were attained without them. The list of the blank data and machine tool settings 
are tabulized in the attached tables. 

The TCA results with V-H check are shown through Fig. B.l through Fig. B.6. The V and H 
values shown in the figures are of inch. It is shown also that the transmission errors are very 
small and the bearing contact is stable for both side at the three positions, toe, mean and heel. 
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BLANK DATA 



PINION 

GEAR 

NUMBER OF TEETH: 

11 

41 

PRESSURE ANGLE: 

20 ° 


SHAFT ANGLE: 

90 ° 


MEAN SPIRAL ANGLE: 

35 . 0 ° 


HAND OF SPIRAL: 

LF 

RH 

OUTER CONE DISTANCE: 


90.07 

FACE WIDTH: 


27.03 

WHOLE DEPTH: 

8.11 

8.11 

CLEARANCE: 

0.81 

0.81 

ADDENDUM: 

5.24 

2.061 

DEDENDUM: 

2.87 

6.05 

PITCH ANGLE: 

15 ° 1 ’ 

74°59 

ROOT ANGLE: 

13 ° 20 ' 

70°39 

FACE ANGLE: 

19 ° 21 ' 

76°40 


GEAR CUTTER 

SPECIFICATIONS 

BLADE ANGLE: 

20 ° 


CUTTER DIAMETER: 

152.40 


POINT WIDTH: 

2.79 
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GEAR MACHINE TOOL SETTINGS 


RADIAL SETTING (s): 

CRADLE ANGLE (q) : 

MACHINE CENTER TO BACK(X g ): 
SLIDING BASE(Xg): 

RATIO OF ROLL(R a ): 

BLANK OFFSET (E ffl ) : 

MACHINE ROOT ANGLE 0^): 


70.43577 
62.3981° 

0.00 
0.00 
1.032397 
0.0 
70.65° 

PINION MACHINE TOOL SETTINGS 


CONVEX CONCAVE 


CUTTER BLADE ANGLE: 

21.5° 

CUTTER POINT RADIUS: 

80.4876 

RADIAL SETTING (s): 

71.55166 

CRADLE ANGLE (q) : 

59.4638° 

MACHINE CENTER TO BACK(X g ): 

1.08497 

SLIDING BASE (X fi ) : 

-0.25021 

RATIO OF ROLL (R & ) : 

3.898097 

BLANK OFFSET (E m ): 

-2. 56862 (Up) 

MACHINE ROOT ANGLE (yJ : 

13.3333° 


18.5° 

71.7222 
69.04316 
64.0624° 
-1.58960 
0.36659 
3.788604 
2. 19033 (Down) 
13.3333° 
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Fig. 2.1.2 Pinion Generation: Additional Coordinate Systems 
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Fig. 2.3.2 Gear Generation: Instantaneous Axis of Rotation 
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Fig. 2.3.3 Gear Head-Cutter Installments 
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Fig. 3.1.1 Gear Head-Cutter 
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Fig. 3.1.2 Point Diameter 
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Fig. 5.5.1 Pinion Head-Cutter Blades 
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Fig. 6.3.1 Principal Directions of Pinion Head-Cutter Surface with 
Circular Arc Blades 
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Fig.6.3.2 Vis 
Ta T 
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Fig. 6.4.1 Parabolic Function of Transmission Errors 
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Fig. A.3.1 Cam Profiles and Guides 
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Fig. A.3.2 Schematic of Motions of Cam and Guides 
without Modified Settings 
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Transmission Error in Meshing Period 

Gear Convex Side - Mean 



TOOTH CONTACT PATTERN 
Gear Convex Side 


a . Mean V/H -0/0 

Bin 
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Root Cone Generatrix (mm) 


Fig. B.l Contact Pattern and Transmission Errors at Gear Mean 
Position of Convex Side 
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Transmission Error in Meshing Period 

Goar Convex Side Tee 
V - 9; H - 0.0 



TOOTH CONTACT PATTERN 
Goar Convex Side 



Root Cone Generatrix (mat) 


Fig. B.2 Contact Pattern and Transmission Errors at Gear Toe Position 
of Convex Side 
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Transmission Error in Meshing Period 

Geer Convex Side ~ Heel 
V ®“8 ; H ■ 3 



TOOTH CONTACT PATTERN 
Gear Convex Side 



B.3 Contact Pattern and Transmission Errors at Gear Heel Position 
of Convex Side 
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Transmission Error in Meshing Period 

Geer Concave Side MEAN 
v . 0: H - 0 



TOOTH CONTACT PATTERN 


Gear Concave Side 



Root Cone Generatrix Cm a) 


Fig. B.4 Contact Pattern and Transmission Errors at Gear Mean 
Position of Concave Side 
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Transmission Error in Meshing Period 

Gear Coaettve Sid# “ TOB 
V ■-11; H - 3 



TOOTH CONTACT PATTERN 



Fig. B.5 Contact Pattern and Transmission Errors at Gear Toe Position 
of Concave Side 
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Transmission Error in Meshing Period 

Gear Cascade Side - Heel 

V e 9 ; h o~3 



TOOTH CONTACT PATTERN 
Gear Concave Side 
Heel V/H = 9 /-3 



-!3.e 


-».» - 4 .# 0.0 4.0 1.0 

Root Cone Generatrix (no) 


11.0 


u.o 


> 0.0 


Fig. B.6 Contact Pattern and Transmission Errors at Gear Heel 
Position of Concave Side 
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this program is to derive the machine tool settings 

FOR PINION GENERATION & TEST THE RESULTS 

IMPLICIT REAL*8(A-H,0-Z) 

REAL *8 KS , KQ , K1 1 , K1 II , K2I ,K2II ,KFI , KFII , KD, KF, KH, mcd 
REAL *8 M11,M12,M13,L11,L12,L13,L14, M2 1 , M22 , M23 , L21 , L22 , L23 , L24 , 
&N11 ,N12,N21 ,N22 
real*8 xi (5) , x (5) , f (5) 

EXTERNAL FCNl , FCN2, FCN, FCNM , FCNR , FCNMR 

DIMENSION CH (3) , P (3) , ElEF (3) .ESN (3) ,EQN(3) ,W1VT2(3) ,WV12(3) , 
$W2VT1 (3) , EFIH (3) ,EFIIH(3) ,RH(3) , GNH (3) , E2IH (3) ,E2IIH(3) ,PI2P(20) , 
&E1IH(3) , El IIH (3) , EFI (3) ,EFII(3) ,E1I (3) ,E1II (3) , GN (3) , EFE1 (3) , 

&ERR (20) , xcp (20) , yep (20) ,AX1 (20) ,AX2(20) , ANG1 (20) ,ANG2 (20) 
COMMON/Al/CNST, TNI, TN2.C.FW, GAMMA, xl,rl ,mcd 
COMMON/A2/B 1 , RGMAl , FGMA1 , PGMAl , DIR , Dl F , AbDl , DEDl 
COMMON/A3/B2 , RGMA2 , FGMA2 , PGMA2 , D2R.D2F , ADD2 , DED2 , WD , CC , D2P 
COMMON/A4/ SR2 , Q2 , RC2 , PW2 , XB2 , XG2 , EM 2 , GaMA2 , CR2 , ALP2 , PHI2 , PHI2P 
COMMON/ A5/ SG , XM , YM, ZM , XNM , YNM , ZNM , X2M , Y2M.Z2M, XN2M , YN2M, ZN2M , 
&XNH2 , YNH2 , ZNH2 , XH2 , YH2 , ZH2 

COMMON/A6/ES (3) ,EQ(3) ,CN(3) ,W1(3) ,W2(3) ,W12(3) ,VTl(3) ,VT2(3) , 

$V12 (3) , KS , KQ , KF , KH , EF (3) ,EH(3) ,SIGSF,PI21 
COMMON / A7 / SR 1 , Q 1 , R c £ , PW 1 , XB 1 , XG1 , EMI , GaMAl , CR 1 , ALP 1 , PHI 1 , PH 1 1 P 
COMMON/A8/Sf ,XM1 , YMl , ZM1 , XNMl , YNM1 , ZNM1 ,X1M, YlM.ZlM, 

&XN1M, YN1M, ZN1M.XNH1 , YNH1 , ZNH1 , XH1 , YH1 , ZH1 
COMMON/A9/PHI2PO , OX, OZ ,XO , ZO ,RHO , ALP ,V,H,CRlT,PCRir 
COMMON/ A10/K1I, Kill, K2I.K2II, DEL, E1IH,E1IIH,E2IH,E2IIH, GNH, 

&A2P , B2P , TAU1R , TAU2R , A2L, B2L 

COMMON/All /RAM, PSI1 , C2 , D6 , E24 , F120 , CX6 ,DX24 , EX120, RU1 , DELT , RUP , 
$RAl , CPF.DPF, EPF , FPF 
CNST=DARCOS (-1 . 0D00) / 180 . 0D00 
C 

C. . . 

C... INPUT THE CONTROL CODES 
C. . . 

C 

c IF V AND H CHECK IS NOT DESIRED, SET JCN = 1 

C DO NOT SET JCN TO BE 3 

C 

JCL=2 

C 

C... FOR RIGHT HAND GEAR JCH=1, FOR LEFT HAND GEAR JCH =2 
C 

JCH=1 

C 

C... FOR STRAIGHT BLADE JCC=1, FOR CURVED BLADE JCC=2 
C 

JCC=1 

c 

C TLl AND TL2 ARE NUMBER OF EXTRA POINT ON CONTACT PATH 
C WHICH SHOULD NOT BE LARGER THAN 2 

C 

TL1-1.0 
TL2=1 . 0 

C. . . 

C. . . INPUT BLANK DATA OF GEAR AND PINION 
C. . . 

TN1-11.0 
TN2=41 . 0 
C=0.0 
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FW=27 . 03 
GAMMA=90.0*CNST 
MCD=76 . 56 

RGMAl= 13 . 3333*CNST 
Bl=35 . 0*CNST 
B2=35 . 0*CNST 
RGMA2=70.6500*CNST 
FGMA2=76 . 6667*CNST 
PGMA2=74.9833*CNST 
D2R= 0.0 
D2F=0.0 
ADD2=2 . 06 
DED2=6 . 05 
WD= 8.11 
CC=0.81 

DEL=0. 00025*25. 4 

. .. INPUT NORMINAL RADIUS OF GEAR CUTTER AND POINT WIDTH, BLADE ANGLE 

RU2=152. 4000/2.0 
PW2=2 . 79 
ALP2=20.0*CNST 
DC2=2.0*RU2 


. .. INPUT THE SYNTHESIS CONDITION PARAMETERS AND PINION CUTTER BLADE 
ANGLE, ALP 1 (FOR GEAR CONVEX, ALP1>0, FOR GEAR CONCAVE ALP1<0) 

**** GEAR CONVEX SIDE 

FI21=-0.0008 
KD=0. 180 
ETAG=-65.0*CNST 
GAMAl= 13 . 3333*CNST 
RHO= 250.0 
C2= 0.00 

ALP1=18.500*CNST 

TN1I-8.0 

SGN=DS IN (ALP 1) /DABS (DS IN (ALP 1 ) ) 

KSIDE=0 


GOTO 1989 

**** GEAR CONCAVE SIDE 

1990 CONTINUE 

FI21= 0.0008 
KD=0. 180 
ETAG= 65. 0*CNST 
GAMAl=13 . 3333*CNST 
RH0= 200.0 
C2= 0.00 

ALP1— 21.50*CNST 
TNI 1=8 . 0 

C. . . 

SGN=DSIN (ALP 1) /DABS (DSIN (ALP 1) ) 

KSIDE=1 

jcl=2 
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. . . INPUT GEAR MACHINE TOOL SETTINGS 

1989 Q2= 52.6589*CNST 
SR2=3. 8872*25. 4 
XG2=0 . 0 

XB2=-0. 0333*25. 4 
CR2=. 9772974 
RAG=I . 0/CR2 
GAMA2=RGMA2 
EM2=0 . 0 

RC2=RU2-SGN*PW2/2 . 0 
ALP2= SGN*ALP2 


. . . CALCULATE GEAR MACHINE TOOL SETTINGS 

1989 hg=mcd*dcos (pgma2-rgma2) -ru2“ds in (b2) 
vg=ru2*dcos (b2) 
q2=datan (vg/hg) 
sr2=dsqrt (hg**2+vg**2) 
xg2=0. 0 
GAMA2=RGMA2 
xb2=d2r*ds in (gama2) 

EM2=0.0 

rag=dcos (pgma2-rgma2) /dsin(pgma2) 
cr2=l . 0/ rag 
RC2=RU2-SGN*PW2/2.0 
ALP2= SGN*ALP2 

. . . DELT IS THE CAM SETING 

DELT=0 . 0 
C 

C... DEFINE THE MEAN CONTACT POINT 
C 

v=o.ooo 

H=0.000 

FA=FGMA2-PGMA2 

RA=PGMA2-RGMA2 

HM=CC+WD-0 . 5*FW* (DTAN (FA) +DTAN (RA) ) 

DED2R=DED2-0 . 5*FW*DTAN (RA) 

XL=MCD*DCOS (PGMA2) + (DED2R-HM/2 . 0) *DSIN (PGMA2) 
RL=MCD*DSIN (PGMA2) - (DED2R-HM/2 . 0) *DCOS (PGMA2) 

C. . . 

AGL=DATAN (RL/XL) 

OX=-DSQRT (XL**2+RL**2) *DC0S (AGL-RGMA2) 

OY=-D2R*DSIN (RGMA2) 

C WRITE (9,11) OX,OY,XL,RL 
C 
C 
C 

C... FIND SURFACE COORDINATES OF THE MEAN CONTACT POINT 
C 

ERRREL=0. ID- 10 
N=2 

ITMAX=200 

IF (JCH.EQ.l) THEN 
Q2=-Q2 
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XI (1)=270.0*CNST+B2 
ELSE 

XI (1) =B2 
END IF 
XI (2) =0.0 

CALL DNEQNF (FCN1 , ERRREL , N , UMAX , XI , X , FNORM) 

TH=X (1) 

PH=X (2) 

ST=DSIN(TH) 

CT=DCOS (TH) 

SH=DSIN(PH) 

CS=DCOS (PH) 

SP=DSIN (ALP2) 

CP=DCOS (ALP 2) 

SM=DSIN (GAMA2) 

CM=DCOS (GAMA2) 

THIG=TH 

C WRITE (9, 11) Xn2M, Yn2M,Zn2M,ZNM,YNM,ZNM 

C WRITE (9 , 11) XM , YM , ZM , s g , hm 
C 

C... DEFINE VECTORS TO COMPUTER THE SECOND ORDER PROPERTY OF GEAR 
C 

ES(1)— DSIN(TH-PH) 

ES (2) = DCOS(TH-PH) 

ES (3) = 0.0 

EQ(1)=-SP*DC0S(TH-PH) 

EQ (2) — SP*DSIN (TH-PH) 

EQ (3) =-CP 
CN (1) =XNM 
CN (2) =YNM 
CN (3) =ZNM 
KS=CP/(RC2-SG*SP) 

KQ=0 . 0 
Wl(l)— CM 
W1 (2) = 0.0 
W1 (3)=-SM 
W2 (1) = 0.0 
W2 (2) = 0.0 
W2(3)=-CR2 

VT1 (1) = YM*SM+EM2*SM 

VT1 (2)— XM*SM+ (ZM-XB2) *CM 

VT1 (3)=-YM*CM-EM2*CM 

VT2 (1) = YM*CR2 

VT2(2)=-XM*CR2 

VT2 (3) = 0.0 

DO 10 1=1,3 

W12 (I) =W1 (I) -W2 (i) 

V12 (i) =VT1 (i) -VT2 (I) 

10 CONTINUE 

. . . FIND THE PRINCIPAL DIRECTION AND CURVATURES AT MEAN POINT 
C 

PI21=0. 0 
CALL CURVAl 

C WRITE (9, 12) KF.KH.SIGSF 
C 12 FORMAT (3X,3(G14.7, 2X) ) 

K2I=KF 
K2I I=KH 
PHI2=PH/CR2 
sh2=ds in (phi2) 
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ch2=dcos (phi2) 

xX= CM*ef (l)+SM*ef (3) 

yY= ef (2) 

zZ=-SM*ef (1) +CM*ef (3) 
ef (1) =xx 

ef (2) = CH2*yY-SH2*zZ 
ef (3) = SH2*yY+CH2 5! zZ 
C WRITE(9,11) xx,yy,zz 

c 

xX= CM*eh (1) +SM*eh (3) 
yY= eh (2) 

zZ=-SM*eh (1) +CM*eh (3) 
eh (1) =xx 

eh (2) = CH2*yY-SH2 ,t zZ 
eh (3) = SH2*yY+CH2*zZ 
C WRITE (9,11) Ef (1) ,Ef (2) ,Ef (3) 

WRITE (9,11) Eh(l) ,Eh(2) ,Eh(3) 


ERRREL=0. ID-10 
N= 1 

ITMAX=200 
XI ( 1) =0 . 0 

CALL DNEQNF (FCN2 , ERRREL , N , ITMAX , XI , X , FNORM) 
PHI2P0=X (1) 

WR ITE (9 , 11) X (1) 

WRITE (9, ID XH2 , YH2 , ZH2 
WRITE (9 ,1D XNH2, YNH2, ZNH2 

CHP=DC0S(X(1)) 

SHP=DSIN(X(1)) 

CMH=DC0S (GAMMA) 

SMM=DS IN (GAMMA) 

XX= ef (1) 

YY=-ef (2) *CHP+ef (3) *shp 
ZZ=-ef (2) *SHP-ef (3) *chp 
EF (1) = XX*CMM+ZZ*SMM 
ef (2) = YY 

EF (3) =-XX*SMM+ZZ*CMM 
XX= eh (1) 

YY=-eh (2) *CHP+eh (3) *shp 
ZZ=-eh (2) *SHP-eh (3) *chp 
EH (1) = XX*CMM+ZZ*SMM 
eh (2) = YY 

EH (3) =-XX*SMM+ZZ*CMM 
WRITE(9, 11) EF (1) , EF (2) , EF (3) 

WRITE (9, ID EH (1) , EH (2) , EH (3) 

ETAG=90 . 0*CNST+S IGSF+ETAG 

LOCAL SYNTHESIS AT MEAN CONTACT POINT 


RH (1) =XH2 
RH (2) =YH2 
RH (3) =ZH2 
GNH (1) =XNH2 
GNH (2) =YNH2 
GNH (3) =ZNH2 
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E2IH (1) = EF (1) 

E2IH(2)= EF (2) 

E2IH (3) = EF (3) 

E2IIH (1) = EH (1) 

E2IIH(2)= EH (2) 

E2IIH (3) = EH (3) 

K2I=KF 

K2II=KH 

... RELATIVE MOTION PARAMETERS IN GEAR & PINION MESHING PROCESS 

R12=TN1/TN2 
Wl(l)— 1.0D00 
W1 (2) =0. 0D00 
W1 (3) =0. 0D00 
W2(l)= RI2*CMM 
W2 (2) =0.0000 
W2(3)=~R12*SMM 
W12 (1) =W1 (1) -W2 CD 
W12(2)=W1 (2)-W2 (2) 

W12 (3) =W1 (3) -W2 (3) 

VT1 (1) = O.ODOO 
VTl (2) = ZH2 
VT1 (3)=-YH2 

VT2 (1)= R12-'(YH2-C)*SMM 
VT2 (2) — R12* (XH2*SMM+ZH2*CMM) 

VT2 (3)= R12*(YH2-C)*CMM 
V12 (1) = VTl (1) -VT2 (1) 

V12 (2) = VTl (2) -VT2 (2) 

V12(3)= VTl (3) -VT2 (3) 

WRITE (9 , 3) V 1 2 ( 1) ,V12 (2) , V12 (3) 

3 FORMAT (5X, 3G14. 7 , /) 

... CALCULATE THE COEFFICIENT A13,A23,A33 

ESN (1) = GNH (2) *E2IH (3) -GNH (3) *E2IH (2) 

ESN (2) =- (GNH (1) *E2IH (3) -GNH (3) *E2IH (1) ) 

ESN (3) = GNH (1) *E2IH (2) -GNH (2) *E2IH (1) 

C. . . 

EQN (1) = GNH (2) *E2IIH (3) -GNH (3) *E2IIH (2) 

EQN (2) =- (GNH (1) *E2I IH (3) -GNH (3) *E2I IH ( 1) ) 

EQN (3) = GNH (1) *E21IH (2) -GNH (2) *E2IIH (1) 

C. . . 

W1VT2 (1) = W1 (2) *VT2 (3) -W1 (3) *VT2 (2) 

W1VT2 (2) =- (W1 (1) *VT2 (3) -Wl (3) *VT2 (1) ) 

W1VT2 (3) = Wl (1) *VT2(2)-W1 (2) ,v VT2(l) 

C. . . 

W2VT1 (1) = W2 (2) *VTl (3) -W2 (3) *VT1 (2) 

W2VT1 (2) =- (W2 (1) *VT 1 (3) -W2 (3) *VT1 (1) ) 

W2VT1 (3)= W2(1)*VT1(2)-W2(2)*VT1 (1) 

C. . . 

WV 12 (1) = W12 (2)*V12 (3)-Wl2(3)*V12(2) 

WV12 (2) =- (W12 (1) *V12(3) -W12 (3) ,v V 12(D) 

WV12 (3) = W12(1)*V12 (2) -W12 (2) *V12 (1) 

C... 

V12S=0. ODOO 
V12Q=0.0D00 
WNES=0 . ODOO 
WNEQ=0. ODOO 
VWN= O.ODOO 


147 



ooo ooooo 


W1TN=0. 0D00 
W2TN=0. ODOO 
VT2N=0.0D00 

C. . . 

DO 1 1=1,3 

V12S= V12 (i) *E2IH (I) +V12S 
V12Q= V12 (I) *E2IIH (I) +V12Q 
WNES= W12 (I) *ESN (I) +WNES 
WNEQ= W12 (I) *EQN (I) +WNE0 
VWN =GNH (I) *WV12 (i) +VWN 
W1TN=GNH (I) *W1VT2 (i) +W1TN 
W2TN=GNH(l)*W2VTl (l)+W2TN 
VT2N= GNH (I) *VT2 (i) +VT2N 
1 CONTINUE 

WRITE (9 , 6) V12S , V12Q 
6 FORMAT (5X.2G14. 7,/) 

COMPUTER THE COEFFICIENTS A13,A23,A33 

B13=-K2I*V12S-WNES 

B23=-K2II*V12Q-WNEQ 

B33=K2I ,v Vl2S**2+K2lI*V12Q**2-VWN-WlTN+W2TN+VT2N*Fl21*TN2/TNl 

... LOCAL SYNTHESIS OF MESHING AT MEAN CONTACT POINT 

DL=KD*FW 
SIGK2= K2I+K2II 
SIGG2= K2I-K2II 
A«DEL/DL**2 

Tl=-B13*V12Q+(B33+B13*V12S)*DTAN(ETAG) 

T2=B33+B23*(V12Q-V12S*DTAN(ETAG)) 

ETAP=DATAN (T1/T2) 

VS1 =B33/(B13+B23*DTAN(ETAP)) 

AMl=DSIN (ETAP ) **2 
AM2=-DSIN (2 . ODOO" ETAP) /2 . ODOO 

AN1= (B13-B23*DTAN (ETAP) )/((!. 0D00+DTAN (ETAP) **2) *VS1) 

AN2= (B 13*DTAN (ETAP) +B23) / ( (1 . 0D00+DTAN (ETAP) **2) *VS 1) 

SGN=DS IN (ALP 1 ) /DABS (DS IN (ALP 1) ) 

A=A*SGN 

SIGK= (4 . 0D00*A**2- ( AN 1 * * 2 + AN 2 * * 2 ) ) / ( 2 . 0D00*A- (ANl*DCOS (2 . ODOO* 
&ETAP) +AN2*DSIN (2 . 0D00*ETAP) ) ) 

SIGK1= SIGK2-SIGK 

Tl= 2. 0D00*AN2-SIGK*DSIN (2 . 0D00*ETAP) 
T2=SIGG2-2.0D00*ANl+SIGK*DCOS(2.0D00*ETAP) 

SIG12= . 50D00*DATAN (T1/T2) 

SIGG1= (2 . 0D00*AN2-SIGK*DSIN (2 . ODOO*ETAP) ) /DSIN (2. 0D00*SIG12) 
K1I= (SIGK1+SIGG1) /2 . ODOO 
K1II= (SIGK1-SIGG1) /2 . ODOO 
C. . 

C. . 

C. . 

C WRITE (9,11) ETAP, K1I, Kill 
C WRITE (9,11) SIGK,SIG12,SIGK1,SIGG1 
C WRITE (9, 8) T1,T2 
C 8 FORMAT (5X, 3G14 . 7) 

C. . 

C.. 

C.. . 

C... PRINCIPLE DIRECTIONS OF PINION SURFACE AT POINT M 
C. . . 
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DO 15 1=1,3 

El IH (I) = DCOS (SIG12) *E2IH (I) -DSIN (SIG12) *E2IIH (i) 

ElIIH (I) = DSIN (SIG12) *E2IH(I) +DCOS (SIG12) *E2IIH (I) 

15 CONTINUE 

C WRITE (9,11) E1IH(1) ,E1IH(2) ,E1IH(3) 

C WRITE (9,11) ElIIH (1) ,E1IIH(2) ,E1IIH(3) 

C. . . 

C... COINCIDE THE NORMALS OF CUTTER AND THE PINION SURFACES 
C. . . 

SMl=DSIN(GAMAl) 

CMl=DCOS (GAMAl) 

SP=DSIN(ALP1) 

CP=DC0S(ALP1) 

T1=-(XNH2+SP*SM1) 

T2= CP"CM1 
IF (JCH.EQ.l) THEN 
THF=DARCOS (T1/T2) 

ELSE 

THF=DARCOS(Tl/T2) 

THF=360.0*CNST-THF 
END IF 

BAl=~CP*DSIN (THF) 

BA2= CP*DSIN (GAMAl) *DCOS (THF) -SP ''DCOS (GAMAl) 

TT— (YNH2**2+ZNH2**2) 

CSH=- (BAl*YNH2+BA2*ZNH2) /TT 
SNH= (BA2*YNH2-BAl*ZNH2) /TT 
PHIH=2 . 0*DATAN2 (SNH, (1 . 0D00+CSH) ) 

C WRITE (9, 8) THF , PHIH 
C 

C... FIND THE PRINCIPAL DIRECTIONS OF PINION GENERATING SURFACE 
C 

EFI (1)=-DSIN(THF) 

EFI (2) = DCOS (THF) 

EFI (3) = O.ODOO 
EFII (1) = SP*DCOS (THF) 

EFII (2) = SP*DSIN(THF) 

EFII(3)=-CP 

C WRITE (9,11) EFI (1) , EFI (2) , EFI (3) 

C WRITE (9, 11) EFI 1(1) , EFI I (2) , EFII (3) 

. . . FIND THE PINION PRINCIPAL DIRECTIONS IN SYSTEM SMI 

XX= E1IH(1) 

YY= DCOS (PHIH) *E1IH (2) +DSIN (PHIH) *E1IH (3) 

ZZ=-DSIN (PHIH) *E1IH(2) +DCOS (PHIH) *ElIH (3) 

Ell (1) = CM1*XX~SM1*ZZ 
Ell (2) = YY 

Ell (3) = SM1*XX+CM1*ZZ 
XX= ElIIH (1) 

YY= DCOS (PHIH) *E1 I IH (2) +DSIN (PHIH) *E1 IIH (3) 

ZZ=-DSIN (PHIH) *E1IIH (2) +DCOS (PHIH) *E1IIH (3) 

El II (1) = CM1*XX-SM1*ZZ 
ElII (2)= YY 
E1II (3) = SM1*XX+CM1*ZZ 

. . , FIND THE UNIT NORNMAL IN SYSTEM SMI 

XX= XNH2 

YY= DCOS (PHIH) *YNH2+DSIN (PHIH) *ZNH2 
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ZZ=-DSIN (PHIH) *YNH2+DC0S (PHIH) *ZNH2 
GN (1) = CM1*XX-SM1*ZZ 
GN(2)= YY 

GN (3) = S M 1 * XX+ CM 1 * Z Z 
C 

C... EXPRESS THE POSITION VECTOR IN SMI 
C 

XX= XH2 

YY= DCOS (PHIH) *YH2+DS IN (PHIH) *ZH2 
ZZ=-DSIN (PHIH) *YH2+DCOS (PHIH) *ZH2 
RX= CMl*XX~SMl*ZZ 
RY= YY 

RZ= SMl*XX+CMl*ZZ 
C XX=~CP*DCOS (THF) 

C YY=-CP*DSIN(THF) 

C ZZ=-SP 

C. . . 

C WRITE (9,11) Ell(l) ,E1I(2) ,E1I (3) 

C WRITE (9,11) ElII(l) ,E1II(2) ,E1II(3) 

C WRITE (9,11) GN (1) , GN (2) , GN (3) 

C WRITE. (9 , 1 1) XX, YY,ZZ 

C WRITE (9,11) RX,RY , RZ 

C. . . 

DO 20 1=1,3 
XX=-EFI (I) 

EFI (I) = EFII(I) 

EFII (I) =XX 
20 CONTINUE 

C 

C WRITE (9,11) E1I(1) , Ell (2) ,E1I (3) 

C WRITE (9 , 11) E1II(1) ,E1II(2) ,E1II(3) 

C. . . 

C... FIND THE ANGLE FORMED BETWEEN PRINCIPAL CURVATURES 

C. . . 

EFE1 (1)= EFI (2) *E1I (3) -EFI (3) *E1I (2) 

EFE1 (2) =— EFI (1) *E1 I (3) +EFI (3) *E1 I (1) 

EFE1 (3) = EFI (1) *E1I (2) -EFI (2) *E1I (1) 

T1=0 . 0D00 
T2=0 . 0D00 
DO 30 1=1,3 
T2=EFI (I) *E1I (I) +T2 
T1=GN (I) *EFE1 (I) +T1 
30 CONTINUE 

S IGF 1=2 . 0*DATAN2 (Tl , 1 . 0+T2) 

C 

C... FIND THE CURVATURE OF PINION GENERATION SURFACE AT MEAN POINT 
C 

IF (JCC.EQ.l) THEN 
KFI=0. 0 

B12-0.5D00*(K1I-K1II)*DSIN(-2.0D00*SIGF1) 

Bl 1=KFI~K1I*DC0S (SIGF1) **2-KlII*DSIN (SIGFl) **2 
TKK= K1I*DSIN (SIGFl) **2+KlII*DCOS (SIGFl) **2 
KF 1 1 = ( B 1 2 * * 2+B 1 1 * TKK) / B 1 1 
C WRITE (9,11) SIGFl, KFII 

C 

C... FIND THE CUTTER POINT RADIUS AND ITS CENTER 
C 

C SF=SG*DCOS(ALP2) /DCOS (ALP 1) 

SF=DABS (RZ) /DCOS (ALP1) 

RCF=CP/DABS (KFII) -SF*SP 
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DO 40 1=1,3 

P (I) = GN (l) *CP~EFI (I) *SP 
40 CONTINUE 

RCX= RX~SF*EFI (l)+RCF*P(l) 

RCY= RY-SF*EFI(2)+RCF*P(2) 

RCZ= RZ-SF^EFI (3) +RCF*P (3) 

WRITE (9, 11) RCF.SF 
WRITE (9, 11) RCX,RCY,RCZ 
ELSE 

KFI=1 . 0/RHO 

B12=0.5DOO*(K.1I-K1II)*DSIN(-2.0DOO*SIGF1) 
B11=KFI~K1I*DC0S (SIGFl) **2~K1II*DSIN (SIGF1) **2 
TKK= Kl I *DS IN (S IGF 1) **2+Kl I I *DCOS (S IGF 1) **2 
KFI 1= (B 12**2+B 1 1 *TKK) /B 1 1 
WRITE (9,11) SIGFl, KFII 
DBT=-RZ 

RM=CP/DABS (KFII) 

ZO=-(DBT+RHO*SP) 

XO=RM-RHO*CP 

RCF=XO+RHO*DSQRT (1 . 0~ (ZO/RHO) **2) 

RCX=RHO*GN (1) -XO*DCOS (THF) +RX 
RCY=RHO*GN (2) -XO*DSIN (THF) +RY 
RCZ=RHO*GN (3) -ZO+RZ 
WRITE (9 , 777) XO.ZO 

'll FORMAT (3X,’ XO, XO = 1 , 2 (2X , G14 . 7) ) 

WRITE (9,11) RCF,RM,DBT 
WRITE (9, 11) RCX,RCY,RCZ 
END IF 


... THE FOLLOWING IS TO FIND THE CUTTING RATIO 

CSM1=CM1 

SNM1=SM1 

TlX-EFI(l) 

T1Y=EFI (2) 

T1Z=EFI (3) 

T2X=EFII (1) 

T2 Y=EFI I (2) 

T2Z=EFII (3) 

XN=GN (1) 

YN=GN (2) 

ZN=GN (3) 

RXC=RX 

RYC=RY 

RZC=RZ 

C 

C WRITE (9,11) RCX,RCY,RCZ 

C WRITE (9,11) T1X,T1Y,T1Z 

C WRITE (9,11) T2X,T2Y,T2Z 

C WRITE (9,11) XN, YN , ZN 

C WRITE (9,11) RXC,RYC,RZC 

C. . . 

C... THE FOLLOWING IS TO DETERMINE DELTA, EM, AND IFM 

C. . . 

M11=XN*T1Y-YN*T1X 
M12=-CSM1* (YN*T1Z-ZN*T1Y) 

M21=XN*T2Y-YN*T2X 
M22=-CSM1* (YN*T2Z-ZN*T2Y) 

C WRITE (9 , 1 1) Mil ,M12,M21 ,M22 


151 



o o o 


C L11=(B12/B11*M11-M21)/KFII 

C LI 2= (B 12/Bl 1*M12~M22) /KFI I 

C. . . 

C L21=-T2Z/T1Z*L11 

C L22=-T2Z/T1Z*L12-RYC*CSM1/T1Z 

DTT=B 1 2*KF I *T2Z+B 1 1 *KF I I*T 1 Z 
Lll— T1Z*(B11*M21-B12*M11)/DTT 

L12=(-B12*KFI*RYC*CSM1-T1Z*(B11*M22-B12*M12))/DTT 
L21=T2Z* (B1 l*M21-B12*Ml 1) /DTT 

L22= (~Bl 1*KFII*RYC*CSM1+T2Z* (B 1 1*M22-B12*M12) ) /DTT 
C WRITE (9,11) Ll 1 , L12 , L21 , L22 

Xll=L21*TlX+Lll*T2X 
X12=L22*T1X+L12*T2X 
X21=L21*TlY+Lll*T2Y 
X22=L22*T1Y+L12*T2Y 
X31=L21*T1Z+L11*T2Z 
X32=L22*T1Z+L12*T2Z 
E 1 1 = YN*X 1 1 -XN*X2 1 

E12-YN*X12-XN*X22-ZN*X21*CSM1+YN*X31*CSM1 

E13— (ZN*X22-YN*X32) *CSMl 

Y1 1=~XN* (RXC*SNM1-RZC*CSM1) -YN*RYC*SNMl 

Y12=Yll*SNMl 

X13=X12-RYC*SNM1 

X23=X22+RXC*SNM1~RZC*CSM1 

X33=X32+RYC*CSM1 

Y21— XN*X21*SNM1+ZN*X21*CSM1+YN*(X11*SNM1-X31*CSM1) 
Y22=~XN*X23*SNM1+YN*(SNM1*X13-CSM1*X33)+ZN*X23*CSM1 
C 

c. . . 

c 

c. . . 

C... THE EFFECT OF SECOND ORDER RATIO OF ROLL ON A33 
TM1 =XN *X 1 1 + YN*X2 1 +ZN*X3 1 
TM2=XN*X1 3+ YN*X23+ZN*X3 3 
C WRITE (9, 163) TM1.TM2 

C163 FORMAT (2X, 'TH1.TM2 ',2 (2X.G14. 7)) 

ZZ1=C2*TM1 

ZZ2=C2*(TM2+SNM1*TM1) 

ZZ3=C2*SNM1*TM2 

Z 1 =KF I I *L 1 1 **2-E 1 1 
Z2=2.0DOO*KFII*L11*L12-E12-Y21+Y11 
C Z3=KFII ,V L12**2-E13-Y22+Y12 

Z1-KFI*L21**2+KFII*L11**2-E11+ZZ1 

Z2=2.0*KFI*L21*L22+2.0D00*KFII*L11*L12-E12-Y21+Y11+ZZ2 

Z3=KFI*L22**2+KFII*L12**2-E13-Y22+Y12+ZZ3 

N 1 1 =KF I I *L 1 1 +M2 1 

N12-KFII*L12+M22 

N21-KFI*L21+M11 

N22=KFI*L22+M12 

C WRITE (9,11) B12,Z1,M11,N11 

AA=B12*Z1-N21*N1 1 
BB=B12*Z2-N21*N12-N22*N11 
CCC=B12*Z3-N22*N12 
C WRITE (9,11) AA, BB , CCC 

IF (AA.GT. 0.000001) GOTO 1949 

T1=-CCC/BB 

GOTO 1950 

1949 T 1= (-BB+DSQRT (BB**2-4 . 0D00 :V AA*CCC) ) / (2 . ODOO*AA) 

1950 FM1=T1+SNM1 
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CR1=FM1 
RAP=1 . 0/CR1 

VF3=X31*T1+X32+RYC*CSM1 

C. . . 

C. . . 

C. . . 

C... THE DETERMINATION OF EM AND DELTA 
C. . . 

EM1=(X11*T1+RYC*FM1+X13) /FM1 
XG1=(X21*T1-RXC*FM1+X23)/(FM1*CSM1) 

C. . . 

RCX=RCX+XG 1 *CSM 1 

RCY=RCY-EM1 

RCZ=RCZ+XG1*SNM1 

C. . . 

V1=RCY 

H1=RCX 

XB1=RCZ 

SRl=DSQRT (V1**2+H1**2) 

Q1=-DARSIN (Vl/SRl) 

XB1=-XB1 


DETERMINE THE CAM SETTING 


RAl=l . 0/CR1 

RAM-TN1/TN1I*RA1 

PSI 1=DATAN (C2*RAM/ (RAM- 1.0)) 

RL'P=15 . 0*DCOS (PSI 1) / (RAM- 1 . 0) 

RU1=RUP 

C DELT=0 . 0 
C 


CALL CAM 

C WRITE (9 , 191) RUP 

C WRITE (9 , 191) RU1 

C191 FORMAT (2X, 'RU1 DELT = ' , 2 (2X.G14. 7) ) 
C WRITE (9 , 199) RAl , C2 , D6 , E24 , F120 
C WRITE (9 , 199) RAl , CPF , DPF, EPF , FPF 
C 199 FORMAT (2X, 'CAM' ,2X,5(2X,G14.7)) 

C 


C. . . 

C 

C 

C 25 
C 

C44 

C 

C45 
C. . . 


131 


331 


WRITE (9 , 25) FM1 

FORMAT (5X, ' FM1 = ' ,G14.7) 

WRITE (9 ,44) EMI, SRI, Q1 

FORMAT (2X, ' EMI , SRI , Q1 1 , 3 (3X, G14 . 7) ) 

WRITE (9 , 45) XG1 , XB1 , VI , HI 

FORMAT (2X, 1 XG1 , XBl , Vl , HI' ,4 (2X, G14. 7) ) 


IF (KSIDE.EQ.0.0) THEN 
WRITE (9,131) 

FORMAT (/2X, ' * * * * * * * * * * * * * * * * * * •> * * * * * • - - ■ * * •< * * • - * * * * - • * * 1 f 
$ 2X,’* OUTPUT FOR GEAR CONVEX SIDE *',/ 

^ 2X, ’ * * * * * * * * * * * * * * * * * * * * *' f * * * * * * * * * * * '' * * * Vc v? ?v * 

ELSE 

WRITE(9,331) 

FORMAT (/2X, 1 ***************************************** 1 J 
$ 2X , 1 * OUTPUT FOR GEAR CONCAVE SIDE * * , / 

& 2X, 1 * * * * * * * * * * * * * * ?V V? * * Vr Vc 7 V Vr *,V * * Vc * * * Vc * * * * * * * * Vc t 


END IF 
WRITE (9, 13) 
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1 3 FORMAT ( / 2 X * ,c * ,f * * * * * « « Vc * * * 'V * Vc * Vc Vc Vc * * Vc Vc * Vc * * Vc Vc Vc Vc * Vc Vc * * * Vc » j 

$ 2X, 1 * gear cutter spicifications **,/ 

& 2X 1 * ,f ,f * * ,c * ' c * ,e * * * * * ,lf ,f * * * * ,c * * ’ * * * * * * >lf ,lf * * ******* ^ 

WRITE (9 , 1 15) DC2 , PW2 , ALP2 
115 FORMAT (/2X, 1 GEAR CUTTER DIAMETER 
$ 2X, ' CUTTER POINT WIDTH : 

S, 2X, ' CUTTER BLADE ANGLE : 

WRITE (9, 3) 

3 FORMAT (/2X 1 * * * * * * * * * ,c * * ,,f ,c * * ' c * * ~ r ~ c ,,f * ,,f ,,c ' ,f ,f ,c “ ~’ c ~' f ' r 1f ,lr ' lf ' c :f ^ ** f , / 
$ 2X,'* BASIC GEAR MACHINE-TOOL SETTINGS * ' , / 

& 2X * * * c ,c ,f ,c * * * * * * * “ * * ,f * ' r * * * * * * * * * * * * * * * * * * * * * * * * 1 j ) 


DC 2 =' , G14. 7 , / 
PW =’ , G14. 7 , / 
PHI2 =' , G14 . 7 , //) 


WRITE (9 

4) 

Q2 , SR2 , XG2 , XB2 , EM2 , GAMA2 , 

RAG 



FORMAT (/2X, 

1 BASIC CRADLE ANGLE : 

Q2 

= ' 

,G14.7 ,/ 

& 

2X, 

’ RADIAL SETTING : 

SR2 

= ' 

, G14 . 7 , / 

# 

2X, 

’ MACHINE CENTER TO BACK 

: XG2 

= ' 

, G14. 7 , / 

// 

2X, 

’ SLIDING BASE : 

XB2 

= 1 

, G14. 7 , / 

$ 

2X, 

1 BLANK OFFSET : 

EM2 

= ' 

,G14. 7 , / 

# 

2X, 

’ MACHINE ROOT ANGLE : 

GAMA2 

= ' 

, G14. 7 , / 

// 

2X, 

1 RATIO OF ROLL : 

RAG 

= ' 

, G14 . 7 , // 


WRITE (9, 6) 

FORMAT (/ / 2X , ' ****************************************** . f / 

& 2X , ' * BASIC PINION MACHINE-TOOL SETTINGS *',/ 

WRITE (9 , 7) ALP 1 , RCF , Q1 , SRI , XG1 ,XB 1 , EMI , GAMA1 ,RAP 


FORMAT(/2X,' BLADE ANGLE 


2X, ' POINT RADIUS : 

2X, ' BASIC CRADLE ANGLE : 
2X, ' RADIAL SETTING : 

2X, ' MACHINE CENTER TO BACK 
2X, ' SLIDING BASE : 

2X, ' BLANK OFFSET : 

2X, ’ MACHINE ROOT ANGLE : 
RATIO OF ROLL : 


ALP1 
RCF 
Q1 
SRI 
: XG1 
XB1 
EMI 
GAMAl 
RAP 


, G14 . 7 , / 

, G14. 7 , / 

, G14. 7 , / 

, G 1 4 . 7 , / 

, G14 . 7 , / 

, G14 . 7 , / 

, G14 . 7 , / 

, G14. 7 , / 

, G14 . 7 , //) 


61 

71 

1919 

16 

17 


2X, 

IF (JCC.EQ.2) THEN 
WRITE(9,61) 

FORMAT (//2X , ’ **************************************** 

& 2X, '* COORDINATES OF THE CENTER OF THE ARC 

& 2X, ' **************************************** 

WRITE (9,71) XO.ZO 
FORMAT (/2X, 1 RADIAL COORDINATE 
AXIAL COORDINATE 




XO «’ 
ZO =' 


, G14 . 7 , / 

, G14 . 7 , //) 


& 2X, 

ELSE 

GOTO 1919 
END IF 
CONTINUE 
WRITE (9, 16) 

FORMAT (/ /2X * 5f 5 * 5r ,r ?,f ;,f jr 5r ' f 5f 5f ,lf ;,f 5r 5r ' f ,f 5 ' r * * * ,f ,f * * '* ,c ,f ,,f ,f * c ,c * ,c ' f 5f ' c ,f ' r ,f ir ,c :: ,r ' r * ’ , / 

& 2X,'*CAM SETTINGS AND COEFFICIENTS OF TAYLOR SERIES*',/ 

& 2 X ' “ ~' f “ , c , c ' c * *' * ' f * “ * * Vc * * Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc » ^ 

WRITE (9,17) PSI1, RUP, DELT, RAl ,C2 ,D6 , E24 , F120 


FORMAT (/2X, ' 

GUIDE ANGLE : 

PSI1 

= ' 

, G14 . 7 , / 

& 

2X, ' 

CAM 

PITCH 

RADIUS : 

RUP 

= ' 

,G14. 7, / 

& 

2X, ' 

CAM 

SETTING : 

DELT 

= ' 

, G14 . 7 , / 

& 

2X, ' 

1ST 

ORDER 

COEFFICIENT 

RAl 

= ' 

, G14 . 7 , / 

& 

2X, ' 

2ND 

ORDER 

COEFFICIENT 

C2 

= ' 

, G14 . 7 , / 

& 

2X, ' 

3RD 

ORDER 

COEFFICIENT 

D6 

= ' 

, G14 . 7 , / 

& 

2X, ' 

4TH 

ORDER 

COEFFICIENT 

E24 

= 1 

, G14 . 7 , / 

& 

2X, ' 

5TH 

ORDER 

COEFFICIENT 

F120 

= 1 

, G14 . 7 , //) 


C 

C. 


CALL TCA 
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DEFINE THE INITIAL POINT 


XI (1) =THIG 
XI (2) = 0.000000 
XI (3) =THF 
XI (4) =0.0 
XI (5) = 0.00 

. . . FIND THE INITIAL CONTACT POINT 


5555 N=5 

ERRREL=0. ID-10 
ITMAX=200 
PHI2P=PHI2P0 
IF (JCC.EQ.l) THEN 

CALL DNEQNF (FCN, ERRREL , N , ITMAX.XI ,X, FNORM) 
ELSE 

CALL DNEQNF (FCNR , ERRREL , N , I TMAX , X I , X , FNORM) 
END IF 
PHI1P0=X(5) 


PHI2P1=PHI2P0-180.0*CNST/TN2-TL1*180.0*CNST/(6.0*TN2) 
PHI2P2=PHI2P0+ 180 . 0*CNST/TN2+TL2* 180 . 0*CNST/ (6 . 0*TN2) 
KK=1 

PHI2P=PHI2P1 
333 CONTINUE 

IF (JCC.EQ.l) THEN 

CALL DNEQNF (FCN , ERRREL , N , I TMAX . XI , X , FNORM) 

ELSE 

CALL DNEQNF (FCNR, ERRREL, N, ITMAX, XI , X , FNORM) 

END IF 
XI (1) =X (1) 

XI (2) =X (2) 

XI (3) =X (3) 

XI (4) =X (4) 

XI (5) =X (5) 

. .. find the transmission error 

ERRR=PHI2P-PHI2P0-TN1/TN2* (X (5) -PHI IPO) 
ERRR=PHI2P-PHI2P0+TN1/TN2* (X (5) -PHI1P0) 

ERR (KK) =3600 . 0*ERRR/CNST 
PI2P (KK)=PHI2P 

. .. computer the contact path 


xlc= x2m 

rlc= dsqrt (y2m**2+z2m**2) 

xcp(KK)= xlc*dcos (rgma2)+rlc*dsin(rgma2)+ox 
yep (KK) =-xlc*dsin (rgma2) +rlc*dcos (rgma2) +oy 

. .. COMPUTER THE PRINCIPAL DIRECTIONS AND CURVATURES OF GEAR 


TH=X(1) 
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PH=X (2) 

ST=DSIN(TH) 

CT=DCOS (TH) 

SH=DSIN (PH) 

CS=DCOS (PH) 

SP=DSIN(ALP2) 

CP=DCOS (ALP 2) 

SM=DSIN (GAMA2) 

CM=DCOS (GAMA2) 

. . . DEFINE VECTORS TO COMPUTER THE SECOND ORDER PROPERTY OF GEAR 

ES(1)=-DSIN(TH-PH) 

ES (2) = DCOS(TH-PH) 

ES (3) = 0.0 

EQ(1)=-SP*DC0S(TH-PH) 

EQ (2) =-SP*DSIN (TH-PH) 

EQ(3)=-CP 
CN (1) =XNM 
CN (2) =YNM 
CN (3) =ZNM 
KS=CP/(RC2-SG*SP) 

KQ=0 . 0 
W1 (1)=-CM 
W1 (2) = 0.0 
W1 (3)=-SM 
W2 (1) = 0.0 
W2 (2) = 0.0 
W2(3)=-CR2 

VT1 (1) = YM*SM+EM2*SM 

VT1 (2) =-XM*SM+ (ZM-XB2) *CM 

VT1 (3) = - YM*CM~ EM2 *CM 

VT2 (1) = YM*CR2 

VT2(2)=-XM*CR2 

VT2 (3) = 0.0 

DO 110 1=1,3 

W12 (I) =W1 (I)-W2(I) 

V12(l)=VTl(l) -VT2 (I) 

110 CONTINUE 


PI21=0. 0 

CALL CURVAl 

K2I=KF 

K2II=KH 

PHI2=PH/CR2 

sh2=ds in (phi2) 

ch2=dcos (phi2) 

xX= CM*ef (D+SM-'ef (3) 

yY= ef (2) 

zZ=-SM ! -ef (l)+CM*ef (3) 
ef (1) =xx 

ef (2) = CH2*yY-SH2*zZ 
ef (3) = SH2*yY+CH2*zZ 

xX= CM*eh(l)+SM*eh(3) 
yY= eh (2) 

zZ=-SM*eh(l)+CM*eh(3) 
eh(l) =xx 

eh (2) = CH2*yY-SH2*zZ 
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eh (3) = SH2*yY+CH2*zZ 

C. . . 

CHP=DC0S(PHI2P) 

SHP=DSIN (PHI2P) 

CMM=DCOS (GAMMA) 

SMM=DS IN (GAMMA) 

XX= ef (1) 

YY— ef (2) *CHP+ef (3) *shp 
ZZ=-ef (2) *SHP-ef (3) *chp 
E2IH (1) = XX*CMM+ZZ*SMM 
E2IH (2) = YY 

E2IH (3) =-XX*SMM+ZZ ,v CMM 

C.. . 

XX= eh (1) 

YY=-eh (2) *CHP+eh (3) *shp 
ZZ=-eh (2) *SHP-eh (3) *chp 
E2IIH(1)= XX*CMM+ZZ*SMM 
E2IIH (2) = YY 
E2IIH (3) =-XX*SMM+ZZ*CMM 


. .. COMPUTER THE PRINCIPAL DIRECTIONS AND CURVATURES OF PINION 

TH1=X (3) 

PH1=X (A) 

STP=DSIN (TH1+PH1) 

CTP=DCOS (TH1+PH1) 

IF(JCC.EQ.l) THEN 
SP1=DSIN (ALP1) 

CPl=DCOS (ALP1) 

ELSE 

SGN=ALP1/DABS(ALP1) 

ALP=SGN*ALP 
SP 1=DSIN (ALP) 

CP l=DCOS (ALP) 

END IF 

SM1=DSIN (GAMAl) 

CMl=DCOS (GAMAl) 

. .. DEFINE VECTORS TO COMPUTER THE SECOND ORDER PROPERTY OF PINION 

ES(1)=-STP 

ES (2) = CTP 

ES (3) = 0.0 

EQ (1) = SP1*CTP 

EQ (2) = SP1*STP 

EQ(3)=-CP1 

CN (1) =XNM1 

CN(2)=YNM1 

CN (3) =ZNM1 

IF (JCC.EQ.l) THEN 

KS=CP1/(RCF+SF*SP1) 

KQ=0.0 

ELSE 

KS=DCOS (ALP) / (RHO*DCOS (ALP) +XO) 

KQ=1 . O/RHO 
END IF 
W1 (1) = CM1 
W1 (2) = 0.0 
W1 (3) = SMI 
W2 (1) = 0.0 
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W2(2)= 0.0 
W2 (3) = CR1T 

VT1(1)=-YM1*SM1-EM1*SM1 

VT1 (2) = XM1*SM1-(ZM1-XB1)*CM1 

VT1 (3) = YMl*CMl+£Ml*CMl 

VT2 (1) =-YMl*CRlT 

VT2 (2) = XM1*CR1T 

VT2 (3) = 0.0 

DO 210 1=1,3 

W12 (I) =Wl (I)-W2(I) 

V12 (I) =VTl (I) -VT2 (I) 

210 CONTINUE 
C 
C 

PI21=PCR1T 
CALL CURVAl 

C WRITE(9, 12) KF,KH,SIGSF 

K1I=KF 
K1II=KH 

C PHI1=PH1/CR1 

SH1=DSIN (PHI1) 

CHl=DCOS (PHI 1) 

XX= CMP'EF (1) +SM1*EF (3) 
yY= ef (2) 

ZZ=-SM1*EF (1) +CM1*EF (3) 
ef (1) =xx 

EF (2) = CHP'YY+SHP'ZZ 
EF (3) =-SHl* YY+CHP-ZZ 
c 

XX= CMP'EH (1) +SM1*EH (3) 
yY= eh (2) 

ZZ=-SM1*EH (1) +CMP-EH (3) 
eh (1) =xx 

EH (2) = CH1*YY+SH1*ZZ 
EH(3)="SH1*YY+CH1*ZZ 

C. . . 

CH1P=DC0S (X (5) ) 

SH1P=DSIN (X (5) ) 

El IH (1) =EF (1) 

ElIH (2) = CH1P*EF (2) -SH1P*EF (3) 

E1IH (3) = SH1P*EF (2) +CH1P*EF (3) 

E1IIH (1) =EH (1) 

El IIH (2) = CH1P*EH (2) -SH1P*EH (3) 

E1IIH (3) = SH1P*EH (2) +CH1P*EH (3) 

DO 109 1=1,3 
E1IH(I)=-E1IH(I) 

El IIH (I) =~El IIH (I) 

109 CONTINUE 

. . . COMPUTER THE DIMENSION AND ORIENTATION OF THE CONTACT ELLIPSE 

GNH (1) =XNH2 
GNH (2) =YNH2 
GNH (3) =ZNH2 
CALL ELLIP 
AXl (KK) =A2L 
AX2 (KK)=B2L 
ANG1 (KK) =TAU1R 
ANG2 (KK) =TAU2R 
C 
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c. . . 
c. . . 

441 


C 


555 

444 

C 

551 


747 

666 

C 

661 


889 

888 

C 

C 

C 


C 

C. . . 

C 


c 

c. . . 

c 


KK=KK+1 

PHI2P=PHI2P+180.0*CNST/(TN2*6.0) 
IF(PHI2P.LE. (PHI2P2+0. 0001) ) GOTO 333 


WRITE (9 , 441) 

FORMAT ( / , ' * * * * * * * * * * ' c ' c * * * * * * * * * * * * * * * * * * * * * * * * V? * Vc Vc * Vc * * Vc * * Vc * Vc • 

& /,'* TRANSMISSION ERROR IN A MESHING PERIOD 

cj j i Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc * * Vc V? Vc * Vc Vc Vc Vc Vc Vc Vc V? Vc * Vc Vc V; * Vc Vc Vc * Vc Vc Vc Vc * Vc * * Vc it Vc Vc » 


DO 444 1=1 ,KK~1 
PI2P(I)=PI2P(I)/CNST 
WRITE (9 , 555) PI2P (I) , ERR (I) 
FORMAT (3X , 3 (G1 4 . 7 , 3X) ) 
CONTINUE 


WRITE (9 , 551) 

F 0 R MAT ( / * * * * * '' ,f ,f * * * * * * * * * ,f ” ” * * * * * * * * “ * * Vc * * Vc Vc * * * * Vc Vc Vc Vc Vc V? Vc Vc Vc Vc Vc ' 

&, / , ' * CONTACT PATH FOR A PAIR OF TEETH IN MESH * ' , 

DO 666 1=1 ,KK-1 

WRITE (9, 747) XCP (I) , YCP (I) 

FORMAT (3X, 2 (G14. 7 , 3X) ) 

CONTINUE 


WRITE (9 , 661) 
FORMAT (/ , ' ** 
S* / , ' " 

$ /,'** 


DIMENSION AND ORIENTATION OF CONTACT ELLIPSE 

Vf Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc Vc ; 


DO 888 I=1,KK-1 

WRITE (9, 889) AX1 (I) , ANG1 (I) ,AX2(I) ,ANG2(I) 
FORMAT (3X, 4 (G14 . 7 , 3X) ) 

CONTINUE 




/) 


IF ( JCL. EQ . 1) GOTO 1111 
IF ( JCL. EQ . 3) GOTO 1113 

V AND H CHECK FOR TOE POSITION 

HMT=WD+CC-3 . 0/4 . 0*FW" (DTAN (FA) +DTAN (RA) ) 
DED2T=DED2-3 . 0/4 . 0*FW*DTAN (RA) 

TMCD=MCD-0.25*FW 

XL=TMCD*DCOS (PGMA2) + (DED2T-HMT/2 . 0) *DSIN (PGMA2) 
RL=TMCD*DS IN (PGMA2) - (DED2T-HMT/2 . 0) *DCOS (PGMA2) 

FIND THE MEAN CONTACT POINT ON THE GEAR SURFACE 

ERRREL=0. ID-7 
N=2 

ITMAX=200 

IF (JCH.EQ.l) THEN 
XI (1) =270 . 0*CNST+B2 
ELSE 

XI (1) =B2 

XI (1) =90 . 0*CNST-B2 
END IF 
XI (2) =0 . 0 

CALL DNEQNF (FCN1 , ERRREL, N , UMAX , XI , X , FNORM) 
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149 


139 


TH=X (1) 
PH=X (2) 
ZY1=X(1) 
ZY2=X(2) 
N=3 


ERRREL=0. ID-10 

ITMAX=200 

XI(1)=0.0 

XI (2) =THF 

XI (3) = 0.0 

IF (JCC.EQ.l) THEN 

CALL DNEQNF (FCNM, ERRREL , N , UMAX, XI , X, FNORM) 
ELSE 

CALL DNEQNF (FCNMR, ERRREL, N, UMAX, XI, X, FNORM) 

END IF 

PHI2PO=X(l) 

XI (1) =ZY1 
XI (2) =ZY2 
XI (3) =X (2) 

XI (4) =X (3) 

XI (5) =PHI IP 
WRITE(9, 149) 

FORMAT (//6X, 

& 6X , 

& 6X , 

WRITE (9,139) 

FORMAT (//4X, 


«■- 


V AND H CHECK AT TOE POSITION 




V , H 

* 


V = ' ,G14.7, ” 


H =' ,G14. 7//) 


JCL=3 

GO TO 5555 


,/ 

/ 

,/) 


V AND H CHECK FOR HEEL POSITION 

13 HMH=WD+CCC- 1.0/4. 0*FW* (DTAN (FA) +DTAN (RA) ) 
DED2H=DED2- 1.0/4. 0*FW*DTAN (RA) 

HMCD=MCD+0. 25*FW 
1113 HMH=WD+CC-0 . 16*FW* (DTAN (FA) +DTAN (RA) ) 
DED2H=DED2-0. 16*FW*DTAN (RA) 

HMCD=MCD+0. 16*FW 

XL=HMCD*DCOS (PGMA2) + (DED2H-HMH/2 . 0) *DSIN (PGMA2) 
RL=HMCD*DSIN (PGMA2) - (DED2H-HMH/2 . 0) *DCOS (PGMA2) 
ERRREL=0. ID- 7 
N=2 

ITMAX=200 

IF (JCH.EQ.l) THEN 
XI (1) =270 . 0*CNST+B2 
ELSE 

C XI(1)=90.0*CNST-B2 
XI (1) =B2 
END IF 
XI (2) =0 . 0 

CALL DNEQNF (FCN1, ERRREL, N, ITMAX, XI, X, FNORM) 
TH=X(1) 

PH=X (2) 

ZY1=X (1) 

ZY2=X (2) 

C... FIND THE V AND H VALUE FOR HEEL POSITION 
N=3 

ERRREL=0. ID-10 
ITMAX=200 
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159 


169 


XI(1)=0.00 
XI (2) =THF 
XI (3) = 0.0 
XI (2) =THF+0. 2 
XI (3) =— 0 . 2 


IF (JCC.EQ.l) THEN 

CALL DNEQNF (FCNM , ERRREL , N , UMAX, XI ,X, FNORM) 

ELSE 

CALL DNEQNF (FCNMR , ERRREL , N , UMAX , XI , X , FNORM) 

END IF 
PHI2P0=X (1) 

XI (1) =ZY1 
XI (2) =ZY2 
XI (3) =X (2) 

XI (4) =X (3) 

XI (5) =PHI IP 

WRITE (9,11) PHI2P0 , PHI IP 
WRITE(9, 159) 

FORMAT (//6X, ******************************************* 

& 6X , ' * V AND H CHECK AT HEEL POSITION 

5, 6X, ' ****************************************** 

WRITE(9,169) V , H 

FORMAT (//4X, ' *** V = ' , G14. 7 , ' *** H =',G14.7//) 


” I 
;*C * 


/ 

/ 

/) 


JCL=1 
GOTO 5555 
1111 CONTINUE 

IF (KSIDE . EQ . 0) GOTO 1990 

STOP 

END 

. . . FCN1 IS TO FIND THE MEAN CONTACT POINT 

SUBROUTINE FCN1(X,F,N) 

IMPLICIT REAL*8 (A-H.O-Z) 

INTEGER N 

REAL*8 X(N) ,F(N) ,mcd 

COMMON/ Al/CNST , TNI , TN2 , C , FW , GAMMA , xl , r 1 , mcd 
COMMON/ A3/B2 , RGMA2 , FGMA2 , PGMA2 , D2R , D2F , ADD2 , DED2 , WD , CC , D2P 
COMMON/ AA/ SR2 , Q2 , RC2 , PW2 , XB2 , XG2 , EM2 , GaMA2 , CR2 , ALP2 , PHI2 , PHI2P 
COMMON/ A5 /SG , XM, YM, ZM, XNM, YNM , ZNM ,X2M, Y2M.Z2M, XN2M , YN2M, ZN2M , 
S.XNH2 , YNH2 , ZNH2 , XH2 , YH2 , ZH2 
TH=X (1) 

PH=X (2) 

SP=DSIN(ALP2) 

CP=DCOS (ALP2) 

SM=DSIN (GAMA2) 

CM=DCOS (GAMA2) 

STP=DSIN (TH~PH) 

CTP=DCOS(TH-PH) 

XNM=-CP*CTP 
YNM=-CP*STP 
ZNM= SP 

AA1=RC2*STP+SR2*DS IN (-Q2-PH) 

AA2=RC2*CTP+SR2*DCOS (-Q2-PH) 

AX=-EM2*SM 
AY= XB2*CM 
AZ= EM2*CM 


C 

C. . . 


FIND SG 
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C 


Tl= XNM* (AX-AAl* (SM-CR2) ) +YNM* (AY+AA2* (SM-CR2) ) +ZNM* (AZ+AAl*CM) 

T2=-XNM*(SM-CR2)*SP*STP+YNM*((SM-CR2)*SP*CTP-CP*CM)+ZNM*CM*SP*STP 

SG=T1/T2 

XM= (RC2-SG*SP) *CTP+SR2*DC0S (-Q2-PH) 

YM= (RC2~SG*SP) *STP+SR2*DSIN(-Q2-PH) 

ZM=-SG*CP 
XM=-SG*SP*CTP+AA2 
YM— SG*SP*STP+AAl 
ZM=-SG*CP 

xX= CM*XM+SM*ZM-XG2-XB2*SM 
yY= YM+EM2 

zZ=~SM*XM+CM*ZM-XB2*CM 
XN= CM*XNM+SM*ZNM 
YN= YNM 

ZN— SM*XNM+CM*ZNM 
PHI2=PH/CR2 
sh2=ds in (phi2) 
ch2=dcos (phi2) 

X2M= xX 

Y2M= CH2*yY~SH2*zZ 
Z2M= SH2*yY+CH2*zZ 
XN2M= XN 

YN2M= CH2*YN-SH2*ZN 
ZN2M= SH2*YN+CH2*ZN 
F(1)=X2M-XL 

F (2) =Yy**2+Zz**2-RL**2 
F (2) = Y 2 M • ' * 2 + Z 2 M * * 2 - R L * * 2 
RETURN 
END 

. . . SUBROUTINE CURVAl IS TO COMPUTER THE CURVATURE OF THE 
GENERATED SURFAFE 

SUBROUTINE CURVAl 
IMPLICIT REAL ; -8(A-H,0-Z) 

REAL*8 KS,KQ,KF,KH 

DIMENSION ESN (3) ,EQN(3) ,W1VT2(3) ,WV12(3) ,W2VT1 (3) 

COMMON/ A6/ES (3) , EQ (3) ,CN(3) ,Wl(3) , W2 (3) , W12 (3) ,VT1 (3) ,VT2(3) , 

$V12 (3) , KS , KQ , KF,KH , EF (3) ,EH (3) , SIGSF, PI21 

C. . . 

ESN (1) = CN (2) *ES (3) -CN (3) *ES (2) 

ESN (2) =- (CN ( 1 ) *ES (3) -CN (3) *ES (1) ) 

ESN (3) = CN (1) *ES (2) -CN (2) *ES (1) 

C. . . 

EQN (1) = CN (2) "EQ (3) -CN (3) ,V EQ (2) 

EQN (2) =- (CN (1) *EQ (3) -CN (3) *EQ (I) ) 

EQN (3) = CN (1) *EQ (2) -CN (2) *EQ (1) 

C. . . 

W1 VT2 (1) = W1 (2) *VT2 (3) -W1 (3) *VT2 (2) 

W1VT2 (2) =- (W1 (1) *VT2 (3) -W1 (3) *VT2 (1) ) 

W1VT2 (3) = W1 (1) *VT2 (2) -W1 (2) *VT2 (1) 

C. . . 

W2VT1 (1) = W2 (2) ,v VTl (3) -W2 (3) *VT1 (2) 

W2VT1 (2) =- (W2 (1) *VT 1 (3) -W2 (3) *VT1 (1) ) 

W2VT1 (3) = W2 (1) *VT1 (2) -W2 (2) *VTl (1) 

C. . . 

WV12 (1) = W12 (2) *V12 (3) -W12 (3) *V12 (2) 

WV12 (2) =- (W12 (1) *V12 (3) -W12 (3) *Vl2 (1) ) 

WV12 (3) = W12 (1) *V12 (2) -W12 (2) *V12 (1) 
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V12S=0. 0 
V12Q=0. 0 
WNES=0. 0 
WNEQ=0 . 0 
VWN= 0.0 
W1TN=0 . 0 
W2TN=0. 0 
VT2N=0.0 

C. . . 

DO 1 1=1,3 

V12S= V 1 2 ( I) *ES (I) +V12S 
V12Q- V12 (I) *EQ (I) +V12Q 
WNES= W12 (i) *ESN (i) +WNES 
WNEQ* W12 (i) ,C EQN (I) +WNEQ 
VWN = CN (I) *WV12 (i) +VWN 
W1TN= CN (i) ’ C W1VT2 (I) +W1TN 
W2TN= CN(I)*W2VT1 (I)+W2TN 
VT2N= CN (I) *VT2 (I) +VT2N 

1 CONTINUE 

. . . COMPUTER THE CURVATURE OF THE GENERATED SURFACE 

A13=-KS*V12S-WNES 
A23=-KQ*V 1 2Q- WNEQ 

A33=KS*Vl2S**2+KQ*V12Q**2-VWN-WlTN+W2TN+PI21*VT2N/W2(3) 
Tl=2.0D00*Al3*A23 
T2=A23* , ' t 2-Al3**2+ (KS-KQ) *a33 
SIG1F=0.5D00*DATAN2(T1,T2) 

KF=0 . 50D00* (KS+KQ) -0 . 5D00* (A13**2+A23**2) /A33 
&+Al3*A23/ (A33*DS IN (2 . 0D00*S IG 1 F) ) 

KH= KF-2.0D00*A13*A23/ (A33*DSIN (2 . 0D00- : SIG1F) ) 

SIGSF=SIG1F 
DO 2 1=1,3 

EF (I) = DCOS (SIG1F) *ES (I) -DSIN (SIG1F) *EQ (i) 

EH (i) = DSIN (SIG1F) *ES (I) +DCOS (SIG1F) *EQ (I) 

2 CONTINUE 
RETURN 
END 

. . . FCN2 IS TO FIND THE INITIAL GEAR ROTATIONAL ANGLE 

SUBROUTINE FCN2(X,F,N) 

IMPLICIT REAL *8 (A-H,0-Z) 

INTEGER N 
REAL *8 X (N) , F (N) 

COMMON/ Al/CNST, TNI , TN2 , C , FW , GAMMA, XL , RL , MCD 
COMMON/A5/SG , XM , YM , ZM , XNM , YNM , ZNM , X2M, Y2M.Z2M, XN2M , YN2M , ZN2M, 
&XNH2 , YNH2 , ZNH2 , XH2 , YH2 , ZH2 
CM=DCOS (GAMMA) 

SM=DS IN (GAMMA) 

CHP=DCOS(X(l)) 

SHP=DSIN(X(1)) 

XX= X2M 

YY=-Y2M*CHP+Z2M*SHP 
ZZ=-Y2M*SHP-Z2M*CHP 
XH2= XX*CM+ZZ*SM 
YH2= YY+C 
ZH2=-XX*SM+ZZ*CM 
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XX= XN2M 

YY— YN2M*CHP+ZN2M*SHP 
ZZ=-YN2M*SHP-ZN2M*CHP 
XNH2= XX*CM+ZZ*SM 
YNH2= YY 

ZNH2=-XX*SM+ZZ*CM 

R12=TNl/TN2 

V12X=- (YH2-C) *SM*R12 

V12Y= XH2*SM*R12+ (1 . 0+R12*CM) *ZH2 

V12Z—YH2* (1 . 0+Rl2*CM) +C*CM*R12 

F (1) =XNH2*V12X+YNH2"'V12Y+ZNH2*V12Z 

RETURN 

END 

. . . THE FOLLOWING IS THE TCA SUBROUTINE FOR CURVED BLADE 

SUBROUTINE FCNR(X,F,N) 

IMPLICIT REAL*8(A-H,0~Z) 
real *8 x (N) , f (N) 

DIMENSION CH (3) , P (3) , E1EF (3) ,ESN(3) ,EQN(3) ,W1VT2(3) ,WV12(3) , 
$W2VT1 (3) , EFIH (3) ,EFIIH(3) ,RH(3) ,GNH(3) ,E2IH(3) ,E2IIH(3) , 
&E1IH(3) ,E1IIH(3) ,EFI (3) , EFII (3) , Ell (3) , El 1 1 (3) ,GN(3) , EFEl (3) , 
&ERR (20) , XP (20) , YP (20) 

COMMON/Al /CNST , TNI , TN2 , C , FW , GAMMA, x 1 , r 1 , mcd 
COMMON/A2/B 1 , RGMAl , FGMAl , PGMAl , DlR , D1 F , ADDl , DED1 
COMMON/A3/B2 , RGMA2 , FGMA2 , PGMA2 , D2R.D2F, ADD 2 , DED2 , WD , CC , D2P 
C0MM0N/A4/SR2 , Q2 , RC2 , PW2 , XB2 , XG2 , EM 2 , GaMA2 , CR2 , ALP2 , PHI2 , PHI2P 
COMMON /A5/SG , XM, YM, ZM,XNM, YNM, ZNM , X2M, Y2M, Z2M ,XN2M, YN2M, ZN2M, 
&XNH2 , YNH2 , ZNH2 , XH2 , YH2 , ZH2 

COMMON/A6/ES (3) ,EQ(3) ,CN(3) ,W1 (3) ,W2(3) ,W12(3) ,VT1 (3) , VT2(3) , 

$ V 12(3) , KS , KQ , KF , KH , EF (3) ,EH(3) ,SIGSF,PI21 
COMMON/ A7/SR 1 , Q 1 , Rc £ , PW1 , XB 1 , XG1 , EMI , GaMAl , CR 1 , ALP 1 , PHI 1 , PHI IP 
COMMON/A8/Sf ,XM1 , YM1 , ZM1 ,XNM1 , YNM1 , ZNM1 , XlM, Y 1M, Z 1M , 

&XN1M, YN1M, ZN1M, XNH1 , YNH1 , ZNH1 ,XH1 , YH1 , ZH1 
COMMON/A9/PHI2PO , OX , OZ , XO , ZO , RHO , ALP , V , H , CR1T , PCR1T 
COMMON/Al 1/RAM, PSI 1,C2,D6,E24, F120, CX6.DX2 A, EX120, RU1 ,DELT,RUP, 
$RAl , CPF, DPF , EPF , FPF 
TH=X (1) 

PH=X (2) 

SP=DSIN(ALP2) 

CP=DCOS(ALP2) 

SM = DS IN (GAMA2) 

CM=DCOS (GAMA2) 

STP=DSIN (TH-PH) 

CTP=DCOS (TH-PH) 

XNM=-CP*CTP 
YNM=-CP*STP 
ZNM= SP 

AAl=RC2*STP+SR2*DSIN (-Q2-PH) 

AA2=RC2*CTP+SR2*DCOS (-Q2-PH) 

AX=-EM2*SM 
AY= XB2*CM 
AZ= EM2*CM 


. .. FIND SG 


Tl= XNM* (AX-AAl* (SM-CR2) ) +YNM* (AY+AA2* (SM-CR2) ) +ZNM* (AZ+AAl *CM) 
T2=-XNM* (SM-CR2) *SP*STP+YNM* ( (SM-CR2) *SP*CTP-CP*CM) +ZNM*CM*SP*STP 
SG=Tl/T2 
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XM= (RC2-SG*SP) *CTP+SR2*DC0S (-Q2-PH) 

YM= (RC2-SG*SP) *STP+SR2*DSIN(-Q2-PH) 

ZM=~SG*CP 

C XM=~SG*SP*CTP+AA2 

C YM=~SG*SP*STP+AAl 

C ZM=~SG*CP 

xX= CM*XM+SM*ZM-XG2-XB2*SM 
yY= YM+EM2 

zZ— SM*XM+CM*ZM-XB2*CM 
XN= CM*XNM+SM*ZNM 
YN= YNM 

ZN— SM*XNM+CM*ZNM 
PHI2=PH/CR2 
sh2=ds in (phi2) 
ch2=dcos (phi2) 

X2M= xX 

Y2M= CH2*yY-SH2*zZ 
Z2M= SH2*yY+CH2*zZ 
XN2M= XN 

YN2M= CH2*YN~SH2*ZN 
ZN2M= SH2*YN+CH2*ZN 
CMM=DCOS (GAMMA) 

SMM=DS IN (GAMMA) 

CHP=DC0S(PHI2P) 

SHP=DSIN(PHI2P) 

XX= X2M 

YY— Y2M*CHP+Z2M*SHP 
ZZ=-Y2M*SHP-Z2M*CHP 
XH2= XX*CMM+ZZ*SMM 
YH2= YY+C+V 
ZH2=~XX*SMM+ZZ*CMM 

C. . . 

XX= XN2M 

YY— YN2M*CHP+ZN2M*SHP 
ZZ— YN2M*SHP~ZN2M*CHP 
XNH2= XX*CMM+ZZ*SMM 
YNH2= YY 

ZNH2=-XX*SMM+ZZ*CMM 

.... DEFINE THE PINION SURFACE 

TH1=X (3) 

PH1=X (4) 

SM1=DSIN (GAMAl) 

CM1=DC0S (GAMAl) 

STP=DSIN (TH1+PH1) 

CTP=DCOS (TH1+PH1) 

... FIND CR1T,PF,PPF,PCR1T 

DDD=DABS (PHI) 

IF (DDD.LE. 0.001) GOTO 6 

PHIl=RAl*(PHl-CPF*PHl**2-DPF*PHl**3~EPF*PHl**4-FPF*PHl**5) 
PF-RA1 * (1 . 0-2 . 0*CPF*PHl-3 . 0*DPF*PH1**2 
$-4 . 0*EPF*PH1**3~5 . 0*FPF*PH1**4) 

PPF=-RAl* (2 . 0*CPF+6 . 0*DPF*PH1+12 . 0*EP F*PH 1 **2+20. 0*FPF*PH 1**3) 
CR1T=1 . 0/PF 
PCR1T=-PPF/PF**3 
GOTO 7 

6 PHI 1=RA1*PH1 
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CR1T=CR1 

PCR1T-2.0*CPF/(RA1**2) 

7 CONTINUE 
CR1T=CR1 
PCR1T=0. 000 

FIND THE NOMAL OF THE EQUIDISTANCE SURFACE 

XMO= XO^'CTP+SRl^DCOS (-Q1+PH1) 

YMO= X0*STP+SR1*DSIN(-Q1+PH1) 

ZMO= ZO 

V1X=-YM0*SM1-EM1*SM1 
V1Y= XM0*SM1-(ZM0-XB1)*CM1 
V1Z= YMO *CM 1 +EM 1 *CM 1 
V2X— YM0*CR1T 
V2Y= XM0*CR1T 
V2Z= 0.0 
VX=V1X-V2X 
VY=V1Y-V2Y 
VZ=V1Z-V2Z 
TX=-CTP 
TY=-STP 
TZ=0. 0 
FX= STP 
FY=-CTP 
FZ=0. 0 

XNN= FY*VZ~FZ*VX 
YNN= FZ*VX~FX*VZ 
ZNN= FX*VY-FY*VX 
DDD=DSQRT(XNN**2+YNN**2+ZNN**2) 
XNM1=XNN/DDD 
YNM1=YNN/DDD 
ZNMl=ZNN/DDD 

DT=TX*XNM1+TY*YNM1+TZ*ZNM1 
IF(DT.GE.O.O) GOTO 10 
XNM1=-XNM1 
YNM1=-YNM1 
ZNM1»-ZNM1 
10 CONTINUE 

XM1- XMO-RHO*XNMl 
YM1= YMO-RHO*YNMl 
ZM1- ZMO-RHO*ZNMl 

ALP=DARCOS (TX*XNM1+TY*YNM1+TZ*ZNM1) 
xX= CM1*XM1+SM1*ZM1-XG1-XB1*SM1 
yY= YM1+EM1 

zZ=-SMl*XMl+CMl*ZMl-XBl*CMl 

XN1=CM1*XNM1+SM1*ZNM1 

YN1=YNM1 

ZN1— SM1*XNM1+CM1*ZNM1 
C PHI l=PHl/CRl 

shl=ds in (phi 1) 
chl=dcos (phi 1) 

X1M- xX 

Y 1M= CHl*yY+SHl*zZ 
ZlM=-SHl*yY+CHl*zZ 
XN1M= XN1 

YN1M= CH 1 * YN 1 +SH 1 *ZN 1 
ZN1M—SH1*YN1+CH1*ZN1 
PHI1P=X(5) 
shlP=ds in (phi IP) 
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chlP=dcos (philP) 

XH1= X1M+H 

YH1- CH1P*Y1M~SH1P*Z1M 
ZH1= SH1P*Y1M+CH1P*Z1M 
XNH1= XN1M 

YNH1= CH1P*YN1M~SH1P*ZN1M 

ZNH1= SH1P*YN1M+CH1P*ZN1M 

F(1)=XH2-XH1 

F(2)=YH2-YH1 

F(3)=ZH2-ZH1 

F(4)=XNH2-XNH1 

F(5)=ZNH2-ZNH1 

RETURN 

END 

... THE FOLLOWING IS THE SUBROUTINE FOR STRAIGHT BLADE 

SUBROUTINE FCN(X,F,N) 

IMPLICIT REAL*8(A-H,0-Z) 
rear-8 x (N) , f (N) 

DIMENSION CH(3) ,P (3) ,E1EF(3) ,ESN(3) ,EQN(3) ,W1VT2(3) ,WV12(3) , 

$W2VTl (3) ,EFIH(3) ,EFIIH(3) ,RH(3) ,GNH(3) ,E2IH(3) ,E2IIH(3) , 

&E1IHC3) , E1IIH (3) , EFI (3) ,EFII (3) , Ell (3) , El II (3) ,GN(3) ,EFEl(3) , 

&ERR (20) , XP (20) ,YP(20) 

COMMON/Al/CNST , TNI , TN2 , C , FW , GAMMA, xl , r 1 , mcd 

C0MM0N/A2/B 1 , RGMAl , FGMAl , PGMAl , DIR , DIF , ADD1 , DEDl 

COMMON/ A3 / B2 , RGMA2 , FGMA2 , PGMA2 , D2R , D2 F , ADD 2 , DED2 , WD , CC , D2P 

COMMON/A4/SR2 ,Q2,RC2, PW2 , XB2 ,XG2 , EM2 , GaMA2 , CR2 , ALP 2 , PHI2 , PHI2P 

COMMON/ A5 / SG , XM , YM , ZM , XNM , YNM, ZNM, X2M , Y2M , Z2M , XN2M , YN2M , ZN2M, 

&XNH2 , YNH2 , ZNH2 , XH2 , YH2 , ZH2 

C0MM0N/A6/ES (3) ,EQ(3) ,CN(3) ,Wl(3) ,W2(3) ,W12(3) ,VT1(3) ,VT2(3) , 

$ V 12(3) , KS , KQ , KF , KH , EF (3) ,EH(3) ,SIGSF,PI21 
COMMON/ A7/SR1 ,Q1 , Ref , PWl ,XB1 ,XG1 , EMI ,GaMAl , CR1 , ALP1 , PHI 1 , PHI IP 
C0MM0N/A8/Sf ,XM1 , YM1 , ZM1 , XNM1 , YNM1 , ZNMl ,X1M, Y1M.Z1M, 

&XN1M, YN1M, ZN1M, XNH1 , YNHI , ZNH1 , XH1 , YH1 , ZH1 
COMMON/A9/PHI2PO , OX , OZ , XO , ZO , RHO , ALP , V , H , CR1T , PCR1T 
COMMON/Al 1/RAM, PSI 1 , C2 , D6 , E24 , FI 20 , CX6 , DX24 , EX120 , RU1 , DELT , RUP , 
$RAl , CPF , DPF , EPF , FPF 
TH=X (1) 

PH=X (2) 

SP=DSIN(ALP2) 

CP=DCOS (ALP2) 

SM=DSIN(GAMA2) 

CM=DCOS (GAMA2) 

STP=DSIN(TH-PH) 

CTP=DCOS (TH-PH) 

XNM=~CP*CTP 
YNM=-CP*STP 
ZNM= SP 

AA1-RC2*STP+SR2*DSIN (-Q2-PH) 

AA2=RC2*CTP+SR2*DCOS (-Q2-PH) 

AX=~EM2*SM 
AY= XB2*CM 
AZ= EM2*CM 
C 

C. . . FIND SG 
C 

Tl= XNM* (AX-AA1* (SM-CR2) ) +YNM* (AY+AA2* (SM-CR2) ) +ZNM* (AZ+AAl*CM) 
T2=-XNM* (SM-CR2) *SP*STP+YNM* ( (SM-CR2) *SP*CTP-CP*CM) +ZNM*CM*SP*STP 
SG=T1/T2 
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XM= (RC2-SG*SP) *CTP+SR2*DC0S (-Q2-PH) 

YM= (RC2-SG*SP) *STP+SR2*DSIN (-Q2-PH) 

ZM=-SG*CP 

XM=-SG*SP*CTP+AA2 

YM=~SG*SP*STP+AAl 

ZM=~SG*CP 

xX= CM*XM+SM*ZM-XG2~XB2*SM 
yY= YM+EM2 

zZ=-SM*XM+CM*ZM~XB2*CM 
XN= CM*XNM+SM*ZNM 
YN= YNM 

ZN=-SM*XNM+CM*ZNM 
PHI2=PH/CR2 
sh2=ds in (phi2) 
ch2=dcos (phi2) 

X2M= xX 

Y2M= CH2*yY~SH2*zZ 
Z2M= SH2*yY+CH2*zZ 
XN2M= XN 

YN2M= CH2*YN-SH2*ZN 
ZN2M= SH2*YN+CH2*ZN 
CMM=DC0S (GAMMA) 

SMM=DS IN (GAMMA) 

CHP=DCOS (PHI2P) 

SHP=DSIN (PHI2P) 

XX= X2M 

YY=-Y2M*CHP+Z2M*SHP 
ZZ=-Y2M*SHP-Z2M*CHP 
XH2= XX*CMM+ZZ*SMM 
YH2= YY+C+V 
YH2= YY+C-V 
ZH2=-XX*SMM+ZZ*CMM 

XX= XN2M 

YY=-YN2M*CHP+ZN2M*SHP 
ZZ=-YN2M*SHP-ZN2M*CHP 
XNH2= XX* CMM+ Z Z * SMM 
YNH2= YY 

ZNH2=-XX*SMM+ZZ*CMM 


. . DEFINE THE PINION SURFACE 

TH1=X (3) 

PH1=X (4) 

SP1=DSIN (-ALP1) 

CP1=DC0S (-ALP1) 

SM1=DSIN (GAMAl) 

CM1=DC0S (GAMAl) 

STP=DSIN (TH1+PH1) 

CTP=DCOS (TH1+PH1) 
XNM1=-CP1*CTP 
YNM1=-CP1*STP 
ZNM1= SP1 

AB1=RCF*STP+SR1*DSIN (-Q1 + PH1) 

AB2=RCF*CTP+SR1*DC0S (-Q1+PH1) 

AXX=-EM1*SM1 

AYY= XB1*CM1 

AZZ= EM1*CM1 


. . FIND SF,CR1T,PF,PPF,PCR1T 
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c 


DDD=DABS (PHI) 

IF (DDD.LE. 0.001) GOTO 6 

PHIl=RAlMPHl-CPF*PHl**2-DPF*PHl , '” v 3-EPF*PHr'”‘4-FPF*PHl**5) 
PF-RAl* (1 . 0-2 . 0*CPF*PHl-3 . 0*DPF*PH1**2 
$ - 4 . 0 * E P F * P H 1 * * 3 - 5 . 0 * F P F * P H 1 * * 4) 

PPF— RA1* (2. 0*CPF+6 . 0*DPF*PH1+12. 0*EPF*PHl**2+20. 0*FPF*PH1**3) 

CR1T-1.0/PF 

PCR1T— PPF/PF**3 

GOTO 7 

6 PHI1=RA1*PH1 
CR1T=CR1 

PCRlT=2 . 0*CPF/ (RAl**2) 

7 CONTINUE 

C CR1T=CR1 

C PCR1T=0. 000 

Tl= XNM1* (AXX-AB1* (SM1-CR1T) )+ 

&YNM1* (AYY+AB2* (SM1-CR1T) ) +ZNM1* (AZZ+ABl*CMl) 

T2— XNM1* (SM1-CR1T) *SP1*STP+ 

&YNM1* ( (SM1-CR1T) *SP1*CTP _ CP1*CM1) +ZNM1*CM1*SP1*STP 
SF=T1/T2 
C 

XM1- (RCF~SF*SP1) *CTP+SR1*DC0S (-Q1+PH1) 

YM1= (RCF~SF*SP1) *STP+SR1*DSIN (-Q1+PH1) 

ZM1— SF*CP1 

xX= CM 1 *XM 1 +SM1 *ZM 1 -XG 1 -XB 1 *SM 1 
yY= YMl+EMl 

zZ=- SM 1 *XM 1+CM1 *ZM 1 -XB 1 * CM 1 
XN 1 =CM 1 *XNM 1+SM1*ZNM1 
YN1=YNM1 

ZN1=-SM1*XNM1+CM1*ZNMI 
C PHI 1=PH1 /CR1 

shl=ds in (phi 1) 
chl=dcos (phi 1) 

X1M= xX 

Y1M= CHl*yY+SHl*zZ 
ZlM=-SHl*yY+CHl*zZ 
XN1M= XN1 

YN1M= CH1*YN1+SH1*ZN1 
ZN1M=-SH1*YN1+CH1*ZN1 
C WRITE (9,111) XlM, YlM, ZlM 

PHI1P= X (5) 
shlP=ds in (phi IP) 
chlP=dcos (philP) 

XH1= X1M+H 

YH1= CH1P*Y1M~SH1P“Z1M 
ZH1= SH1P*Y1M+CH1P*Z1M 
XNH1= XN1M 

YNH1= CH1P*YN1M-SH1P*ZN1M 

ZNH1= SH1P*YN1M+CH1P*ZN1M 

F(1)=XH2-XH1 

F(2)=YH2-YH1 

F(3)=ZH2-ZH1 

F(4)=XNH2-XNH1 

F(5)=YNH2-YNH1 

RETURN 

END 


* SUBROUTINE ELLIP IS TO DETERMINE THE SIZE AND ORINTATION * 
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. .. * OF THE CONTACT ELLIPSE * 

* >V Vc Vc Vc Vc * Vc * Vc Vc Vc * * * Vc Vc * Vc * * Vc * Vc * Vc Vc V? Vc iV Vc Vc Vc Vc Vc Vc * * Vc * Vc Vc Vc * * Vc Vc Vc Vc * * * Vc * * Vc * * V? Vc 


SUBROUTINE ELLIP 
IMPLICIT REAL*8(A-H,0-Z) 

REAL* 8 KS , KQ ,K1I , K1 II ,K2I , K2II 

DIMENSION RO (3) ,ETA2(3) ,ZETA2(3) ,E1E2 (3) , ETA(3) ,ZETA(3) 
DIMENSION El IH (3) ,ElIIH(3) ,E2IH(3) ,E2IIH(3) ,GNH(3) 

COMMON/ Al/CNST , TNI , TN2 , C , FW , GAMMA, x 1 , r 1 , mcd 
COMMON/A3/B2 , RGMA2 , FGMA2 , PGMA2 , D2R , D2F , ADD2 , DED2 , WD , CC , D2P 
COMMON/ A4/ SR2 , Q2 , RC2 , PW2 , XB2 , XG2 , EM2 , GaMA2 , CR2 , ALP2 , PHI2 , PHI2P 
COMMON/A5 / SG , XM , YM , ZM , XNM , YNM , ZNM , X2M , Y 2M , Z 2M , XN2M , YN2M , ZN2M , 
&XNH2 , YNH2 , ZNH2 , XH2 , YH2 , ZH2 
COMMON/ A9/PHI2P0, OX , OZ , XO , ZO , RHO , ALP , V , H , CR IT , PCR IT 
COMMON/AlO/KlI ,K1II,K2I,K2II, DEL , ElIH , E1IIH , E2IH, E2IIH, GNH, 
&A2P , B2P , TAU1R, TAU2R , A2L , B2L 
CNST=DARCOS (-1 . ODOO) / 180 . 00 

E1E2 (1) = El IH (2) *E2IH(3)-E1IH(3)*E2IH(2) 

E1E2 (2) =- (El IH (1) *E2lH (3)-ElIH(3) *E2IH (1) ) 

E1E2 (3) = E1IH(1) *E2IH(2) -ElIH (2) *E2IH(1) 


T1=0 . 0 
T2=0 . 0 
DO 1 1=1,3 

Tl= ElIH (I) *E2IH (I) +T1 
T2= GNH (i) *E1E2 (i) +T2 
1 CONTINUE 

Tl=Tl+l . ODOO 

SIG12=2 . 0D00*DATAN2 (T2 , Tl) 

SK1= K1I+K1II 
SK2= K2I+K2II 
SG1= K1I-K1II 
SG2= K2I-K2II 

T1=SG1-SG2*DC0S (2 . 0D00*SIG12) 

T2=SG2*DSIN(2.0D00*SIG12) 

T3=DSQRT (SGl**2+SG2**2-2 . ODOO*SG1*SG2*DCOS (2 . ODOO*SIG12) ) 

TX=T2/T3 
TY=T1/T3+1 .ODOO 
ALP12=DATAN2 (TX,TY) 

. . . THE DIRECTION AND LENGTH OF THE AXES OF CONTACT ELLIPSE 

DEL=0 . OO7OODO0 

AL=0. 25D00* (SK1-SK2-T3) 

BL=0. 25D00* (SK1-SK2+T3) 

WRITE (9 , 5) SIG12,AL,BL 
5 FORMAT(2X, 'SIG12, AL,BL= ' , 3 (2X , G14 . 7) ) 

AL=DABS (AL) 

BL=DABS(BL) 

A2L=2 . ODOO*DSQRT (del/al) 

B2L=2 . 0D00*DSQRT (DEL/BL) 


DO 2 1=1,3 

ETA(I) = DCOS (ALP 12) *E1 IH (l) -DSIN (ALP1 2) *E1IIH (i) 
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ZETA (I) = DSIN (ALP 12) *E1 IH (i) +DCOS (ALP 12) *E1 IIH (I) 

2 CONTINUE 

. .. DETERMINE THE PROJECTION OF CONTACT ELLIPS IN AXIAL SECTION 

CHP=DCOS (PHI2P) 

SHP=DSIN (PHI2P) 

CMM=DCOS (GAMMA) 

SMM=DS IN (GAMMA) 

XX= ETA(1)*CMM-ETA(3)*SMM 
YY= ETA (2) 

ZZ= ETA(l) *SMM+ETA(3) 5t CMM 
ETA2 (1) = XX 
ETA2 (2) =-YY*CHP-ZZ*SHP 
ETA2 (3) = YY*SHP~ZZ*CHP 

XX= ZETA(1)*CMM~ZETA(3)*SMM 
YY= ZETA(2) 

ZZ= ZETA(1)*SMM+ZETA(3)*CMM 
ZETA2 (1) = XX 
ZETA2 (2) =-YY*CHP-ZZ*SHP 
ZETA2 (3) = YY*SHP-ZZ*CHP 

RO (2) =Y2M/DSQRT (Z2M**2+Y2M**2) 

RO (3) =Z2M/DSQRT (Z2M**2+Y2M**2) 

R0(1)=0.0D00 

T1 1=0 . 0D00 
T12=0. 0D00 
DO 3 1=1,3 

T12= ETA2 (I) *R0 (I) +T12 
Tll=ZETA2(l)*R0(l)+Tll 

3 CONTINUE 

TAU 1 = DAT AN 2 (Til, ZETA2 (1) ) 

TAU2=DATAN2 (T 1 2 , ETA2 ( 1 ) ) 

A2P=A2L*ZETA2 (1) /DCOS (TAU1) 

B2P=B2L*ETA2 (1) /DCOS (TAU2) 

TAU1R= (TAU1-RGMA2) /CNST 
TAU2R= (TAU2-RGMA2) /CNST 
RETURN 
END 
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. . THE FOLLOWING IS THE V-H CHECK PROGRAM FOR CURVED BLADE 

SUBROUTINE FCNMR (X , F , N) 

IMPLICIT REAL*8(A-H,0-Z) 
real *8 x (N) , f (N) 

COMMON/ Al /CNST , TN 1 , TN2 , C , FW , GAMMA , x 1 , r 1 , mcd 
COMMON/ A5/ SG , XM , YM , ZM , XNM , YNM , ZNM , X2M , Y2M , Z2M , XN2M , YN2M, ZN2M , 
&XNH2 , YNH2 , ZNH2 , XH2 , YH2 , ZH2 

COMMON/ A7/SR1 , Q1 ,Rcf , PW1 , XBl , XG1 , EMI , GaMAl , CR1 , ALP1 , PHI 1 , PHI IP 
COMMON/ A9/PHI 2P0 , OX , OZ , XO , ZO , RHO , ALP , V , H , CR IT , PCR IT 
COMMON/A1 1/RAM, PSI 1 , C2 , D6 , E24 , F120 , CX6 , DX24 , EX1 20 , RU1 , DELT , RUP , 
$RAl , CPF , DPF , EPF , FPF 
CM=DCOS (GAMMA) 
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SM=DS IN (GAMMA) 

C CHP=DCOS (PHI2P0) 

C SHP=DSIN (PHI2P0) 

CHP=DC0S(X(1)) 

SHP=DSIN(X(1)) 

XX= X2M 

YY=-Y2M*CHP+Z2M*SHP 
ZZ=-Y2M*SHP-Z2M*CHP 
XH2= XX*CM+ZZ*SM 
YH2= YY+C 
ZH2=-XX*SM+ZZ*CM 

XX= XN2M 

YY— YN2M*CHP+ZN2M*SHP 
ZZ=-YN2M*SHP-ZN2M*CHP 
XNH2= XX*CM+ZZ*SM 
YNH2= YY 

ZNH2=-XX*SM+ZZ*CM 

... DEFINE THE PINION SURFACE 

TH1=X (2) 

PH1=X (3) 

SM1=DSIN (GAMAl) 

CM1=DC0S (GAMAl) 

STP=DSIN (TH1+PH1) 

CTP=DCOS (TH1+PH1) 

... FIND CR1T , PF , PPF , PCR1T 

DDD=DABS (PHI) 

IF (DDD.LE. 0.001) GOTO 6 

PHI l=RAl !< (PH1-CPF*PH1**2-DPF*PH1**3“EPF*PH1**A-FPF*PH1**5) 
PF=RAl* (1 . 0-2 . 0*CPF*PH1~3. 0*DPF*PH1 **2 
$-4 . 0*EPF*PHl**3-5 . 0*FPF*PHl**4) 

PPF—RAl* (2 . 0*CPF+6 . 0*DPF*PH1 + 12. 0*EPF*PHl* 5 '2+20 . 0*FPF*PH1**3) 

CR1T-1.0/PF 

PCR1T— PPF/PF**3 

GOTO 7 

6 PHI1-RA1*PH1 
CR IT=CR 1 

PCR1T=2 . 0*CPF/ (RA1**2) 

7 CONTINUE 

FIND THE NOMAL OF THE EQUIDISTANCE SURFACE 

XMO= XO*CTP+SRl*DCOS (-Q1+PH1) 

YMO= XO*STP+SRl*DSIN (-Q1+PH1) 

ZMO= ZO 

V1X— YM0*SM1-EM1*SM1 

V1Y= XM0*SM1- (ZMO-XB1) *CM1 

VIZ- YM0*CM1+EM1*CM1 

V2X=-YM0*CR1T 

V2Y= XMO*CR IT 

V2Z= 0.0 

VX=V1X-V2X 

VY=V1 Y-V2Y 

VZ=V 1Z-V2Z 

TX=-CTP 

TY=-STP 
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TZ=0. 0 
FX= STP 
FY=-CTP 
FZ=0 . 0 

XNN= FY*VZ-FZ*VX 

YNN= FZ*VX-FX*VZ 

ZNN= FX*VY~FY*VX 

DDD=DSQRT (XNN**2+YNN 5 '*2+ZNN**2) 

XNMl=XNN/DDD 

YNM1=YNN/DDD 

ZNM1=ZNN/DDD 

DT=TX*XNM1+TY*YNM1+TZ*ZNM1 
IF(DT.GE.O.O) GOTO 10 
XNMl=-XNMl 
YNMl=-YNMl 
ZNM1=-ZNM1 
10 CONTINUE 

XMl = XM0-RH0*XNM1 
YM1- YMO-RHO*YNMl 
ZM1= ZM0-RH0*ZNM1 

ALP=DARC0S(TX*XNM1+TY*YNM1+TZ*ZNM1) 
xX= CM 1 *XM 1 + SM 1 *ZM 1 -XG 1 -XB 1 *SM 1 
yY= YM1+EM1 

ZZ=-SM1*XM1+CM1*ZM1-XB1*CM1 

XN1=CM1*XNM1+SM1*ZNM1 

YN1=YNM1 

ZN 1 =- SM 1 * C XNM 1 +CM 1 *ZNM 1 
C PHI1-PH1/CR1 

shl=ds in (phi 1) 
chl=dcos (phi 1) 

X1M= xX 

Y1M= CHl*yY+SHl*zZ 
ZlM=-SHl*yY+CHl*zZ 
XN1M= XN1 

YN1M= CH 1 * YN 1 +SH 1 *ZN 1 
ZN1M=-SH1*YN1+CH1*ZN1 


TT=YN1M**2+ZN1M**2 

SH1P= (-ZN1M*YNH2+YN1M*ZNH2) /TT 

CH1P= ( YN 1M* YNH2+ZN 1M*ZNH2) /TT 

PHI 1P=2 . 0D00*DATAN2 (SHIP , (1 . 0D00+CH1P) ) 


XH1= XlM 

YH1- CH1P*Y1M-SH1P*Z1M 
ZH1- SH1P*Y1M+CH1P*Z1M 
XNH1= XN1M 

YNHl= CH1P*YN1M-SH1P*ZN1M 
ZNH1= SH1P*YN1M+CH1P*ZN1M 
V=- (YH2-YH1) 

H=XH2-XH1 

F(1)=ZH2-ZH1 

F(2)=XNH2-XNH1 

F(2)=YNH2**2+ZNH2**2-TT 


R12=TN1/TN2 
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V12X=0.0-YH2*SM*R12 

V12Y= ZH1+R12* (XH2 , 'SM+ZH2 , ' t CM) 

V12Z=-YH1-R12*YH2*CM 
C V12X=- (YH2- (C~V) ) *SM*R12 

C V12Y= XH2*SM*R12+(1.0+R12*CM)*ZH2 

C V12Z—YH2* (1 . 0+R12*CM) + (C-V) *CM*R12 

F (3) =XNH2*V 12X+ YNH2*V 1 2 Y+ZNH2*V 1 2Z 
RETURN 
END 

... THE FOLLOWING IS THE V-H CHECK SUBROUTINE FOR STRAIGHT BLADE 

SUBROUTINE FCNM(X,F,N) 

IMPLICIT REAL*8(A-H,0-Z) 
real *8 x (N) , £ (N) 

COMMON/Al /CNST , TNI , TN2 , C , FW , GAMMA, xl , r 1 , me d 

COMMON/A5/SG , XM, YM, ZM,XNM, YNM , ZNM , X2M, Y2M, Z2M, XN2M , YN2M , ZN2M , 
&XNH2 , YNH2 , ZNH2 , XH2 , YH2 , ZH2 

COMMON/ A7/SR1 , Q1 , Ref , PW1 , XB1 ,XG1 , EMI , GaMAl , CR1 , ALP1 , PHI1 , PHI1P 
COMMON/ A9/PHI2P0 , OX, OZ , XO, ZO, RHO , ALP ,V,H,CRlT,PCRlT 
COMMON/Al 1 /RAM , PS 1 1 , C2 , D6 , E24 , F120 , CX6 , DX24 , EX120 , RU1 , DELT , RUP , 
$RAl , CPF, DPF , EPF , FPF 
CM=DCOS (GAMMA) 

SM=DS IN (GAMMA) 

CHP=DCOS (X(l)) 

SHP=DSIN (X (1) ) 

XX= X2M 

YY— Y2M*CHP+Z2M*SHP 
ZZ— Y2M*SHP-Z2M*CHP 
XH2= XX*CM+ZZ*SM 
YH2= YY+C 
ZH2— XX*SM+ZZ*CM 

C. . . 

XX= XN2M 

YY— YN2M*CHP+ZN2M*SHP 
ZZ— YN2M*SHP-ZN2M*CHP 
XNH2= XX*CM+ZZ*SM 
YNH2= YY 

ZNH2— XX*SM+ZZ*CM 

C. . . 

c 

C. . . DEFINE THE PINION SURFACE 
C 

THl=X (2) 

PH1=X (3) 

SPl=DSIN (-ALP1) 

CPl=DCOS (-ALP1) 

SM1=DSIN (GAMAl) 

CMl=DCOS (GAMAl) 

STP=DSIN (TH1+PH1) 

CTP=DCOS (THl+PHl) 

XNM1— CP1*CTP 
YNM1— CP1*STP 
ZNMl= SP1 

ABl=RCF*STP+SRl *DSIN (-Q1+PH1) 

AB2=RCF*CTP+SR1*DC0S(-Q1+PH1) 

AXX— EM1*SM1 
AYY= XB1*CM1 
AZZ= EM1*CM1 
C 


174 



non noon o no 


C... FIND SF,CR1T,PF,PPF, PCR1T 
C 

PHI1=RA1*(PH1-CPF*PH1**2-DPF*PH1**3-EPF*PH1**4-FPF*PH1**5) 
PF=RAl* (1 . 0-2. 0*CPF*PHl-3 . 0*DPF*PH1**2 
$-4.0*EPF*PHl**3-5.0*FPF*PHl**4) 

PPF=-RAl* (2 . 0*CPF+6 . 0*DPF*PH1 + 12 . 0*EPF*PHl**2+20 . 0 , ’ t FPF*PHP' c *3) 
CR1T= 1 . 0/PF 
PCR1T=-PPF/PF**3 
C 

Tl= XNM1*(AXX-AB1*(SM1-CR1T))+ 

&YNM1* (AYY+AB2* (SM1-CR1T) ) +ZNM1* (AZZ+AB1*CM1) 
T2=-XNMl*(SMl-CRlT)*SPl sV STP+ 

&YNM1* ( (SM1-CR1T) *SP1*CTP~CP1*CM1) +ZNM1*CM1*SP1*STP 
SF=T1/T2 

XM1= (RCF-SF*SP1)*CTP+SR1*DC0S(~Q1+PH1) 

YH1= (RCF~SF*SPl) *STP+SRP V DSIN (-Q1+PH1) 

ZM1=-SF*CP1 

xX= CM1*XM1+SM1*ZM1-XG1-XB1*SM1 
yY= YM1+EM1 

zZ=-SH1*XM1+CM1*ZM1-XB1*CM1 
XN 1 = CM 1 *XNM 1+SM1*ZNM1 
YN1=YNM1 

ZN1=-SM1*XNM1+CM1*ZNM1 
C PHI l=PHl/CRl 

shl=ds in (phi 1) 
chl=dcos (phi 1) 

X1M= xX 

Y1M= CHP'yY+SHP'zZ 
ZlM=~SHl*yY+CHl*zZ 
XN1M= XN1 

YN1M= CH1*YNH-SH1*ZN1 
ZN1M— SH1*YN1+CH1*ZN1 
TT=YN1M**2+ZN1M**2 
SH1P= (~ZN1M*YNH2+YN1M*ZNH2) /TT 
CH1P= ( YN1M*YNH2+ZN1K*ZNH2) /TT 
PHI1P-2. 0D00*DATAN2 (SHIP , (1 . 0D00+CH1P) ) 

SH1P=DSIN (PHI IP) 

CHlP=DCOS (PHIlP) 

XH1= X1M 

YH1= CH1P*Y1M-SH1P*Z1M 
ZH1= SH1P !C Y1M+CH1P"Z1M 
XNH1= XN1M 

YNH1= CH1P*YN1M-SH1P*ZN1M 
ZNH1= SH1P*YN1M+CH1P*ZN1M 
V=-(YH2-YH1) 

H=XH2-XH1 

F(1)=ZH2-ZH1 
F(2)=XNH2-XNH1 
F (2) =YNH2**2+ZNH2**2-TT 


R12=TN1/TN2 

V12X=0.0-YH2*SM*R12 

V12Y= ZH1+R12*(XH2*SM+ZH2*CM) 

V12Z=-YH1-R12*YH2*CM 

V12X=- (YH2- (C-V) ) *SM*R12 

V 1 2 Y= XH2*SM*R 1 2+ ( 1 . 0+R 1 2*CM) *ZH2 

V12Z=-YH2* (1 . 0+Rl2*CM) + (C~V) *CM*R12 
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F(3)=XNH2*V12X+YNH2*V12Y+ZNH2*Vl2Z 

RETURN 

END 

C 

C SUBROUTINE CAM IS FOR THE COEFFICIENTS OF GENERATION MOTION 
C 

SUBROUTINE CAM AA 00070 

IMPLICIT REAL*8(A-H,0-Z) AA 00080 

COMMON/ Al 1 /RAM, PSI 1 , C2 , D6 , E24 , F120 , CX6 , DX24 , EX120 , RU1 , DELT.RUP , 

$RAl ,CPF , DPF , EPF , FPF 
Tl=l . 0+3 . 0*C2*DTAN (PSIl) 

&+ (1 . 0-RAM) **3* (RU 1**3/ 15 ,0**2+DELT) / (15 . 0*DCOS (PS II)) 

T2=l. 0+ (RU1+DELT) / (15 . 0*DCOS (PSI 1) ) 

CX6=T1/T2 

Tl= 6 . 0*C2*DCOS (PSIl) + (A. 0*CX6+3 . 0*C2**2-1 . 0) *DSIN (PSIl) 

&+6 . 0*C2* (1 . 0-RAM) * * 2 * (RU 1 * * 3 / 1 5 . 0**3+DELT/15 . 0) 

T2= DCOS (PSI 1) + (RU1+DELT) / 15 • 0 
DX24= T1/T2 

T 1 = ( 1 0 . 0*CX6+ 15 . 0*C2**2- 1 . 0) *DCOS (PSIl) 

S, + (5 . 0*DX24+10 . 0*C2*CX6-10.0*C2) *DSIN (PSIl) 

& + (10. 0*CX6*(1. 0-RAM) **2 

& +15 . 0*C2**2* (1 . 0-RAM) ) * (RU 1**3/ 15 . 0**3+DELT/15 . 0) 

& - (1 . 0-RAM) **5* (RU 1**5/ 15 . 0**5+DELT/15 . 0) 

T2=DCOS (PSI 1) + (RU1+DELT) / 15 . 0 

EX120=T1/T2 

D6=CX6-3.0*C2**2 

E24=DX24+C2* (15. 0*C2**2-10 . 0*CX6) 

F120=EX120-15 . 0*C2*DX24+ 105 . 0*C2**2* (CXo-C2**2) -10 . 0*CX6**2 

CFF=C2/2 . 0 

DPF=D6/6.0 

EPF=E24/ 24 . 0 

FPF=F 120/ 120.0 

RETURN 

END 
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